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Abstract 

Propagating fronts are seen in varieties of non-equilibrium pattern forming systems 
in Physics, Chemistry and Biology. In the last two decades, many researchers have 
contributed to the understanding of the underlying dynamics of the propagating 
fronts. Of these, the deterministic and mean-field dynamics of the fronts were mostly 
understood in late 1980s and 1990s. On the other hand, although the earliest work on 
the effect of fluctuations on propagating fronts dates back to early 1980s, the subject 
of fluctuating fronts did not reach its adolescence until the mid 1990s. From there 
onwards the last few years witnessed a surge in activities in the effect of fluctuations 
on propagating fronts. Scores of papers have been written on this subject since then, 
contributing to a significant maturity of our understanding, and only recently a full 
picture of fluctuating fronts has started to emerge. This review is an attempt to 
collect all the works on fluctuating (propagating) fronts in a coherent and cogent 
manner in proper perspective. It is based on the idea of making our knowledge in 
this field available to a broader audience, and it is also expected to help to collect 
bits and pieces of loose thread-ends together for possible further investigation. 
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1 Front Propagation in Far from Equilibrium Systems 

1.1 Propagating Fronts into Unstable States: Deterministic Systems 

In pattern forming systems quite often situations occur where patches of dif- 
ferent bulk phases get separated by fronts or interfaces. In such cases, the rel- 
evant dynamics of the system is usually dominated by the dynamics of these 
fronts. When the interface separates two thermodynamically stable phases, as 
in crystal-melt interfacial growth problems, the width of the interfacial zone 
is usually of atomic dimensions. For such systems, one often has to resort to a 
moving boundary description, in which the boundary conditions at the inter- 
face are determined phenomenologically or by microscopic considerations. A 
question that naturally arises for such interfaces is the influence of stochastic 
fluctuations on the motion and scaling properties of such interfaces. 

4>t 




Fig. 1. A typical deterministic front propagating into a linearly unstable state cj) u 
and moving with speed v as : at the leading edge, the deviation of the front field 
<f>(x, t) from its unstable state value is infinitesimal and the nonlinear terms in the 
deterministic equation that describes the front can be neglected. For pulled fronts, 
the perturbation around the unstable state at the leading edge ahead of the front 
grows linearly and spreads with speed v* , thereby pulling the rest of the front along 
with it at = v* (hence the nomenclature). For pushed fronts, it is the bulk phase 
that drives the front propagation mechanism, as if the bulk phase (whose description 
requires nonlinear terms in the equation describing the front) pushes the leading 
edge from behind to achieve a front speed > v*. 

At the other extreme is a class of fronts that arise in pattern forming systems, 
where the occurrence of fronts or transition zones is fundamentally related to 
their inherent nonequilibrium nature, as they do not connect two thermody- 
namic equilibrium phases that are separated by a first order phase transition. 
In such cases — for example, chemical fronts [78], the temperature and den- 
sity transition zones in thermal plumes [121], the domain walls separating do- 
mains of different orientations in rotating Rayleigh-Benard convection [117], 
streamer fronts in discharges [39], or the aggregate fronts in diffusion limited 
aggregation [15] - - the fronts are relatively wide and are therefore described 
by the same continuum equations that describe nonequilibrium bulk patterns. 
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The lore in nonequilibrium pattern formation is that when the relevant length 
scales are large, (thermal) fluctuation effects are relatively small [26]. For this 
reason, in the first approximation, the dynamics of many pattern forming 
systems can be understood in terms of the deterministic dynamics of the ba- 
sic patterns and coherent structures. The first questions to study for them 
are properties like the existence and the asymptotic speed of propagation of 
the front solutions in deterministic equations. For example, in one dimension, 
the existence of an asymptotic propagation speed v as for deterministic fronts 
means that in the comoving co-ordinate x — t> as t, moving w.r.t. to the labora- 
tory co-ordinate a; at a speed t> as , the front profile approaches a fixed shape as 
t — > oo. Given the dynamical equations, the interest is in the magnitude of v as , 
and as well as in how the asymptotic front shape and speed are approached in 
time from a given initial configuration. In most cases of deterministic fronts, 
these equations are partial differential equations. 

Instances of such fronts, where fluctuation effects are small (and therefore 
they can be neglected in favour of a deterministic description) are abundant 
in physics [1,25,30,42,48,51,63,84,103,104,113,114], chemistry [49,50,75,120] 
and biology [17,41,58,82]. 1 As a result of detailed studies carried out in the 
last decade, it has emerged that in these systems, the dynamics is described 
by propagation of (deterministic) fronts into unstable states, i.e., the state of 
the system in the region far ahead of the front is linearly unstable [102]. These 
studies have classified deterministic fronts propagating into unstable states in 
two categories in a broad sense: the so-called pulled and pushed fronts (see Fig. 
1). Pulled fronts are the fronts that propagate into a linearly unstable state, 
and for which the asymptotic front speed t> as is the linear spreading speed v* 
of infinitesimal perturbations around the unstable state [12,29,38,100]. The 
name pulled fronts refer to the intuitive picture that at the leading edge of 
these fronts, 2 the perturbation around the unstable state ahead of the front 
grows and spreads with speed v* , while the rest of the front gets "pulled along" 
by the leading edge. On the other hand, fronts that propagate into a linearly 
unstable state and whose asymptotic speed is > v* are referred to as pushed 
fronts, as it is the nonlinear growth in the region behind the leading edge 
that pushes their front speeds to higher values. 3 If the state is not linearly 

1 On account of the fact that this article is about fluctuating fronts, these citations 
are representative and by no means complete. For a more comprehensive set of 
references, see Ref. [102], a review article on deterministic fronts. 

2 The term leading edge, which will appear in this review article over and over 
again, is (and will be) used to denote the front region where the value of the front 
field 4> is very close to its value at the the state it propagates into. In other words, 
in this region, the front evolution equation can be linearized around the value of <f> 
at the the state it propagates into. 

3 By definition, therefore, the asymptotic speed of a pushed front can be obtained 
only by solving the full nonlinear equation; in general it is not possible to do so 
except for special sets of parameter values. 
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unstable, then v* is trivially zero, and in addition, there can also be fronts 
propagating into unstable states with v* = 0; in such cases front propagation 
is always dominated by the nonlinear growth in the front region itself, and 
hence fronts in these cases are in a sense pushed too. 4 

For a front propagating into linearly unstable state, the linear spreading speed 
v* is obtained from the time evolution of a localized initial perturbation around 
the unstable state. A beautiful analysis of this can be found in Sec. 2.1 of 
Ref. [102]. The knowledge of the existence or the magnitude of v* alone how- 
ever does not answer the important questions like how and under what condi- 
tions one can expect a pulled front in the system, or when to expect a pushed 
front in the a given model. Addressing these questions satisfactorily is an in- 
volved process, and as already mentioned in footnote 4, from the point of view 
of this review article, in Sec. 1.1 we present the basic necessary results by con- 
sidering an example system. The equation that we choose for the illustration of 
the properties of pulled fronts is the so-called Fisher-Kolmogorov-Petrovsky- 
Piscounov equation (we will refer to it as Fisher-Kolmogorov equation here- 
after), which was at first used to model the spreading of advantageous genes 
in a population [43,59]. In this model, the density of the advantageous genes 
is denoted by (f>(x,t), a non-negative quantity, whose dynamics is described 
by the equation 



d(j) (J) 

— = 7—^- + — d) n . n > 1, for example 2 or 3 (1) 

at ox 2 

Equation (1) has two stationary states, of which <p u = is (linearly) unstable 
and <f) s = 1 is stable. Therefore, if the system is prepared in a way such that 
these two states coexist in a certain region of space, then the stable state 
invades the unstable one and propagates into it. The profile of the resulting 



4 At this point, I must warn the reader that I am being much more than just a little 
bit naive and simplistic in describing the broad classification of fronts propagating 
into unstable states in this manner. The actual issue of how v* comes out naturally 
for fronts propagating into unstable states and there onwards when one can expect 
a pulled front (possibly a further subclassification of coherent or incoherent pulled 
fronts) or a pushed front is a fairly involved subject in itself. For a proper analysis 
of it, the reader is encouraged to go through Sec. 2 of Ref. [102]. The point, once 
again, is that this is a review article on fluctuating (and propagating) fronts. From 
this point of view, it is enough to understand the basic issues that involve the 
classification in terms of pulled or pushed fronts in a broad sense. Such an attitude 
will be reflected all along Sec. 1.1, where I will highlight these basic features of pulled 
and pushed fronts in terms of examples and pictorial representations. Furthermore, 
in Sec. 1.1, I would completely leave out any discussion on deterministic pattern 
forming fronts such as Complex Ginzburg-Landau equation or Swift-Hohenberg 
equation etc., simply because literature on the effect of fluctuations in such pattern 
forming fronts do not exist in the literature. 
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front is similar to that of Fig. 1. 

To obtain the front solution admitted by Eq. (1), we rewrite it in a frame that 
moves w.r.t. the laboratory frame at a constant speed v, by means of a change 
of variables from x to the comoving co-ordinate £ = x — vt as 



~dt 



d(f) d 2 4> 



+ 



(2) 



and look for a stationary solution of in this comoving frame. The crucial 
relevance of the growth and spreading of infinitesimal perturbations enters 
naturally in this front solution, as the propagating infinitesimal perturbations 
around the unstable state in the leading edge ahead of the front sets on the 
instability making way for further growth. At the leading edge of the front, 
the 0- values are very close to the unstable state value, i.e., ^ <C 1, and one can 
neglect the nonlinear term <p n compared to in Eq. (2). The stationary solution 
of the resulting linear equation can then be solved by using (/>(£) ~ exp[— A£], 
yielding the relation 



v(X) = A + 



A 



(3) 



between v and A. The curve for the dispersion relation between v(X) and A is 



v(X) 



V 












Fig. 2. The dispersion relation (3) is schematically shown above, 
schematically shown in Fig. 2. It has a minimum at A*, and v* = v(X*) 



2. 



For any front propagating into a linearly unstable state, there is a solution 
of the leading edge of the form exp[— A£], and consequently, there exists a 
dispersion relation between A and the front speed v(X). Of course, the actual 
dispersion relation depends on the model that one studies. However, irrespec- 
tive of the model, the universality of pulled fronts lies in the fact that although 
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Fig. 2 indicates that one has a front solution for all values of A (and corre- 
spondingly all possible front speeds), from which it might a priori seem that 
the quantities A* and v* are not special in any way, the actual selection of 
the asymptotic front speed is obtained only after a proper stability analysis 
of the front profile in the comoving frame. Such a stability analysis yields the 
result that with an initial condition that (f>(x,t)\ t =o that decays faster than 
exp[— X*x] for x — > oo 5 , for long times, the front speed converges uniformly 6 
to v* as [13,38] 



v(t) =v*-^- t + 0(r^) , (4) 

as the front shape relaxes to its asymptotic configuration <j)*(x — v*t). How- 
ever, if the initial shape of the front is such that the leading edge is given by 
<f>(x,t)\t=o ~ exp[— Xx] for x — > oo with A < A*, then the front speed v re- 
mains fixed at v(X), given by the dispersion relation, while its shape remains 
unchanged at exp[— X(x — vt)]. 

(a) ^>^* (b) l< t 




x 



Fig. 3. An intuitive illustration of front speed selection: (a) for A > A*, and (b) for 
A < A*. The dashed and the solid profile propagate with their respective speeds. The 
kink (or the discontinuity in slope is only a symbolic representation of the crossover 
region [100]. A more precise analysis can be found in Sec. 2 of Ref. [102]. 

Although obtaining the power law (4) is mathematically quite involved, a 
beautiful qualitative illustration of the different speed selection behaviour of a 
pulled front, for A < A* and A > A* with an initial configuration of the leading 
edge <j)(x,t = 0) ~ exp[— Xx] as x — > oo, has been provided in Ref. [100] (a 
more precise analysis can be found in Sec. 2 of Ref. [102]). Imagine that we 
start with an initial configuration of the leading edge, shown in the leftmost 
curve of Fig. 3(a): an exponential decay exp[— Xx] with A > A* on the right 
denoted by the solid curve, and the exponential decay exp[— X*x] on the left 
denoted by the dashed curve. The two curves are joined together at a kink 
shown by a small filled circle. One should keep in mind that a kink in the initial 
front profile is not propagated as it is for fronts that involve higher than first 

5 this condition is also known as "sufficiently steep initial condition". 

6 Uniform convergence means that the convergence behaviour (4) of the front speed 
is the same irrespective of the value of <f> at which the speed is being measured. 
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order spatial derivative, but nevertheless, one can get a clear intuitive picture 
of the evolution of the leading edge by evolving the solid curve exp[— X*x] 
and dashed curve exp[— Xx] profile separately and by following the position of 
the kink in a pictorial representation. The U-shape of the dispersion relation 
v(X) vs. A dictates that the speed of dashed curve profile v(X) is higher than 
the solid curve profile v*, as shown in three later snapshots of the evolving 
leading edge, taken at times ti < £3 < £4. As a result of the differences in 
speeds of the dashed and the solid curves, the height of the kink from the 
a;-axis keeps decreasing, which shows that front profile at the leading edge is 
being taken over by the dashed curve, corresponding to an exponential decay 
with exponent A*. A similar picture is shown in Fig. 3(b), where the solid 
and dashed curves correspond to exp[— Xx] and exp[— X*x] respectively, and 
A < A*. The situation is reversed in this case — the front profile at the leading 
edge is being taken over by the solid curve profile. 7 

The discussions below Eq. (3), so far, is relevant only with an underlying 
understanding that we are considering a front (propagating into a linearly 
unstable state) that asymptotically yields a pulled front. In that sense, the 
usage of Fisher-Kolmogorov equation as an example is very well-placed, since 
it asymptotically admits only a pulled front solution. As for the question 
whether an equation describing a propagating front into an unstable state 
gives rise to a pulled front or not, a necessary condition is not known, but 
for a sufficient condition, it is known that if all the nonlinear terms in the 
time evolution of the front suppress growth for a front that propagates into a 
linearly unstable state, then the resulting front is a pulled front. In particular, 
for the nonlinear diffusion equation d t (f) = Dd^cf) + /(</>) (an equation that 
has been extensively studied in the front propagation literature), the above 
sentence implies that one can expect pulled fronts if f{u)/u < f'(0) [102]. 

A similar general analysis [i.e., Eqs. (2-3) and related discussions] for the 
asymptotic speed selection mechanism for pushed fronts does not exist. This 
is not a surprise in itself, as for pushed fronts, the front speed is really deter- 
mined in the nonlinear bulk phase of the front (which one cannot solve except 
for special sets of parameter values). As a result, the leading edge of the front 
does not play any role in the front speed selection. Nevertheless, it is known 
that for sufficiently steep initial conditions, as opposed to the \jt convergence 
to the asymptotic front speed v* for pulled fronts, the convergence to asymp- 

7 The discussion in the above two paragraphs may, at first glance, seem to downplay 
any special significance of v* — after all, any other front speed > v* can be reached 
with an appropriately chosen A in the initial configuration of the front. However, 
in most physical situations, one is interested in the dynamics of "localized initial 
conditions", for which the spatial decay of the leading edge is steeper than that of 
exp[— X*x]. From this point of view, the exp[— Ax] (with A < A*) initial configuration 
of the leading edge, which leads to the eventual front speed v(X), is very special, 
and therefore is not interesting. 
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totic front speed for pushed fronts 8 is exponential in time. For the front 
speed, a more mathematical representation of the exponential convergence for 
pushed fronts vs. the non-exponential convergence for pulled fronts can be 
traced to the stability criteria of the asymptotic front solutions [37,101,102]. 
The idea is that at long times, an intermediate front profile evolving towards 
its final asymptotics can be decomposed as the asymptotic front profile and 
infinitesimal localized perturbations around it. The convergence properties to 
the asymptotic front speed is then determined by how fast, in the comoving 
frame of the asymptotic front profile, these infinitesimal perturbations decay 
in time. 



To elucidate the relation between the stability of the asymptotic pulled or 
pushed front solutions and their convergence (in time) behaviour, let us con- 
sider a front solution propagating with a speed v, and decompose an 
intermediate front profile (f>(x, t) as 



<KS,t) = $,(0 + v(t,t), (5) 

such that r/(£, t) is of infinitesimal magnitude everywhere. At the leading order, 
the dynamics of r)(£,t) is described by 9 



^ = C vV + 0(rj 2 ) : (6) 

where the linear operator C v is obtained by linearizing the equation for t) 
in the comoving frame of the asymptotic solution. 10 In Eq. (6), r)(£,t) can 
be expanded as a linear combination of the eigenf unctions of the operator of 
C v . 11 In this expansion, there is an eigenfunction with zero eigenvalue that 
corresponds to the translational invariance of the asymptotic front solution, 
while the nature of its lowest-lying eigenvalues decides the decay properties of 
time convergence properties of the asymptotic front speed [37,101]. For pushed 
fronts with asymptotic front speed v\ the eigenspectrum of L v \ for these 
lowest-lying eigenf unctions is gapped — this indicates that the convergence 
to the asymptotic front solution <${x — vH) is exponential in time. On the 
other hand, for pulled fronts, the spectrum of C v * is gapless, indicating that 
the convergence to the corresponding asymptotic front solution <p*(x — v*t) 



8 Just like v* is a standard notation for the linear spreading speed, v* is also a 
standard notation for the asymptotic speed of pushed fronts. 

9 In this example, we use an equation that involves only a first order derivative in 
time. A similar analysis can also be worked out for higher time derivatives. 

10 E.g., for Eq. (1), C v = Ddf + vd^ + (1 - 

11 In general, the operator C v is not a Hermitian operator, and therefore its left and 
right eigenvectors are different. 
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has to be non-exponential in nature [101]. 



1.2 Fluctuating (and propagating) Fronts: A Separate Paradigm 

From the point of view that the thermal fluctuations are rather small for the 
pattern forming systems that relate to the inherent non-equilibrium nature 
of the bulk phases, it is less of a surprise that the effect of fluctuations on 
propagating fronts has not attracted the attention of physicists until relatively 
recently. The understanding that the fluctuations are indeed important was 
not only motivated by the realization that matter is composed of discrete 
particles, but also from the fact that there are situations, where fronts are 
naturally made of discrete constituents on a lattice (see for example, Refs. 
[16,28,55,62,96,118,119,124]). In addition, it is also of general interest to see 
the effect of externally added noise on the otherwise deterministic fronts, to 
see how severely this properties are affected, or if these noise terms give rise 
to new phenomena. Although in some cases, studies concerning the effect of 
externally added fluctuations on the fronts were motivated by such theoretical 
interests, e.g. in Refs. [3,4,79,94,97,99,105,107,108,115], in some others, 
the noise terms were added to model the effect of discreteness of constituent 
particles and the lattice, e.g. in Refs. [35,36,67,70,95]. In this review article, 
we will review both of these cases. In Sec. 2, our focus will be on propagating 
fronts made of discrete particles on a lattice, there we will review all the 
known results of such fronts. In Sec. 3, we will review propagating fronts 
where fluctuations are introduced by means of externally added noise terms. 
The discussion regarding to what extent the externally added fluctuations 
correspond to discrete particle models will be discussed subsequently. 

The presence of fluctuations, be it as a result of discreteness of particles and 
the lattice, or be it as a result of the externally noise terms to otherwise 
deterministic equations, immediately implies that (i) at a contrast to the de- 
terministic fronts propagating with a fixed shape at a fixed speed at t — > oo, 
strictly speaking, different snapshots of one particular realization of a fluctu- 
ating front, taken at different times, will be microscopically different from each 
other; although in an average sense, each of these realizations will have the 
same well-defined shape. The consequences of this is the following: unlike the 
deterministic fronts, the front position defined by the location x of the point, 
where the number of particles per correlation volume is a fixed number, say Nq 
(or by the point where the front field reaches a fixed value 0o), does not move 

12 The continuous spectrum of L v * for pulled fronts is known to be responsible for 
the breakdown of the so-called solvability analysis; for a detailed study, see Ref. [37]. 
Moreover, the continuity of the spectrum of C v * is a mathematical representation 
of the so-called marginal stability criterion of pulled fronts, see for example Refs. 
[12,22,23,37,92,101,102] etc. and references cited therein. 
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with a constant speed for a given realization of a fluctuating front even for 
large times. It is only when the movement of this point over a long period of 
time is considered (at large times), the speed of this point, calculated as a long 
time average, can then be defined as the asymptotic front speed. 13 (ii) More- 
over, just like any stochastic quantity having a variance around its average, 
if one follows Xo as a function of time, then one also expects ([xo(t) — f as £] 2 ) 
to increase over time. The angular brackets denote an averaging over an en- 
semble, each individual member realization of which has reached their steady 
shapes in an average sense. This means that the displacements of individual 
front realizations w.r.t. each other as a whole front keeps increasing with time 
- this a phenomenon known as the wandering of fluctuating fronts. 

10 4 





20 40 

X 

Fig. 4. Adopted from Ref. [16], this figure depicts the front shapes for the reac- 
tion-diffusion process X ^ 2X for discrete particles on a lattice, indexed by A. The 
rate of forward reaction X — > 2X is N/2 with N = 10 4 , while rate for the backward 
reaction 2X — > X has been normalized to unity. Particles diffuse to their nearest 
neighbour sites with a unit rate. Solid curves: front profiles for seven different re- 
alizations at a certain time. Solid curve with error bars: front profile obtained by 
averaging over number of particles on each lattice site for 66 realizations. 

An example of the phenomena in (i) and (ii) of the previous paragraph for 
fluctuating fronts is shown in Fig. 4. It has been adopted from Ref. [16], 
and it depicts the shapes of fronts made of discrete particles in a reaction- 
diffusion process X ^ 2X on a lattice indexed by A. The rate of forward 
reaction X — > 2X is N/2 with N = 10 4 , while the rate for the backward 
reaction 2X — > X has been normalized to unity. The seven solid curves are 
the front shapes obtained by taking simultaneous snapshots of seven different 
realizations at a certain time. All these realizations started with the same 
initial microscopic configuration. At a microscopic level, the shapes of these 
snapshots are different, but in an average sense, each of them has a shape 

13 Clearly, for the asymptotic front speed to be well-defined, it should be independent 
of Nq (or <fi ), which in itself is a consequence that in an average sense, the front 
realization has a well-defined shape. 
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identical to that of the dashed curve. The displacements of the individual 
front realizations as a whole w.r.t. each other demonstrate front wandering. 
On the other hand, the solid line with error bars in it depicts the front shape 
obtained by averaging over number of particles on the lattice sites for 66 
different front realizations. As one can expect, the front obtained by such an 
averaging process is much wider than any of the seven realizations, since such 
an averaging fails to filter out the effect of front wandering — in future, we 
will refer to this averaging process as "simple averaging" . 

The above discussion, therefore, points us to the following conclusions: (a) to 
obtain the average asymptotic shape of fluctuating fronts, one has to separate 
out the wandering effects, i.e., one has to first realign the simultaneous real- 
izations and then calculate the average shape — from now on, we will refer to 
such an averaging as "conditional (ensemble) averaging", 14 and ((Nk)) will 
denote the conditionally averaged number of particles in the k-th lattice site, 
(b) For measuring the front speed, a conditional averaging is not necessary; 
one can obtain the front speed by tracking how fast a given value of the sim- 
ply averaged front profile 0, say <f> , moves. However, we have seen before, for 
deterministic fronts, that the front shape and the speed are related to each 
other, and therefore, from that angle, for a theoretical prediction/analysis of 
front speed, one does need the conditionally averaged front profile. A good 
description of fluctuating fronts must take both (a) and (b) into account. In 
addition, (c) a comprehensive description of fluctuating fronts must not only 
include the expression of the asymptotic front speed and as well as its approach 
in time towards a steady (conditionally averaged) front shape starting from an 
initial configuration, but also must it identify the nature of its wandering (i.e., 
whether it is diffusive, sub- or super- diffusive), along with the relevant charac- 
terizing exponents. 15 Such a comprehensive description for fluctuating fronts 
propagating into unstable states is the central theme of this review article. In 
one spatial dimension, our convention will be to consider fronts propagating 
from the left to the right, and keeping in line with the propagation direction 
of the front, we will often use forward for rightward and backward for leftward 
interchangeably. 

Although I claim to provide a comprehensive description for fluctuating fronts 
in this review article, its precise meaning has to be properly laid out at this 
point — after all, a front is nothing but an interface that demarcates the 
boundary between two different phases within a system; and the properties of 

14 It is precisely this conditionally averaged front profile that is described by an ap- 
propriate deterministic mean-field equation; in future, these conditionally averaged 
front profile will be denoted by (f>^ . To find out how to align different snapshots of 
fluctuating fronts in the sense of (a) above, see Eq. (78) and footnote 65. For now, 
the mere assumption that (f>^ exists will be enough. 

15 In one spatial dimension, fluctuating fronts are always diffusive; in higher dimen- 
sions, the situation is more complicated. 
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such interfaces for, e.g., fluid flows in porous media, flame fronts in a burn- 
ing paper, atom deposition processes in molecular beam epitaxy experiments, 
aggregate fronts in diffusion-limited-aggregation etc. have already been very 
well-studied in the literature. Analyzing the properties of fronts in such vastly 
different systems is by no means the intention of this review. Instead, my 
purpose here is the following: in spite of the existence of a large variety of dy- 
namical equations that describe moving interfaces in different systems, there 
is an overall underlying framework for the propagation mechanism of deter- 
ministic fronts; thus, I want to review the effect of fluctuations on propagating 
fronts under the same framework in a comprehensive manner. 

In general, providing such a comprehensive description for fluctuating fronts 
is by no means an easy task. As mentioned in the first paragraph of Sec. 1.2, 
two major approaches have been used to this end — the first one being that of 
noisy dynamical equations, where one simply adds external noise terms to the 
otherwise deterministic equations to introduce fluctuations. In the second ap- 
proach, many researchers in the last few years have taken the route of studying 
fluctuating fronts that are constituted of discrete particles on a lattice, where 
stochasticity of the microscopic dynamics of the particles themselves gives 
rise to fluctuations. In a way, the philosophy of both these approaches is to 
first have a deterministic front as the underlying structure, and then to study 
the effect of fluctuations superimposed thereupon: there exists a parameter 
in these approaches that controls the fluctuation strength. For example, in 
the second approach, this parameter is 1/N, where N is the conditionally av- 
eraged number of particles per lattice site in the stable phase of the front, 
and the deterministic mean- field results are obtained when N — > oo. In the 
recent past, there has also appeared a new approach (along the lines of quite 
the opposite philosophy of the two approaches discussed above) to extract the 
mean-field behaviour at small N [81]. Therein, one defines a variable N* that 
measures the number of particles per correlation volume in the stable phase 
of the front, and depending on the system parameters, one can have N* — > oo 
for correlation volume — > oo and small N; in this limit, one also recovers the 
deterministic mean- field description for the front. 

Nonetheless, it has to be noted that the majority of the studies on fluctuating 
fronts have been carried out on fluctuating "pulled" 16 fronts, i.e., on fronts 
whose deterministic limit yield pulled fronts, as opposed to on fluctuating 
pushed fronts 17 . This is not very surprising — first of all, for the determin- 
istic counterparts, it is the pulled fronts, for which there is a solid theoretical 



16 The use of the quotes is motivated by the fact that fluctuating fronts, whose 
deterministic limit yield pulled fronts, are actually weakly pushed. See points IV 
and V of Sec. 2.1 for more details. 

17 In view of footnote 16, what we mean by fluctuating pushed fronts is that the 
deterministic limit of these fronts yield pushed fronts. 
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understanding. Secondly (and more importantly), it is the sensitivity of the 
pulled fronts to the dynamics of the leading edge, which is severely affected 
by the discreteness of particles and the lattice, or by the presence of external 
noise terms — as a result, fluctuating "pulled" fronts exhibit many surprising 
characteristics. 
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1.3 A Note on the Notations 



Before we proceed further on, it will be very helpful to have the convention 
of the notations used summarized. Below we provide a table for this purpose. 
Any notation that does not appear in Table I is explained before its usage in 
this review article. 



x, k 


Notations for continuum space and for lattice index respectively. 


V* 


Notation for the linear spreading speed for deterministic fronts 
propagating into unstable states, and also for the selected asymp- 
totic front speed for (deterministic) pulled fronts with a steep 
enough initial condition. Has already been introduced in Sec. 1.1. 
Applies to deterministic pulled fronts both in continuum space and 
on a lattice. 


«t 


Notation for the selected asymptotic front speed for (determin- 
istic) pushed fronts. Also introduced already in Sec. 1.1. Applies 
to deterministic pushed fronts both in continuum space and on a 
lattice. 


Was 


General notation for the asymptotic front speed of any front 
(pulled and pushed) in any system (deterministic and stochastic; 
discrete and as well as continuum space). 


V N 


Notation for front speed for a discrete particle and lattice system 
of fronts, where there are N particles on (conditional) average per 
lattice site. 


c 

? 


Notation for comoving co-ordinate for pulled (£ = x — v*t, or 
£ = k — v*t), or fluctuating "pulled" fronts (£ = k — v^t), or even 
for fronts whose deterministic mean-field limit yield pulled fronts 
but whose speeds have not been determined yet (in these cases, £ 
will be used to denote x — v^t or k — v^t or k — v^t as appropriate). 


c 


Notation for comoving co-ordinate for pushed (£ = x — vH, or 
( = k — v't), or fluctuating pushed fronts (( = k — vjyt), or even 
for fronts whose deterministic mean-field limit yield pushed fronts 
but whose speeds have not been determined yet (in these cases, ( 
will be used to denote x — v^t or k — v^t or k — v^t as appropriate). 




Notation for the front profile for pulled fronts; either £ — x — v*t 
or £ = k — v*t. Already introduced in Sec. 1.1. 


f (C) 


Notation for the front profile for pulled fronts; either ( = x — vH 
or £ = k — vH. Already introduced in Sec. 1.1. 



Table I: Convention behind the notations used in this review article. 
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2 Fluctuating Fronts in Discrete Particle and Lattice Systems 



2. 1 Summary of Known Results 



In Sec. 2.1, we summarize all the known results for fluctuating (propagating) 
fronts, made of discrete particles on a lattice. Among the physical systems, 
where front propagation has been studied under this scheme, various sorts 
of reaction-diffusion systems constitute the majority [5, 14, 16, 34, 54-56, 66, 
70,71,73,80,86,87,96,116,118,119]. In addition, propagating fronts were also 
studied in the context of directed polymers in random media [19-21,24,27,28], 
and the calculation of the largest Lyapunov exponent in a gas of hard spheres 
[122-126]. In all these models, fronts propagate from a stable to an unstable 
state. The unstable state is constituted of empty lattice sites, whereas in our 
notation, the (conditionally) averaged number of particles on the A;-th lattice 
site ((N k )) in the stable state and the corresponding asymptotic front speed 
are respectively denoted by N from now on. Also, unless otherwise stated, in 
Sec. 2, fronts will be assumed to propagate in one spatial dimension; fronts in 
higher than one spatial dimensions will be considered only in Sec. 2.6.6. 

For these models, the following results have been obtained in the last decade: 

I. The limit of deterministic pulled front propagating with speed v* is 
reached for N — > oo. However, the convergence of to its determinis- 
tic limit v* for increasing value of N is extremely slow. The approach of 
vn to v* is from below, and for asymptotically large N values, the leading 
order form of v N is model independent [19-21,56,86,87,122,124-126], given 
by 



* d 2 v(\) 
Vn — v ■ 



d\ 2 



2 \*2 



7l 2 X 



ln 2 N 



(7) 



II. For intermediately large values of N, the subdominant corrections to the 
leading order result (7) are very strong, and they depend on the details of 
the model. For some models, these corrections can be so strong that one 
does not observe the asymptotic scaling (7) unless N is extremely large 
[16,56,66,70,86,87]. 

III. Unlike the semi-infinite leading edge for deterministic pulled fronts, in 
any snapshot of these discrete particle and lattice systems for any value 
of iV at any time, there exists a "foremost occupied lattice site (f.o.l.s.)", 
on the right of which no lattice site has ever been occupied before. In the 
hopes of providing a theory that would yield front speeds at intermediately 
large values of N for these fronts, a stochastic dynamics of the f.o.l.s., and a 
(standard) deterministic and uniformly translating front solution reasonably 
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far behind has been proposed in Ref. [86,87]. The two solutions are then 
matched in the "tip region" of the front, spatially extended over a few lattice 
sites behind the f.o.l.s. (and also including the f.o.l.s.). Despite the fact that 
this formalism yields a consistent (and verifiable) picture of these fronts, 
it is not predictive, and a first-principle based predictive theory for vn at 
intermediately large values of N is still lacking. 

IV. Unlike pulled fronts, where the asymptotic front speed v* is reached only 
if the initial front profile decays faster than exp[— \*x] for x — > oo, it can be 
argued that the conditionally averaged front profile for fluctuating "pulled" 
fronts is reached uniquely, independent of the initial microscopic configura- 
tion. This observation has been made in Ref. [21], and it is consistent with 
all the known simulation results so far. 

V. In the limit of N — > oo, the spectrum of the stability operator C VN 
[c.f. Eq. (6)] is gapped, and for the lowest-lying eigenvalues, the gap is 
oc 1/ In 2 N [56,89]. Although we will get into more details later on, here we 
briefly mention that points IV and V together bear the signature of these 
fronts being weakly pushed. 

VI. In the limit of N — > oo, the diffusion coefficient of the front Df that char- 
acterizes front wandering, approaches zero. For fluctuating "pulled" front 
in a model that closely resembles the so-called clock model [124], Df has 
been shown to scale as 1/ In 3 N in the large N asymptotic regime [21]. 

VII. As N is reduced "too much" , the fluctuations in the number of particles 
per lattice site in the stable state become stronger and stronger. Perhaps in 
view of II above, in the sense that due to the lack of first-principle based 
predictive theory for finite-iV corrections, works on such values of N are 
somewhat rare [71, 118, 119]. The extreme limit of this, namely the case 
when there is at most one particle is allowed per lattice site, however, is 
very well studied (we will refer to this scenario as "N < 1"). In view of the 
fact that there are simply too many results to summarize here, the readers 
are referred directly to Sec. 2.6 for further details. 

VIII. As far as the discrete particle (and lattice) systems of fluctuating 
pushed fronts are concerned (c.f. footnote 17), at the time of writing this 
review article, I have not seen any work on the behaviour of v n for N — > oo. 
Nevertheless, having combined the knowledge that one musters from points 
I through VI above and the solvability analysis for pushed fronts [37], for 
N — > oo, the convergence of the asymptotic front speed to (c.f. footnote 
8) for fluctuating pushed fronts to their deterministic limit behaves as a 
power law of N. 

In the rest of Sec. 2, we will elaborate points I- VIII, and provide a unified 
picture of discrete particle systems of fluctuating fronts. Details of which point 
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is elaborated where are given below: point I in Sec. 2.2, points II and III in 
Sec. 2.3, points IV and V in Sec. 2.4, point VI in Sec. 2.5, point VII in Sec. 
2.6, and point VIII in Sec. 2.7. Finally, in Sec. 2.8 we put the results of Sees. 
2.2, 2.4-2.5 and 2.7 in a unified perspective. 

2.2 Derivation of 1/ In 2 N Convergence of the Asymptotic Front Speed to v* 
for Fluctuating "Pulled" Fronts as N — > oo 

From its definition in Sec. 1.1, pulled fronts can be clearly seen to be realized in 
a continuum description, for otherwise the "infinitesimal perturbation around 
the unstable state" does not quite make sense. Through its definition (and 
propagation mechanism), pulled fronts and infinitesimal perturbations around 
the unstable state are intimately intertwined — ahead of the leading edge of 
a pulled front, the perturbation around the unstable state is infinitesimal, as 
it grows and propagates itself with speed v* and pulls the rest of the front 
along with it. From this perspective, the discrete particle and lattice systems 
of fluctuating "pulled" fronts belong to a different category — at the leading 
edge of these fronts, due to discreteness of the particles, the smallest amount of 
the "infinitesimal perturbation" for growth to start on any lattice site cannot 
be less than that of one quantum of particle, a fact that is more than likely 
to yield a front speed different from v*. 

In hindsight, the above argument sounds simple, and it is appreciated easily 
when the conditionally averaged number of particles per lattice site is trans- 
lated into the language of of Sec. 1.1. Notice that in the stable phase of the 
front, {{N k )) = N corresponds to = 1. As a result, the lowest amount of 
corresponding to one quantum of particle, needed on any lattice site ahead of 
the leading edge for growth to begin is 1/JV. The pulled front limit in these 
discrete particle systems is therefore expected to be reached only for N — > oo. 

These subtleties were first realized by Brunet and Derrida [19], shortly after 
simulations revealed that for models of fluctuating "pulled" fronts the devi- 
ations of the asymptotic front speed from v* can be significant [16], even for 
N ^> l. 18 To mimic the lowest quantum of particle that is needed on any 
lattice site ahead of the leading edge for the growth to begin, Brunet and 
Derrida conjectured that the dynamics of a fluctuating "pulled" front at large 
N should be described by an equation with a growth cutoff for <f> < e, where 
e = l/N. 

To verify their ideas, they used modified Eq. (1): 

18 Severe deviations of the front speeds from v* for N not so larger than 1 have also 
been observed in simulations, but we leave those aside until we discuss point VII at 
length in Sec. 2.6. 
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The stationary front solution of Eq. (8), propagating with speed v e , was then 
divided into three regions: I, II and III. In region II and III, <fi <C 1 and hence 
the nonlinear term 0™ can be ignored; in addition, in region III, <fi < e, which 
leads to a further simplification of Eq. (8). On the other hand, in region I, one 
has to consider the full nonlinear equation. This leads one to 
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At the boundary between regions II and III, the continuity of the value of the 
field as well as of its comoving derivative </>'(£) has to be maintained. With 
the understanding that £ < 1, the leading order solutions that obey these 
boundary conditions were shown to be [19] 19 
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in region III, and 
e -A*(£-£i) in re gi on n, (10) 



such that the boundary between regions II and III lies at £o — |hie|/A*, 
whereas the solution in region II crosses over to that in region I at £i, and A 
is a quantity of 0(1); yielding 
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d 2 v(X) 
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(11) 



Notice that to obtain the expression (11) of v e , the solution of Eq. (9) in region 
I is not required. Moreover, although from solution (10) might seem to 

go to zero at £ = violating the monotonicity of the front profile, in reality, 
there exists no sharp boundary between regions I and II, like the one at £o- 
Instead, the solution (10) crosses over to that of region I, and the crossover 
happens in a way that maintains the monotonicity of the front profile. At the 
leading order, the characteristic position £ = £i for the crossover is given by 

6-0. 



19 Notice the notation <p^: the idea is motivated by the fact that as N — ► oo, the 
conditionally averaged front shape is given by the effective deterministic mean-field 
equation (9). See footnote 14 in this context. 
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Few points in Brunet and Derrida's conjecture must be noted: (a) One arrives 
at the leading order large N asymptotic behaviour of the front speed (7) by 
setting e = 1/N. The form of is independent of the details of the micro- 
scopic dynamics of the constituent particles, although the actual dispersion 
relation v(X) vs. A is obtained from the equation obeyed by the deterministic 
mean-field limit of the microscopic model at hand, (b) Region III is charac- 
terized by an effective (finite) width l/v £ due to the exponential decay of the 
front profile. On the other hand, with the substitution e = 1/N for asymp- 
totically large N, the length (£ — £x) of region II increases logarithmically 
with N, while the length of region I remains finite — this yields an effective 
front width, scaling as IniV. (c) Furthermore, the leading order form of Vn is 
independent of whether one substitutes e = 1/N or e = c/N, with c of 0(1). 
In addition to amplitude A, c simply provides a subdominant correction to 
the r.h.s. of Eq. (7), if c ^ 1. 
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Fig. 5. The plot of v* — vn vs. N for the discrete particle and lattice model of 
Refs. [19-21] relating to the study of directed polymers in random media. In the 
above figure, v m i n should be read as v*. The solid curves represent Eq. (7) for various 
p values [20]. A corresponding graph for the same model for much higher values of 
N for p = 0.25 can be found in Fig. 2 of Ref. [21]. 

In the history of this subject, Brunet and Derrida's insight has turned out 
to have had an enormously significant impact. Their prediction of the N- 
dependence of the front speed was immediately confirmed in a fluctuating 
"pulled" front model (where time was also discretized) [19-21] 20 relating to 
directed polymers in random media (see Fig. 5) [24,27]. Subsequently, in an- 
other model of a fluctuating "pulled" front - - namely, the so-called clock 
model relating to the calculation of the largest Lyapunov exponent in a gas of 
hard spheres [122, 124-126] 21 - the idea of a cutoff was employed; and the 

20 The deterministic mean-field limit for this model is described by the equation 
<j>(k,t + l) = l-p[l- <p(k - l,t)] 2 - (1 -p)[l - 4>(k,t)} 2 , where k is the lattice site 
index and p is a fixed number in the interval (0,0.5). The corresponding U-shaped 
dispersion relation v(X) vs. A has the form v(\) = In [2pe x + 2(1 — p)] /A. 

21 The deterministic mean-field limit of the so-called clock model is described by 
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sinusoidal profile of Eq. (10) at the leading edge, along with the result that 
v * — vn oc 1/ln 2 N for asymptotically large N, was obtained (see Fig. 6). 




N 



Fig. 6. The plot of v* - v N vs. N for the so-called clock model [122, 124-126]. 
Symbols: simulation data; solid curve in the main graph: expression (7); inset: the 
same data as of the main graph plotted differently, with 5v = v* — vjy. 



Although a continuum-space equation was used in Ref. [19] for the deviation 
of Eq. (7), it turns out that for asymptotically large N, one can motivate the 
sinusoidal profile of the leading edge along the lines of the discussion in the 
first two paragraphs of Sec. 2.2 and the positivity of (in these models 4> k 
represents the scaled number of particles on the lattice site k, and hence it 
can never be negative) [86,87, 122, 125, 126], from which one arrives at Eq. 
(7) very simply for any given microscopic model. The idea is the following: 
since the growth of particles ahead of the leading edge of these fronts does 
not start unless there is at least one particle on any lattice site, one can 
expect the discreteness of particles to inhibit the spreading of the leading 
edge, resulting in a front speed Vn smaller than v*. However, v* is the minimum 
value of v (A) for any real A in the deterministic mean- field description of the 
microscopic model at hand, and therefore the only way to obtain vn < v* 
from the dispersion relation v{\) vs. A is to consider Im(A) ^ 0. 

Equation (10) originates from the consideration of non-zero Im(A), which en- 
ters the solution as the inverse of the wavelength of the sin-function. The 
positivity of <fi then restricts the value of its argument to lie within (0, n) 
(modulo a constant phase factor, which can be set to zero by shifting the 
position of the origin in the comoving frame). For iV — > oo, Re(A) — > A*, 
Im(A) — > and — > v*, and as a result, for asymptotically large N, one can 
expand the front solution around these asymptotic values. Such an expansion 

-q— = —<Pk + leading to the dispersion relation v(X) = [2e x — l] /A. 
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confines the values of Re(A), Im(A) and vn close to the minimum of the dis- 
persion relation v (A) vs. A on the Re(A)-plane, and at the leading order, it 
yields the following results: (a) Im(A) = In A/7rA*, just like in the argument of 
the sin-function in Eq. (10), and (b) the large N asymptotics (7), where the 
second derivative v (A)| A * appears due to the fact that near the extremum of 
v(X) on the Re(A)-plane, the variation of the front speed for a small variation 
of Re(A) from A* is oc t>"(A)| A ». 

We conclude Sec. 2.2 with a short note on the subdominant corrections to the 
oc l/\n 2 N deviation of v N from v*. As mentioned earlier, according to the 
derivations in Ref. [19], along with other explicit dependence on N, both A 
and c enter the expression of the subdominant corrections. For asymptotically 
large values of N, the behaviour of the subdominant corrections is dominated 
by these explicit A-dependent terms rather than the ones involving A and 
c. For the microscopic model relating to directed polymers in random media, 
a 1/ In 13 N form has been suggested for the leading order behaviour of the 
subdominant corrections in Ref. [21], where (3 has been estimated to lie some- 
where between 2.5 and 3. However, for the reaction-diffusion system 2X in 
Ref. [86,87], at the leading order, the subdominant corrections seem to behave 
as In [In iV]/ln 3 N. Although this issue has not yet been settled, it is only a 
theoretical and a marginal one — after all, at A values when the subdominant 
corrections start to become stronger, a very different description of fluctuating 
"pulled" fronts is called for, as we will find out next. 



2.3 The Case of Intermediately Large Values of N, Foremost Occupied Lat- 
tice Site, Tip Region of the Front and All That 

The large A asymptotics (7) describes v N as an explicit function of A at the 
leading order. However, from comparison with the simulation data, it seems 
that to observe this asymptotic result, one often has to go to very high values 
of A. In support of their theoretical prediction (7) in Ref. [19], Brunet and 
Derrida took A up to 10 16 (Fig. 5) [19, 20] in the lattice model relating to 
directed polymers in random media; and in a later work [21], by means of a 
clever simulation algorithm, they further took A to as high as 10 150 in the 
same model. In these studies, the range where A is intermediately large but 
not asymptotically large, was not examined carefully. Amid this backdrop, a 
reinterpretation of the numerical results of Ref. [16] for the reaction-diffusion 
system X±=? 2X suggests that there are strong deviations of the data from the 
prediction (7) when A is of (9(10 4 -10 6 ). Such deviations were also noted by 
Kessler and co-workers [56] — their simulations for fluctuating "pulled" fronts 
in the reaction-diffusion system X+Y^ 2X on a lattice suggested that even 
for A as large as 10 12 , differences between simulation data and the large A 
asymptotics (7) remain significant (see Fig. 7). In view of these, one's worry 
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naturally shifts to the subdominant corrections to the r.h.s. of Eq. (7). 
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Fig. 7. Simulation data from Ref. [56] for the reaction-diffusion system X+Y^ 2X. 
The rate of diffusion D = 0.5, the reaction rate is 0.1/N, and Sv = v* — vjy. Note 
the increasing trend of the data with increasing N. For these parameter values, 
according to Eq. (7), 5v In 2 N/v* should be ~ —9.6. The figure is an altered version 
of Fig. 5 in Ref. [56] . The alteration was carried out for greater clarity. 



The common denominator of all the derivations of Eq. (7) [19,86,87,122,125, 
126] for fluctuating "pulled" fronts is that they were carried out to describe 
a conditionally averaged uniformly translating profile all over the front. In 
such a description, for the large N asymptotics, the leading order expression 
of v jy is independent of specific details of the microscopic model. The specific 
model dependence 22 of v n enters in the subdominant corrections through two 
parameters, A and c (e — c/N). In that sense, Brunet and Derrida's derivation 
of Eq. (7) is quite instructive - - it shows that the speed of a fluctuating 
"pulled" front is actually fully determined from the overall shape of the front 
(thereby taking into account the specific details of the model in a conditionally 
averaged manner), and not only from the property of the leading edge — after 
all, to obtain the value of A, one needs to know the overall front shape, even 
in the region where the nonlinearities are important. 23 On the other hand, 
although c was set equal to 1 by Brunet and Derrida, one might wonder what 
its true value is, and to what extent it influences the subdominant corrections 



22 The dependence of vn on the specific details of the model is apparent from a 
comparison of Figs. 6 and 7. While in Fig. 7, even at N as large as 10 12 , the 
l/ln 2 iV scaling of v* — vn seems rather far off for the reaction-diffusion system 
X+Y^ 2X, Fig. 6 finds confirmation of the 1/ln 2 N scaling at N ~ 10 3 for the 
so-called clock model. 

23 This stands in sharp contrast to pulled fronts. To obtain the front speed v* for 
pulled fronts, one only needs the dispersion relation v(X) vs. A. The dispersion 
relation is obtained from the consideration of the leading edge of the front, where 
the nonlinearities are not important. 
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for a given value of N in a microscopic model. 24 It is in fact conceivable that 
these parameters have further dependence on N; nevertheless, the feasibility 
of ever obtaining an explicit TV-dependent prediction of for any value of 
iV [56] remains an open question. 

In view of these difficulties, if one were to provide a description of fluctuating 
"pulled" fronts to yield front speeds vn for moderately large values of N as 
well, one must certainly look elsewhere to address a pressing concern; namely 
the effect of stochasticity in the microscopic dynamics of the individual parti- 
cles on the propagating front, an effect that disappears when the conditional 
averaging of a large number of front realizations is performed all over the front 
region. 



Fig. 8. An instantaneous snapshot of a front realization at large times. Note the 
existence of the foremost occupied lattice site, on the right of which no lattice site 
has been previously occupied. 

Such a description, put forward in Ref. [86,87], is motivated by the micro- 
scopic picture of the front dynamics at the far end of the front. The idea is 
based on the observation that in every discrete particle and lattice system of 
a fluctuating front, at all times there exists a foremost occupied lattice site 
(f.o.l.s.), which is defined as the one, on the right of which no lattice site has 
ever been occupied before (see Fig. 8). As front propagation in this framework 
is effectively tantamount to forward (rightward) movement of the f.o.l.s., the 
crucial role played by the discreteness effects of the particles and the lattice is 
reflected in the mechanism of the front propagation at the f.o.l.s., where the 
mechanism is not that of uniform translation, but instead, is of "halt- and- g o" . 
A lattice site, which has never been occupied before, attains the status of the 
f.o.l.s. as soon as one particle hops into it from behind (left). In reference to 
the lattice, the position of the f.o.l.s. remains fixed at this site for some time, 
i.e., after its creation, an f.o.l.s. remains the f.o.l.s. for some time. During this 
time, however, the number of particles on and behind the f.o.l.s. continues to 
grow, and as a consequence, the chance of one of them making a hop on to 

24 In this context, we draw the reader's attention to the quantity a in the derivation 
of the large N asymptotics (7) in Ref. [86,87]. The quantities a and c are in fact 
equivalent, both are introduced by hand, and the values of none of them can be 
obtained from a first principle derivation. 
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the right of the present f.o.l.s. also increases. At some instant, a particle from 
behind hops over to the right of the present f.o.l.s.: as a result of this hop, the 
position of the f.o.l.s. advances, or, viewed from another angle, a new f.o.l.s. 
is created on the right of the present one. 
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Fig. 9. A snapshot of a realization of a fluctuating front (jaggered curve) at large 
times, and how a theorist might picture such a front (smooth curve). In this picture, 
as the smooth curve shows, a uniformly translating solution travelling with speed 
vn, given by <f>k(t) = 4>^°\k — V]sft) and obeying the deterministic Eq. (14). However, 
the uniformly translating profile is suitable everywhere but few lattice sites at the 
tip of the front leading up to the f.o.l.s. The uniformly translating region is further 
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subdivided into two parts — in the "linear region" , the nonlinear term 

be neglected. In the "nonlinear region" , however, all the terms of Eq. (5) have to 
be taken into account. 

Having taken this into consideration, the following scenario was proposed in 
Ref. [86,87]: the lattice and finite particle effects lead to a halt-and-go dy- 
namics at the f.o.l.s., while reasonably far behind, the average front "crosses 
over" to a uniformly translating (conditionally averaged) solution (see Fig. 
9). In this formulation, the effect of stochasticity on the asymptotic front 
speed is coded in the probability distribution of the times required for the 
advancement of the f.o.l.s. The goal therein has been to develop a separate 
probabilistic theory for the particle hops to create the new f.o.l.s., and then 
to demonstrate by matching the front solution on the f.o.l.s. with the more 
standard one (of growth and uniform translation) reasonably far behind it in a 
mean-field type approximation at the "tip region" of the front, 25 one obtains 
a self-consistent description of the stochastic and discreteness effects on front 
propagation. Although this formalism is not predictive, it allows one to deal 
with much smaller values of N than it is required for the ln~ N asymptotics 
to be applicable. Moreover, it has the right asymptotics — it shows that as 
the value of iV is increased, one does approach to a uniformly translating front 
profile all the way to one lattice site behind the f.o.l.s., justifying the validity 



25 The tip region of the front spans a few lattice sites right behind the f.o.l.s. and 
also includes the f.o.l.s. (see Fig. 9). 
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of the use of the conditionally averaged uniformly translating profile all over 
the front region, used in Refs. [19,86,87, 122, 125, 126] for the large- N asymp- 
totics (7). Below we provide few mathematical details of this procedure and 
present the key results, while the full details can be found in Ref. [86,87]. We 
note here that the formalism therein was developed for the reaction-diffusion 
system 2X, but in principle it can be worked out for any other microscopic 
model of fluctuating "pulled" fronts. 

The difficulty in prediction of Vn in this formalism is the following. Micro- 
scopically, the selection process for the length of the time span between two 
consecutive f.o.l.s. creations is stochastic, and having followed the movements 
of the f.o.l.s. of one single front realization over a long time at large times, 
and thereafter having denoted the j successive values of the duration of halts 
of the f.o.l.s. by At± : At 2 , . . . , Atj (Atj> > OVj'), one defines the front speed 
as the inverse of the long time average of this time span, i.e., 

-l 

(12) 

Simultaneously, the amount of growth of particle numbers on and behind the 
f.o.l.s. itself depends on the time span between two consecutive f.o.l.s. creations 
(the longer the time span, the longer the amount of growth). As a consequence, 
on average, the selection mechanism for the length of the time span between 
two consecutive f.o.l.s. movements, which determines the asymptotic front 
speed, is nonlinear. 

This inherent nonlinearity makes the prediction of the asymptotic front speed 
difficult. One might recall the difficulties associated with the prediction of 
pushed fronts due to nonlinear terms in this context, although the nature of the 
nonlinearities in these two cases is completely different. In the case of pushed 
fronts, the asymptotic front speed is determined by the mean-field dynamics of 
the fronts, and the nonlinearities originate from the nonlinear growth terms in 
the partial differential equations that describe the mean-field dynamics. On the 
other hand, for fronts consisting of discrete particles on a discrete lattice, the 
corresponding mean-field growth terms in the tip region of the front are linear. 
However, since the asymptotic front speed is determined from the probability 
distribution of the time span between two consecutive f.o.l.s. movements, on 
average, it is the relation between this probability distribution and the effect 
of the linear growth terms that the nonlinearities stem from (and makes it 
very difficult to theoretically deal with). 

If we now denote the probability that a f.o.l.s. remains the f.o.l.s. for time At 
by V(At), then the asymptotic front speed of Eq. (12) is given by 26 

26 In this form, one starts the clock at zero as soon as a lattice site attains the status 
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(13) 

As outlined before, the theoretical expression of V(At) is obtained from the 
crossover of the mean-field type front solution at the tip region to the uni- 
formly translating front solution behind. The front solution in the uniformly 
translating region (of Fig. 9) is given by <f>(°\k — v N At), such that 



v N = {Aty 



d(At) At V(At) 



ML 

d(At) 



= D 



6 {0) + 6 {0) 



(0) 



4°> 



(14) 



where D is the rate of diffusion of particles to their nearest neighbour lat- 
tice sites in the reaction-diffusion system 2X. The uniformly translating 
region can be further subdivided into two parts, a "linear region" and a "non- 



linear region" , where the nonlinear 



do) 



term can and cannot be neglected 



respectively. On the other hand, in the mean-field type description of the 
tip region, the front solution is written as (f^°\k — v^At) + 5(j>k(At). The 
quantities 5<pk{At) are expected to be strictly positive due to the halt-and-go 
motion of the f.o.l.s. — during the time the f.o.l.s. has halted at a given lattice 
site, the accumulation of particles flowing in from behind should result in an 
excess of "particle density buildup" over the uniformly translating solution 
(f><®(k-v N At) at the tip region. Finally, in this formalism, the crossover point 
from the tip region to the uniformly translating region is located at the site, 
where 5<J)k(At) / '<p^\k — v^At) remains negligibly small VAi. 
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Fig. 10. The agreement between the self-consistent theoretical curve of P(At) [86,87] 
and computer simulation results. 

It is clear that to obtain the expression of V(At), one needs the expressions of 
4>(°\k — VNAt) and S(pk(At). On the other hand, the expression of vn itself is 
needed to solve for <j)(°\k — VNAt), and vn can be determined only from V(At) 



of the f.o.l.s. Thereafter, as soon as the next new f.o.l.s. is created, the clock reading 
is reset to zero. 
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as Eq. (9) shows. This indicates that the only way to obtain the expression of 
V(At) is to solve a whole system of equations self-consistently [86,87]. We also 
note here that in this self-consistent formalism, there is an effective parameter. 




k,-3 k,-2 k,-\ k m ' k,-3 k,-2 k,-\ k„, 



lattice index, k lattice index, k 

Fig. 11. Comparison between the 5Nk(At = 0) = N5<j>(At = 0) and the 
(At = 0) = iV0 (o) (k) values obtained from the mean- field type theory of [86,87] 
for four foremost occupied lattice sites (the f.o.l.s. is indexed by k m ) and for N = 10 2 , 
10 3 , 10 4 and 10 5 . Note that as N increases, the 5Nk(At = 0) values compared to 
the N^(At = 0) profile become less and less important. 

The corresponding self-consistent set of equations are highly nonlinear and 
complicated, but due to the presence of the effective parameter in this for- 
malism, the procedure to obtain the V(At) in this formalism is not predictive. 
However, the fact that it generates a probability distribution that agrees so 
well with numerical simulations is indicative of the essential correctness of 
such a description of a fluctuating "pulled" front. The self-consistent theoret- 

oo 

ical curves of P(At) = J dt'V(t') for D = 1 and N = 10 4 , 10 2 , 10 3 and 10 5 (in 

At 

that order) are shown in Fig. 10. The corresponding numerical comparison of 
front speeds are shown in Table I. First, we notice that in the graphs of Fig. 
10, the theoretical curves lie below the simulation histograms at At ~ 2/vn 
- this is an artifact of the matching that one has to carry out for the ex- 
pressions of P(At) below and above At ~ 2/vn. This difference occurs due 
to certain fluctuation and correlation effects [86,87]. Secondly, the agreement 
between the simulation data and the self-consistent formalism is not very good 
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for N = 10 5 — at this value of N, the simulation gets very slow and one has 
to continuously remove particles from the stable phase of the front to gain 
program speed, which affects the -P(Ai) histograms significantly for large At 
values. 



N 


v n (simulation) 


v n (theoretical) 


v N [Eq. (7)] 


10 2 


1.778 


1.808 


1.465 


10 3 


1.901 


1.899 


1.803 


10 4 


1.964 


1.988 


1.925 


10 5 


2.001 


2.057 


1.976 



Table II: Comparison of vjy values: simulation, self-consistent formalism [86, 87] 
[denoted by v n (theoretical)], and that of the large N asymptotics (7). 

io- 2 
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10 6 10 8 10" IO 20 10 30 10 60 10 100 10 180 

N 

Fig. 12. Comparison between the actual stochastic model and the corresponding 
"simplified model" relating to directed polymers in random media [21]. For the 
large N asymptotics of vn and the leading order iV-dependent behaviour of the 
subdominant corrections to it, note the agreement between the results of the two 
models. Here, v m i n is v*, and K approx is the coefficient of 1/ln 2 N in Eq. (7). The 
figure has been modified from the original for greater clarity. 

In addition to good agreements between the self-consistent formalism and 
simulations for the -P(At) curves in Fig. 10, a significant observation of the 
formalism is the following: as the value of N is increased in the self-consistent 
formalism, it is found that the 5<pk/ values at the tip region of the front 
gradually reduces [86,87] (see Fig. 11). The stochastic halt-and-go character of 
the movement of the f.o.l.s., which is usually occupied by 0(1) number of par- 
ticles, however, continues to remain valid for any value of N. This implies that 
for very large N, one approaches the picture of a fluctuating "pulled" front, 
where a uniformly translating mean-field description (5) holds all the way up 
to one lattice site behind the f.o.l.s., while only the dynamics of the f.o.l.s. 
remains a stochastic halt-and-go process. Such a scenario was numerically in- 
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vestigated in Ref. [21] by means of constructing a "simplified model" , in which 
stochasticity was introduced only on the f.o.l.s., while the rest of the front was 
described by a deterministic equation. A comparison of the simulation data of 
the actual stochastic model (the "full model" ) and the simplified model from 
Ref. [21] is shown in Fig. 12. The agreement between the two cases — not only 
at the level of large N asymptotics (7), but also the leading order iV-dependent 
behaviour of the subdominant corrections (c.f. the last paragraph of Sec. 2.2) 

- suggests that the "simplified model" is a very good representation of the 
true stochastic model at asymptotically large iV values. 27 

2.4 The weakly pushed nature of fluctuating "pulled" fronts 

2.4-1 (In) dependence of the front shape and speed on the initial configuration 

One common aspect of propagating fronts in the microscopic models of fluctu- 
ating fronts is that at a given value of At, the more particles there are on the 
f.o.l.s., the more likely it is for the f.o.l.s. to move forward. This indicates that 
if there are more particles put on the f.o.l.s. in these models at short values of 
At, the faster the f.o.l.s. moves. 

In reaction-diffusion models, the number of particles on the f.o.l.s. can in- 
crease by means of two processes: flow of particles from the other lattice sites 
in the tip region behind the f.o.l.s., and growth of particles on the f.o.l.s. by 
themselves without having any input of particles from other lattice sites. In 
the so-called clock model [122-126] or the discrete particle and lattice model 
related to directed polymers in random media [19-21,24,27,28], the number 
of particles on the f.o.l.s. can increase only by means of flow of particles from 
the other lattice sites behind. For the two latter models, when a condition- 
ally averaged steady shape of the front emerges at long times, a complicated 
balance develops between (At) = t^ 1 (i.e., on average for how long the f.o.l.s. 
halts at a given lattice site) and the average flow of particles on to the f.o.l.s. 
from behind, and the balance is established at v N < v*. Such a balance also 
develops for the reaction-diffusion models, but it is more complicated in nature 

- the particles that flow on to the f.o.l.s. can start their own growth, thereby 
further contributing to the increasing number of particles on the f.o.l.s. 

The observation above demonstrates the pushed 28 nature of these fronts, as 

27 Such a statement goes beyond the consideration of simply the front speed, as 
the front diffusion coefficients for the actual stochastic model and the "simplified 
model" show no difference for asymptotically large TV values. We will address the 
issue of front diffusion later in Sec. 2.5. 

28 The above argument is not quantitative, but we have seen earlier that the asymp- 
totic front speed vn differs from v* by only a small amount for large values of N. 
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their speed of propagation depends on the front solution far behind the f.o.l.s. 
The front solution far behind the f.o.l.s. contributes to the movements of the 
f.o.l.s. in the following way: a uniformly translating solution exists behind the 
tip region (see Fig. 9), and for the time the f.o.l.s. halts at a given lattice site, 
the uniformly translating region of the front brings in particles into the tip 
region — ultimately that contributes to how fast the f.o.l.s. moves forward. 
The supply of particles from behind is therefore a very crucial ingredient for 
the front (along with the f.o.l.s.) to propagate at a speed vjy, indicating that if 
the supply of particles on to the f.o.l.s. from behind were cut off, the speed of 
the f.o.l.s. would reduce drastically (this is already illustrated by the discussion 
of the first two paragraphs of Sec. 2.2, where the "inertia" of the f.o.l.s. against 
its forward movement due to the discreteness of particles and the lattice effects 
are mimicked by a cutoff). It also indicates that if the number of particles on 
the f.o.l.s. are increased artificially in these models at small values of At, the 
front should propagate faster. Such a model was investigated in Ref. [88] in 
a reaction-diffusion system with enhanced growth rate of particles in the tip 
region, resulting in a front that propagates with a speed significantly higher 
than v*. 29 

On the basis of the above two paragraphs, the independence of the front speed 
and the conditionally averaged front shape on the initial configuration can be 
argued in the following manner, although no formal proof exists in the litera- 
ture. 30 When a discrete particle and lattice system of a fluctuating "pulled" 
front has not reached its conditionally averaged steady shape, the balance 
between the average value of At and the average flow of particles on to the 
f.o.l.s. from behind does not hold. For example, consider an initial microscopic 
configuration, when the front shape at the leading edge resembles that of the 
functional form exp[— Xk], with A = A> > A*. With such a configuration, "lin- 
ear region" initially propagates with a speed v(X > ) > v*, and bring in more 
particles on to the tip region (and eventually on the f.o.l.s.) than what it does 



This is why we call them weakly pushed. We will see later that such a nomenclature 
is also supported by the spectrum of C VN . 

29 As readers can find out, Ref. [88] studies a continuum space Fisher-Kolmogorov 
equation type model, where the cutoff in the growth function at (j) = e is supple- 
mented by an enhancement in growth between e < 4> < e/r for r < 1. Of course, 
as e — > 0, the evolution equation of the front reaches its pulled front limit (Fisher- 
Kolmogorov equation) . The analysis however shows that first taking the e — > limit 
to obtain the Fisher-Kolmogorov equation and then obtaining the front speed (v*) 
is not right — one has to first calculate the front speed as a function of e and r 
and only then take the limit e — ► 0. Interestingly enough, it turns out that when 
r < r c = 0.283833..., the front speed converges to v*, otherwise the front speed 
can be significantly higher than v*. 

30 As we have seen in Sec. 2.3, the balance between how long the f.o.l.s. halts at a 
given lattice site and the average flow of particles on to the f.o.l.s. from behind is 
very complicated. In view of this, it is unlikely that a formal proof can be obtained. 
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at the steady state, speeding up the f.o.l.s. Subsequently, the "inertia" of the 
f.o.l.s. also gets communicated to the uniformly translating front solution be- 
hind, which then starts to slow down. A combination of these two phenomena 
makes the leading edge less and less steep, until the front finally settles to its 
steady shape (an illustration of it can be found in Fig. 1 of Ref. [16]). Con- 
trast this situation with an initial microscopic configuration that resembles the 
functional form exp[— Xk], with A = A< < A*, such that v(X < ) = i>(A>). In this 
case, the "linear region" initially propagates with the same speed as before, 
but since it is less steep than exp[— A>/c], it does not bring in as many particles 
on to the f.o.l.s. as it does for A = A>, and as a result, the f.o.l.s. does not 
speed up as much as it did before. By means of a simple comparison of couple 
of diagrams contrasting these two cases, one can easily convince oneself that 
the front profile at the leading edge in this case continues to becomes steeper 
and steeper. Although there does not exist any specific simulation data show- 
ing that for arbitrary initial configuration, the front ultimately settles down 
to the same conditionally averaged front shape (and speed), the arguments 
above strongly suggest that it is indeed the case. 



2.4-2 The Large N Asymptotic Spectrum of C VN 

To obtain the large N asymptotic spectrum of C VN , we follow the procedure 
developed in Refs. [56,89], which in turn use the standard route of transforming 
the linear eigenvalue problem into a Schrodinger eigenvalue problem [12,22, 
23,38,92]. To represent a discrete particle and lattice system of a fluctuating 
"pulled" front" as iV — > oo, the idea is to consider the cutoff picture (8) in 
continuum space, developed by Brunet and Derrida. 

In the comoving co-ordinate £ = x — vwt, Eq. (8) in the steady state reads 



aie 



(15) 



with /[0 (O) (O1 = {0 (O) (O - [0 (O) (O]™} @[0 (O) (O - l l N \- N o w consider a func- 
tion (j)(x,t) infinitesimally different from = (/>^(x — v^t) in the comov- 
ing frame, i.e., <fi(x, t) = 0^ (x — v^t) +?7(£, t). Upon linearizing the dynamical 
equation in the comoving frame, one finds that the function r)(x,t) = r}(£,t) 
obeys the following equation: 



dr]_ 

dt ~ L 



V N 



d d 2 5f((f)) 



(16) 



In this formalism, the substitution r)(£,t) = e Et e Vn ^ 2 ^Pe(0 converts Eq. 
(16) to the following one-dimensional Schrodinger equation for a particle in a 
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potential (with h 2 /2m = 1) [12,22,23,38,92]: 



N 



de + 4 



8(j) 



MO = EMO- 



(17) 



In Eq. (17), the quantity V(£) 



\ < 5/(0) 




_ 4 5(f) 


= </>(O)_ 



plays the role of the 



potential. Having denoted the the coordinate of the point where ^ {0) — 1/N 
by £ as before, we have 



v(H) 



4 



■1+n 



\(o) 



n—l 



e(eo-0 — 5(£-£o)+^e(£-£o). (i ? 

uat 4 



The form of the potential V(^) is sketched in the left figure of Fig. 13. Notice 
that </>^(£) is a monotonically increasing function from 1/N at £o towards 
the left, asymptotically reaching the value 1 as £ — > — oo. As a result, for 
£ < £o) a l so increases monotonically towards the left, from — 1 + 
n N l ~ n ~ -7r 2 /ln 2 Af at £ = £ _, to (n - 7r 2 /ln 2 A0 « n as £ -> -oo. At 
£o, there is an attractive 5-function potential of strength 1/vn ~ 1/2 and a 
finite step of height 1. The crucial feature for the stability analysis below is 
the fact that V(£) stays remarkably flat at a value — rr 2 / In 2 N over a distance 
(£o — £i) — IniV, and then on the left of £i, it increases to the value ~ n, over 
a distance of order unity. This is a consequence of the nature of the solution 
(10). 




Fig. 13. Left: The potential V(£) in the Schrodinger operator obtained in the stabil- 
ity analysis. Right: The approximate potential Vb(£) that can be used for calculating 
the low-lying modes for large widths of the bottom well, i.e., for large N. 

The question of stability of the front solution can be answered by studying the 
spectrum of the temporal eigenvalues. If there are negative eigenvalues of the 
above Schrodinger equation, then according to Eq. (16) rj(£,t) grows in time 
in the comoving frame, i.e., the front solution 4> is unstable. On the other 
hand, if there are no negative eigenvalues, then the asymptotic front shape is 
stable, and the spectrum of the eigenvalues then determines the nature of the 
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relaxation of <f>(x,t) to the solution The full spectrum in general de- 

pends on the boundary conditions imposed on the eigenfunctions ipE- Here we 
consider only localized perturbations, for which we need to have r?(£, t) — > as 
£ — > ±00. Because of the exponential factor in Eq. (16), any eigenfunction tpE 
which vanishes as £ — > 00 is consistent with vanishing 77 towards the right 31 . 
However, for £ — > —00, the eigenfunctions ^ need to vanish exponentially fast 
with a sufficiently large exponent, so that when it is combined with the expo- 
nentially diverging term e~ VN ^ 2 , they are still consistent with the requirement 
that 7] vanishes for £ — > —00. For the lowest "energy" eigenvalues investigated 
below, we demonstrate that these requirements are obeyed. 

From the form in the potential, it is clear that the lowest "energy" eigenmodes, 
i.e., the slowest relaxation eigenmodes, are the ones that are confined to the 
bottom of the potential. We notice that among these modes, invariably there 
is a zero mode of the stability operator that is associated with the uniformly 
translating front solution of a dynamical equation like Eq. (8): since </>^(£) and 
(/>(°)(£ + a) are solutions of Eq. (14) for any arbitrary a, we find by expanding 

to first order in a that ipo(0 — e VN ^ 2 is a solution of Eq. (17) with 

eigenvalue E = 0. 32 From the result (19) for the asymptotic front solution, 
we then immediately get, to dominant order, 



^0 ~ sin[7rA*(£ - 6)/ hi AT] ~ & < £ < (19) 

Furthermore, since (£) is a monotonically decreasing function of £, the solu- 

tion ipo(£) = e VN ^l 2 is nodeless. Since we know from elementary quantum 

mechanics that the nodeless eigenfunction has the lowest eigenvalue, this im- 
plies that all the other eigenvalues of Eq. (17) are positive, i.e., the solution 
0(°)(O is stable. 

The spectrum of eigenvalues of Eq. (17) for E > determines the decay prop- 
erty of localized perturbations r)(£, t) in time. We notice that for E > v%/4 m 1 
(the value of the potential on the far right), the spectrum of eigenvalues are 
continuous. However, we are particularly interested in the smallest eigenvalues 
E m > for small m, since these are the eigenmodes that decay the slowest in 
time. These are the eigenvalues associated with bound states in the potential 
well. As iV — > 00, the bottom well of the potential becomes very wide: its 



31 The fact that eigenfunctions ipE which diverge as £ — > 00 are allowed, means 
that there are admissible eigenfunctions which are not in the Hilbert space of the 
Schrodinger operator. See e.g., Ref. [38] for further discussion of this point. 

32 We will later see in Sec. 3.1.2 that — - — is also known as the so-called Goldstone 
mode [94,109,110]. 
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width diverges as In N. As we know from elementary quantum mechanics, the 
lowest "energy" eigenfunctions then become essentially sine or cosine waves 
in the potential well with small wave numbers k and correspondingly small 
"energy" eigenvalues. Based on the fact that the potential V(£) on the left 
rises over length scales of order unity, we now make an approximation. In the 
limit that the bottom well is very wide and the k values of the bound state 
eigenmodes very small, it becomes an increasingly good and an asymptotically 
correct approximation to view the left wall of the well as a steep step, sketched 
in the right figure of Fig. 13 — one can therefore approximate the potential 
by 



V (0 = n [1 - 0(0] - ^©(0 [1 " ©(£ " &)] 

~*(e-eo)+e(e-&). (20) 

On the right hand side, there is an attractive delta-function potential at the 
point where the potential shows a step to a value close to 1. It is easy to check 
that the prefactor of the delta-function of 1/2 is not strong enough to give 
rise to bound states with E < 0, and as a result, for very large values of N, 
the low-lying eigenmodes approach sine waves with nodes at the position of 
the walls of the potential, 33 ip m ~ sin [q m (£ — £i)]. The condition that these 
solutions have nodes at the right edge of the well then yields 



(m+l)7r ^ (m + l)n 



implying that the corresponding eigenvalues are given by 
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E m ^ [{m+1) l ~ 1]?r2 m = 0,1,2,... (22) 
ln 2 iV ' ' ' V ; 

Here, the first term between square brackets comes from the "kinetic energy" 
term q 2 , while the second term originates from the value of the potential at 



33 More explicitly, if we write the solution within the well as Asm[q(£ — £i) + B], and 
for £ > £o as A + e~^~^ (which is correct to lowest order in e), then we get from 
the boundary conditions at £o : 2g , cotg[g(£o — £i) + B] = —1; likewise, at the left 
boundary, of the well, we get qcotgB = const., where the constant is determined 
by the size of potential step. For ^ - 6 » 1, the small-g solutions are those with 
B -> 0. 

34 Notice that for the eigenvalues in Eq. (22), the corresponding eigenmodes iPe{Q 
decay as exp[— y/n |£|] for £ — ► — oo and as exp[— vn£/2] for £ — > oo, which make 
e VN ^/ 2 ^E(Q to go to zero for £ — > ±oo, satisfying the boundary conditions discussed 
previously. 
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the bottom. Note that for m — 0, the eigenmode sin go with eigenvalue Eq 
is indeed the same as the zero eigenmode of Eq. (19) with q = n\*/lnN, 
which we calculated from the shape of the front solution in the leading 
edge. Besides verifying the consistency of our approach, this also confirms 
that there are no corrections to Eq. (22) for m = 0: for m = 0, it will yield an 
eigenvalue zero to all orders in N. Therefore, the smallest nonzero eigenvalue, 
which governs the relaxation of the front velocity and profile to the asymptotic 
ones is E 1 with relaxation time T\ given by 



3tt 2 
In 2 iV 



(23) 



The time convergence of the asymptotic front speed for the reaction-diffusion 
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Fig. 14. The time convergence of the asymptotic front speeds: for Eq. (8) [short 
dashed line], and for the stochastic reaction-diffusion system X+Y^ 2X (solid 
line) [56]. They follow corresponding curve for the deterministic equation (1) with 
n = 2 (long dashed line) for a while, and then they suddenly break off to settle down 
to their asymptotic values. As the long dashed line represents the power law conver- 
gence t~ l for pulled fronts, such a behaviour for the solid and the short dashed curve 
is expected for a very slow exponential decay of Eq. (23). Note that the notations 
have been changed from the original to suit the present article. Here 5v = v* — v. 

system X+Y^ 2X is shown in the simulation data from Ref. [56] (Fig. 14), in 
support of Eq. (23). As can be seen therein, as a function of time, the front 
speed for Eq. (8) [short dashed line] and for the stochastic reaction-diffusion 
system X+Y^ 2X (solid line) follow corresponding curve for the deterministic 
equation (1) with n = 2 (long dashed line) for a while, and then they suddenly 
break off to settle down to their asymptotic values. As the long dashed line 
represents the power law convergence t" 1 for pulled fronts, such a behaviour 
for the solid and the short dashed curve is expected for a very slow exponential 
decay of Eq. (23). 

We have noted earlier in the introduction that for a pushed front, the lowest 
lying eigenspectrum is gapped. For discrete particle and lattice systems of a 
fluctuating "pulled" front with A" — > oo, according to the above analysis, the 
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spectrum is also gapped, but the gap is very small, as it is oc ln~ 2 N. This is 
another confirmation of these fronts being weakly pushed. Moreover, Eq. (22) 
also confirms that as N — > oo, the gap between the spectral lines decreases as 
In -2 N, which is consistent with the fact that for a pulled front N = oo and 
the spectrum becomes gapless (see footnote 12 in this context). 

In connection to the gapless spectrum in the limit N — > oo, a curious obser- 
vation was made in Ref. [56]. Put simply, the observation is the following: in 
the time convergence of the front speed to its asymptotic value v*, the leading 
edge increases in length ~ y/t leading to the uniform t" 1 convergence (4). 
For a discrete particle and lattice system of fluctuating "pulled" fronts with 
N — > oo, the length of the leading edge in Eq. (10) increases as In N, which ul- 
timately leads to the 1/ In 2 N deviation of the asymptotic front speed vn from 
v* (although in this similarity the prefactors do not match correctly). Based on 
these observation, Kessler and co-workers suggested [56] that the initial condi- 
tions, which are responsible for the front speed selection of a pulled front, act 
as an effective time- dependent cutoff' beyond which the steady state equation 
no longer holds, and the time-dependent cutoff goes to zero with increasing 
time 35 . To what extent there is a deep physics underlying this observation, 
however, is not clear. 



2.5 The Large N Asymptotic Scaling of Front Diffusion Coefficient Df for 
Fluctuating "Pulled" Fronts 

So far, we have only dealt with (and observed) how the discreteness of particles 
and the lattice sites affect the asymptotic front speed of fluctuating "pulled" 
fronts. With the idea behind random wandering of fronts has been introduced 
already in Sec. 1.2, we now address the issue of their diffusive wandering. 

It turns out that the theoretical concept of random wandering of fronts is 
actually an extremely tricky issue, the full flavour of which will only be clear 
in Sees. 3 and 4.2. However, operationally, i.e., in a computer simulation, the 
front diffusion for a front consisting of discrete particles on a lattice is rather 
straightforward to measure. Consider for example the fluctuating "pulled" 
front in the so-called clock model [122, 124-126]. This model, which arose in 
the calculation of the largest Lyapunov exponent for a gas of hard spheres 
[122, 124-126], one considers a set of iV clocks with integral clock readings 

k = 0,1,2, The number of clocks with a certain reading k at time t is 

rik(t). Any two clocks can "collide" continuously in time; the post-collisional 
readings of a pair of mutually "colliding" clocks with pre-collisional readings kj 
and kj are both updated to max(/cj, kf) + 1. In terms of <f>k(t) = Hk'=k n k l if)/N 



A similar observation was also made therein about pushed fronts with cutoffs. 
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(with time t redefined to have a unit mean single clock collision frequency), 
the relevant asymptotic solution corresponding to an initial state in which all 
clock readings are finite is a front propagating into the state (fik = for large 
k. In the mean-field limit, the dynamics is that of a pulled front propagating 
with speed v* = 4.31107. . . [122,124-126]. 
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Fig. 15. Simulation data from Ref. [21] for the front diffusion coefficient, for both 
the microscopic stochastic model relating to directed polymers in random media 
("full" model) and the "simplified" model. In the above figure, Dn should be read 
as Df. The figure has been modified from the original for greater clarity. 

A straightforward way to measure the front speed and the front diffusion 
coefficient in a computer simulations of the clock model is by tracking the 
position of the centre of mass of the clock distribution in individual realizations 
as a function of time. With the centre of mass location for the clock distribution 
in an actual realization r at time t given by 

S T (t) = N- 1 £ kn k (t) = £ Mt) , (24) 

k k 

one has v N = (S T (t)), and D f = lim^^ d([J t t+T dtS r (t) - v N T] 2 )/dT [the 
angular brackets denote first an average over all possible evolution sequence 
for a given initial realization (at time t), and then a further averaging over 
the ensemble of initial realizations]. This was indeed the same method used 
by Brunet and Derrida to obtain the front speed and front diffusion coeffi- 
cient [21] for a model which is slightly different from the clock model. The 
superb cleverness in their measurement of the front diffusion lies in the simu- 
lation algorithm they used — this algorithm allowed them to take iV as high as 
10 150 , and they demonstrated that the iV — > oo asymptotic scaling of front dif- 
fusion coefficient Df behaves as 1/ln 3 N [21]. Moreover, a "simplified model", 
where the fluctuations are randomly generated only on the instantaneous fore- 
most occupied lattice site (i.f.o.l.s.), 36 was found to exhibit the same 1/ In 3 TV 

36 The i.f.o.l.s. is simply the location of the lattice site on the right of which all 
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asymptotic scaling of Df as the full stochastic model [21]. However, we note 
that in the absence of any studies on the diffusion coefficient of fluctuating 
"pulled" fronts on other microscopic models, to what extent this scaling is 
independent of the model is not known. 

For the clock model, we will provide a argument of the l/ln 3 iV asymptotic 
scaling of Df in Sec. 4.2, while presently we simply obtain the formal Green- 
Kubo expression for Df. Needless to say, dividing the derivation in this manner 
is complicated for the readers. Nevertheless, the reason behind postponing 
the derivation of the l/ln 3 iV asymptotic scaling of Df till Sec. 4.2.2 is the 
following: it will require certain inputs and comparisons with the stochastic 
differential equation used in the Langevin-type field-theoretical approaches of 
fluctuating fronts, which we will discuss in Sees. 3 and 4. 

The central theme to obtain Df is that there is stochasticity in the front 
evolution at two levels: first, in any snapshot of the clock model realization 
r at time t, its shape {nk(t)} r fluctuates around {n^\t)} = {n^(^)}, where 
n ^°H0 = 0^(0 — — 1). In the subsequent front evolution, in an average 
or mean-fieldish sense, the clocks are chosen from a "distorted distribution" 
{nk(t)} T rather than from {nf\t)}, and this introduces fluctuations in the 
front speed Sv r (t) [measured by tracking the centre of mass of the clocks]. Let 
us denote, by 5v Tjia {(t), this average fluctuations in the instantaneous front 
speed around v N at time t due to the shape fluctuation of the front at time t. 
From the actual clocks that constitute the front, at time t, we have 

v N + Sv TM (t) = E t 2 " <WJ - k i\ + 2 )KK- - <WJ)/2AT 2 • (25) 

Equation (25) is obtained by applying a probabilistic argument in the mi- 
croscopic dynamics of the clock model — in a collision between two clocks 
with readings ki and kj, the mean clock reading increases by an amount 
(\ki — kj\+2)/N, the probability of having such a collision is [2 — 5^^] [ n kX n kj ~ 
Ski,kj)]/W(N — 1)], and finally in a unit time, on average, (iV — l)/2 collisions 
occur. The S^k/s arise due to the indistinguishability of the clocks with the 

same readings. When we write {nk(t)} v = nf\t) + {5nk(t)} T , the expression 
of Vj\f is obtained from Eq. (25) by replacing {nk{t)} T by {n^\t)} in the first 
step of Eq. (25), implying that at the leading order 

, r 1 \h- kj\ + 2 4? [4? -Kk 3 ] 

Sv TM (t) = 2J2- h^} ^ ^ x 

i,3 



lattice sites are instantaneously empty. It is different from the f.o.l.s. of Sec. 2.3. 
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The idea behind writing 5i> r ,mf(0 hi this manner is to emphasize the fact that 
despite several factors of N involved in the expression of vn, its magnitude 
is of 0(1); and therefore, the order of magnitude estimate of 5v T>m f(t) must 
behave as the ratio of 5rik(t) and n^\t) in realization r. 

Secondly, beyond (the mean-fieldish) 5v rjm f(t), there is another source of the 
instantaneous front speed fluctuation that depends on what is the precise 
sequence of random numbers that updates the clock distribution in an in- 
finitesimally short time window after t, i.e., between times t and t + e for 
an infmitesimally small e. Let us denote the front speed fluctuation due to 
the stochasticity in the front evolution at time t by 5v T)S (t + ) [this depends on 
the updating sequence of random numbers s immediately after t}. The total 
instantaneous front speed fluctuation Sv T (t) around vn is therefore given by 

Sv T (t) = S T (t) - v N = 5v rM (t) + 5v TiS (t+) . (27) 

Notice that the definition of 5v TjS (t + ) in this manner automatically warrants 
the condition that the average of 5v TjS (t + ) [but not Sv T (t)] over all possible 
collision sequences after t is identically zero. A further average of Sv T (t) over 
an ensemble of realizations at time t, however, must yield zero. 

The formal expression of Df is then finally obtained using the Green-Kubo 
formula 

T 

D f = \ lim / dt' ((Sv.it) 5v v {t + t')) s ) r . (28) 

2 T—>co J 


In Eq. (28), one first needs to average over the random number sequences 
chosen to update the front after time t (subscript s), and then average over 
the initial realizations at time t (subscript r). 

2.6 The Case of Small Values of N 

By now we have seen that the extent of dependence of the front speed on the 
details of the models increases with decreasing N. This in itself is not surpris- 
ing when one realizes that the relative strength in fluctuations in the number 
of particles increases with deceasing iV and that the fluctuation character- 
istics are indeed controlled by the specific details of the model. In Sec. 2.6, 
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we discuss propagating fronts in discrete particle and lattice systems where 
the average number of particles N per lattice site in the stable phase of the 
front is "small" . Studies on front propagation at small values of N, however, 
have been limited to only the X±=f 2X and the X+Y-^2X reaction-diffusion 
systems [5, 14, 34, 54, 55, 71-74, 80, 90, 96, 116, 118, 119]. 37 The mean-field de- 
scription of all the one-dimensional models studied in these references yield 
the Fisher-Kolmogorov equation on a lattice, given by (with a redefined time 
scale and in full generality of parameters) 



= D [<t> k+l + fc _x - 20 fc ] + r(Y0 fc - 4>\) , (29) 

where r is the rate of reaction, D is the rate of diffusion of the particles to their 
nearest neighbours, and <pk = {{Nk))/N. In the stable phase of these fronts 
there are N particles per lattice site on average, and in the mean-field limit, 
the front becomes a pulled front with speed v* (when the space continuum 
limit is taken, then v* = 2\ / DrN). 



2. 6. 1 Summary of the Models Studied and Known Results 

In the various models studied, the value of iV varies between N <C 1 to iV ~ 20. 
Section 2.6 is divided into several subsections with progressively reduced value 
of N. One common trend in all these (one-dimensional) lattice models is to 
identify the front position with the instantaneous position of the foremost 
particle — do notice that this is identical to the i.f.o.l.s. defined in Sec. 2.5, 
and therefore, it is not the same as the f.o.l.s. defined in Sec. 2.3. 

(i) Reaction-diffusion model X+Y^ 2X on a lattice, where both X and Y 
particles can diffuse to their next nearest neighbour lattice sites with the 
same diffusion rate D; any Y particle coming into contact with an X particle 
on the same lattice site is instantaneously converted to an X particle [71]: In 
this model, initially, the entire lattice contains only Y particles, and one X 
particle is introduced at the left end of the lattice at t — 0. Soon afterwards, 
a front consisting of X particles propagates into a region full of Y particles. 
The mean- field description of this model yields a reaction rate r — > oo, 

37 In general, any model of discrete particle on a lattice with a small number of 
particles per lattice site, relatively strong fluctuations render mean-field descrip- 
tions inapplicable, and there exist numerous accounts of it in the literature. In the 
context of reaction-diffusion systems, on general issues of suitability of mean- field 
descriptions of discrete particle and lattice systems at small number of particles per 
lattice site, readers may be interested to consult Refs. [6-9,31-33] and the references 
cited therein, and also references cited in Ref. [71]. The reaction-diffusion systems in 
these references are of 2X, X+Y^2X, X+X^X, X+X^inert and X+Y^inert 
type. 
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resulting in an infinite mean-field front speed. However, in the actual lattice 
model for iV — > 0, this is not the case: dimensional arguments, confirmed 
by simulation shows that oc DN, while the front diffusion coefficient 
Df is oc D, with proportionality constants of 0(1) [71]. In a discrete-time 
equivalent of this model, where time is measured in discrete units of D 
(thereby rendering D — 1) and the system is parallelly updated, the front 
speed has been shown to behave as [118, 119] 

v N = 1- e- N ' 2 , (30) 
with vn saturating at 1 for large N. 

(ii) Reaction-diffusion model X+Y^ 2X, where on a lattice both X and 
Y particles can diffuse to their next nearest neighbour lattice sites with 
the same diffusion rate D, and a Y particle coming into contact with an X 
particle on the same lattice site is converted to an X particle with a finite 
rate r: although one might expect the front speed at small N to have no 
resemblance with the mean-field expression v*, it turns out that for some 
small value of N and r C 1, surprisingly the mean- field front speed v* is 
recovered again [73]. 

(iii) A reaction-diffusion model similar to the one in (i) above was considered 
in Ref. [74], with the X and Y particles have different diffusion rates, say 
Dx and Dy, then for Dx ^ 0, the front speed depends only on Dx, and 
not on Dy- At Dx = however, there is no uniformly translating front 
solution [74]. 

(iv) With only one particle allowed per lattice site, one of the major variants 
of the model that has been considered in the literature is the reaction- 
diffusion system X^ 2X. In this system, only three basic moves are allowed 
(see Fig. 16): (a) an X particle can diffuse to any one of its neighbour lattice 
sites with a diffusion rate D, provided this neighbouring site is empty; (b) an 
X particle can give birth to another one on any one of its empty neighbour 
lattice site with a birth rate e; (c) any one of two X particles belonging to 
two neighbouring filled lattice sites can get annihilated with a annihilation 
rate W. 

For all values of the parameters D, e and W, there exists an exact relation 
between the front speed vn and Df, namely [90] 

V f + D f = D + e. (31) 

For W = and D/e -> oo [14], and for W ^ and y/De/(e + W) -> oo [81], 
the front speed is given by its mean-field limit, v* = 2\j2Ds, but when D and 
e are of the same order of magnitude, then there are tremendous deviations 
from this mean- field theory results [14,54,55]. In Ref. [55], these differences 
were studied numerically, and it was shown in [54] that the differences from 
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the mean-field theory results can be explained by means of a "self-consistent 
two-particle approach". For W ^ 0, the W = D case was considered in 
Refs. [5,34] in continuum space, where the front speed is obtained as v = e 
and front diffusion coefficient Df — D. 

(a) | y\ J I (b) I JJ»J I I (c) 1 • I • I 



• I • I L 





Fig. 16. The microscopic processes that take place inside the system: (a) a diffusive 
hop with rate D to a neighbouring empty site; (b) creation of a new particle on a 
site neighbouring an occupied site with rate e; (c) annihilation of a particle on a 
site adjacent to an occupied site at a rate W. 

The full lattice model where D, W and e are all different, and D/W and 
D/e are of the same order of magnitude was considered in Ref. [90]; in this 
case, one has a good approximate solution [90] 

e e(e + D) , „ (e + 2W)(e + D) , , 
N 777 v n ^ 777 1 and D f ~ ^ — - — J -\— '- . 32 

e + W e + W s 2(e + W) y ' 

(v) Fluctuating fronts in higher than one spatial dimensions: all the known 
studies on propagating fronts in discrete particle and lattice realizations 
in higher than one spatial dimensions have been on X+X^ 2X reaction- 
diffusion system [80,96,116], or on X+Y— >• 2X. We have seen earlier that for 
fluctuating fronts in one spatial dimension, the front position suffers diffu- 
sive wandering around its mean. In higher than one spatial dimensions, the 
front wandering properties change: it is then possible to formulate the prob- 
lem along the lines of interfacial growth phenomena [10,47,53,60,61] and 
study the associated scaling properties. However, there exists tremendous 
ambiguity regarding how to properly characterize the wandering behaviour 
of fronts in higher than one dimensions [80,116], although the results indi- 
cate that any generic upper critical dimension must be higher than 4 [80]. 

In the following sections 2.6.2 through 2.6.6, we provide detailed descriptions 
of points (i) through (v) above, one point in each section. 



2. 6.2 Front Propagation in the Diffusion- Limited Irreversible Process X+ Y— > 
2X 

In Sec. 2.6.2, we elaborate the model introduced in point (i) of Sec. 2.6.1. It was 
first introduced and studied in the N — > limit in Ref. [71]. In this model, 
the front position is (naturally) identified by the location of the foremost 
(rightmost) X particle, and the only way the front (consisting of X particles) 
can propagate into the region occupied by the Y particles is by diffusion. Front 
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propagation in this lattice model is therefore diffusion-limited, which results 
in a front speed oc N, as opposed to an infinite front speed as the mean-field 
theory predicts. 

To solve for front speed and front diffusion, the following dimensional argu- 
ment was proposed in Ref. [71]: there are only two dimensionful parameters in 
this problem: N, which scales as inverse length, and D. Therefore, a = v^/DN 
is a dimensionless quantity. On the other hand, wandering of the front region 
is diffusive, and if one denotes the front width due to the diffusive broadening 
by w(t), then w 2 (t) ~ Dft, and from the definition of w(t) and using the 
parameters of the problem N and D, one can define another dimensionless 
quantity (3 = w 2 (t)/Dt = Df/D. These arguments suggest that oc DN 
and Dj (x D. These conjectures were confirmed in simulation, and the values 
of a and j3 were numerically obtained to be equal to 1.05 and 3.57 respectively 
(with uncertainties in the last digits), using N values 0.025, 0.05 and 0.1. 38 

In a continuum space formulation, in the comoving frame of the foremost X 
particle (£ = Nx), the probability density Pxy(0 that the first Y particle is at 
a distance £ on the right of the foremost X particle, and the probability density 
PxxH) that the first X particles on the left of the foremost X particle were 
also calculated in Ref. [71]. In this frame, the foremost X particle is stationary 
(say at £ = 0), and the Y particles in the region £ > and the X particles 
in the region £ < undergo uncorrelated diffusive movements towards £ = 0. 
The density of particles n(£) (Y particles for £ > and X particles for £ < 0) 
are then obtained from 39 



dn , T dn „ „ t9 d 2 n , , 

-- VnN - = 2DN* w , (33) 

The stationary state solution of this equation with the appropriate boundary 
conditions then yields 



Notice that the front really propagates due to the diffusion of X particles, in the 
model of Ref. [71], if the two species X and Y have different diffusion rates, then 
one would simply expect vn = aDxN with the same value of a as deduced above. 
That this idea is actually correct was verified by simulations [74] . We will return to 
this point in Sec. 2.6.4. 

39 In Eq. (33) we need to use the diffusion rate of X particles for £ < (or the 
diffusion rate of Y particles for £ > 0) in the comoving frame of the foremost X 
particle. Therefore, these diffusion rates are actually relative diffusion rates of a pair 
of already diffusive particles, and the value for both of them is 2D. Notice however 
that in Ref. [71], D was used as the relative diffusion rate as opposed to 2D. As a 
result, the solutions of Eq. (33) that we obtain here, which are the correct ones, are 
not the same as the ones obtained in Ref. [71]. 
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ny (£ > 0) = N[l - exp(-v N £/[2DN])] — N [1 — exp(-af/2)] and 
n x (£<0) = N. (34) 



From Eq. (34), the probability p(£) of finding a Y particle at a distance £ 
on the right of the foremost X particle, is obtained from the integral equa- 

tion n Y (i)/N = p(g) + J pffldg, leading to the solution p(f) = a(e~ a ^ 2 - 

o 

e~ ? )/(2 — a). The probability 0(f) of finding a X particle at a distance f on the 

left of the foremost X particle, is similarly obtained from the integral equation 
o 

n x {()/N = 0(f) + J 0(g)dg, yielding 0(f) = ae~ a ^/ 2 /2. As already men- 

tioned, these solutions are different from the ones obtained in Ref. [71] (see 
footnote 39). 

Notice that Eq. (33) does not determine the front speed in any way, instead 
Vn enters as an input for the solution of p(f) and 0(f). At best, vn has to 
be determined from somewhere else. In a model closely related with the the 
one analyzed in Ref. [71], where time was discretized in units of -D -1 , this 
incompleteness was removed: in the iV — > limit, was deduced from first 
principles to be equal to N/2 at the leading order of N [118,119]. In this model, 
in a unit time, each particle moves one lattice distance randomly towards its 
right or left, and the method of "parallel updating" 40 was employed in the 
simulation algorithm. The argument to derive the expression of front speed in 
this model in the N — > limit is the following: the front speed is essentially 
determined by considering the two foremost X particles (which are separated 
by an average of 1/N lattice sites). 41 In fact, the only non-zero contribution 
to the front speed occurs when the second X particle is on the lattice site just 
behind the foremost one, or on the same lattice site as the foremost one (see 
footnote 42 for an explanation). In the latter case, one of the X particles at 
the location of the front is tagged to be the foremost X particle and the other 
is then tagged as the second X particle. In the comoving frame of the front, 
let us tag the location of the foremost X particle by k = 0, and denote the 
probability that the second X particle is k lattice sites behind the foremost X 
particle by It is then easily seen that 42 



40 The idea behind parallel updating in this model is that in a unit time, all particles 
move diffusively one lattice distance randomly either towards their left or the right, 
and when this update brings an X particle in contact with any number of Y particles 
on a lattice site, all the Y particles on that lattice site are instantaneously converted 
to X particles. 

41 In the limit N — ► 0, the role of a third (or more) particle near the foremost X 
particle can be neglected; they only provide higher order corrections. 

42 Equation (35) is obtained by noticing that in a unit time, the two foremost X 
particles can execute a total of four moves, each with probability 1/4. When the two 
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The steady state solution of p_i and po are obtained from the master equation 
for p k 's around foremost X particle (i.e., k = —2, —1, . . . , 2 etc.), yielding 43 



p_i = N, po = N/2 and ^ = N/2. 



(36) 



The above derivation clearly does not hold for higher values of N, as one 
needs to consider the case where a lattice site is occupied by multiple number 
(many more than two) of X particles. Such an analysis was also presented in 
Ref. [118, 119], in which, having denoted the probability of having the k-th 
lattice site occupied by n k number of particles by p(n k ), one has [Eq. (37) is 
obtained by following the logic of footnote 42] 
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respectively for the probability of the front to step backward one lattice site, 
to remain where it presently is, and to step forward one lattice site in one 
unit time. To obtain the front speed from Eq. (37) however, one needs the 
probability distributions p(n ) and p(n_i) [118, 119]. In Ref. [118, 119], while 
p{n_i) was obtained from the simulations as a Poissonian with average N, a 
truncated Poissonian 44 was used for p(no), yielding the result [notice that one 
recovers Eq. (36) for from Eq. (38) in the limit N — > 0] 



foremost X particles are separated by one lattice site, then the front moves forward 
one lattice distance with probability 1/2, moves backward one lattice distance with 
probability 1/4 and remains stationary with probability 1/4. This explains the origin 
of the p-i/4 term. On the other hand, when the two foremost X particles on the 
same site, then the front moves forward one lattice distance with a probability 3/4, 
and moves backward one lattice distance with a probability 1/4 — giving rise to 
the po/2 term. When the second X particle is more than one lattice site behind 
the foremost X particle, then the front moves forward or backward by one lattice 
distance with equal probabilities 1/2, which is why p_ k for k > 1 does not contribute 
in Eq. (35). 

43 For the master equation and its solution, we urge the readers to consult Ref. 
[118,119]; although here we point out the similarity between Eqs. (36) and (34) - 
in both of these (very similar) models, the density of X particles behind the foremost 
one is N. 

44 The truncated Poissonian distribution for the probability of uq number of particles 
at k = is given by p T {n ) = (1 - 5„ 0i0 )(l ~ e~ N )~ 1 n e~ n ° . 
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Fig. 17. Equation (38) vs. simulation data for front speed, adopted from Ref. [118, 
119]. Here v should be read as vn- The dotted curve shows the agreement between 
theory and simulation if one uses a Poissonian distribution also for p(no). Notations 
have been modified from the original to maintain consistency with the text. 



2.6.3 Reaction- diffusion Model X+Y—+2X at Small N with Reaction Rate 
r -> 

In Ref. [73], the reaction-diffusion model X+Y— >2X was numerically studied 
at small values of N, (N = 0.1,0.3,1,5 and 10) and reaction rates r < 1 
(r < 0.02). The details of the microscopic model are already described in 
point (ii) of Sec. 2.6.1. 

With the front position identified by the position of the foremost X particle, 
the front speed at different values of iV and as a function of reaction rate r 
is shown in Fig. 18 (p and c should be read as r and N respectively), where 
v(r) = v N (r)/2\ / DN. If the mean-field theory were applicable in this model, 
then v{r) should scale as y/r (dashed curve of Fig. 18). Figure 18 reveals that 
typically most values of N yield strong deviation from mean-field result for 
v(r), except for N = 10 for all values of r. Reference [73] also claimed that 
the trend below r ~ 0.001 indicates that the mean-field result front speed for 
r < 0.001, however there are too few data points below r ~ 0.001 in Fig. 18 
to arrive at a definitive conclusion. 

No theoretical argument for why the N — 10 data for the front speed follows 
the mean-field theory results so well was provided in Ref. [73] (neither is it 
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Fig. 18. Front speed at different values of N and as a function of reaction rate r (p 
and c should be read as r and N respectively), where v{r) = vn(v)/2\/ DN. Dashed 
curve: mean-field result v(r) = \Jr. 

clear from a simple argument). Moreover, since Fig. 18 indicates that the 
mean-field results are closely followed for iV = 10, it is a natural extension to 
investigate, by means of simulations, if the front profile also follows the mean- 
field result. However, this point was not investigated in Ref. [73]. Instead, 
the front profile was investigated for N = 0.1, and from the plots that were 
provided in Ref. [73], the front profile seems to behave as 



(39) 

The plots that were supplied in Ref. [73] in support of Eq. (39) are shown 
in Fig. 19. Figure 19(a) shows the front the front profile (here A should be 
read as </>(£) and x should be read as £) for three values of the reaction rates 
r = 0.0025,0.005 and 0.01 (left above from left to right). The dashed lines 
in Fig. 19(a) correspond to Eq. (39) with numerically determined a, b and 
a. The plots of f3(£) = yTog 10 | log lo [0(£)]| are shown in Fig. 19(b): the right 
above left to right plots supposedly correspond to r values 0.0025, 0.005 and 
0.01 respectively. 
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2.6.4 The Case of Different Diffusion Rates for X and Y Particles in the 
Model of Point (i) 

In Ref. [74], a model similar to the one proposed in Ref. [71] was studied. In 
this model, the two species X and Y particles have two different diffusion rates, 
namely D x and D Y . Both the cases D x 7^ and D x = were considered 
in the ensuing analysis. The lattice is initially filled with Y particles with an 
average number of particles N per lattice site and at t = 0, an X particle is 
introduced at the leftmost lattice site, just like the model discussed in point 
(i) of Sec. 2.6.1. However, only one particle is allowed per lattice site in this 
model, and when two neighbouring lattice sites are occupied by one X and one 
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Fig. 19. Plots in support of Eq. (39): (a) front the front profile (here A should 
be read as </>(£) and x should be read as £) for three values of the reaction 
rates r = 0.0025,0.005 and 0.01 (left above from left to right). The dashed lines 
correspond to Eq. (39) with numerically determined a, b and a. (b) Plots of 
/?(£) = A/logxo | logio[<A(£)]l are shown in Fig. 19(b): the right above left to right 
plots supposedly correspond to r values 0.0025, 0.005 and 0.01 respectively. 

Y particle, the Y particle is instantaneously converted to an X particle. As 
pointed out in Ref. [74], the restriction that only one particle is allowed per 
lattice site is not necessary, other than when D x = 0. 

In this model, the D x ^ does not reveal any surprises (see footnote 38 [74]). 
Moreover, the density profile of Y particles on the right of the foremost X 
particle (located at £ = 0) in a continuum-space formulation (£ = Nx) turns 
out to be n y (0 = N[l - exp{-v N £/[(D x + D Y )N]}] [74] (see Fig. 20) — this 
is also consistent with Eq. (34). 
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Fig. 20. Comparison of the density profile of Y particles on the right of the 
foremost X particle in a continuum-space type theory g(£) = ny(£)/Y = 
1 — exp{— v x ^/[(Dx + Dy)N]} and simulation data for two different N values. 
Triangles: D x = 0.2, D Y = 0.8, circles: D x = 0.8, D Y = 0.2. 



The case of D x = 0, however, is more complicated as we explain below. In 
this case, the front propagates only due to the diffusive flux of Y particles 
towards the foremost X particle, and at any time t, with the front position 
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defined by the location Xf(t) of the foremost X particle in a continuum-space 
formulation, one has [74] 



v N (t) 



D 



Y 



dny(x, t) 
dx 



Xf(t) 



(40) 



while the density profile of Y particles on the right of the foremost X particle 
is obtained from the analogous form of Eq. (33) as (notice that the rate of 
diffusion of the Y particles towards the foremost X particle is Dy) 45 



dn dn d 2 n 

dt Vn dx Y dx 2 



(41) 



Equations (40) and (41) together are then solved with the ansatz ny(x, t) = 
Ng{[x—Xf(t)]/ y/Dyt} and a yet unknown function g [74]. This ansatz actually 
suggests that Xf(t) ~ y/t, and that the front speed i>jv(£) ~ t^ 1 ^ 2 is actually 
time-dependent. At low density of particles per lattice site N, the solution of 
ny(x,t) turns out to be the error function g{[x — Xf{t)]/\/Dyi} = erfjjrr — 
Xf(t)}/[2y/D Y t}}. These aspects were verified in Ref. [74] (see Fig. 18). 




0.15 



0.10 



X 

o 



0.05 



0.00 




2000 



Fig. 21. The plots of Xf(t) vs. t in log-log scale, and that of ny(x,t) at t = 10000 
units. Here c(x) should be read as ny(x). Dashed lines correspond to Xf(t) ~ s/i 
and Eq. (41) respectively. 



5 One must be aware of the fact that the usage of Dy as the rate of diffusion 
of the Y particles towards the foremost X particle is valid at a low density of Y 
particles. In this model, which allows only one particle per site, the diffusion of 
a Y particle is permitted only if its neighbouring lattice sites are empty. Strictly 
speaking, therefore, to obtain the the rate of diffusion of the Y particles towards 
the foremost X particle one should also take into account probability of having two 
next-nearest neighbour Y particles. At low density of Y particles this probability 
can be neglected, permitting the use of Dy in Eq. (41). 
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2.6.5 Reaction- diffusion Model 2X with Only One Particle Allowed per 
Lattice Site 

In Sec. 2.6.5, we discuss the model of point (iv) of Sec. 2.6.1. Below, up to Eq. 
(46) we present the general features of this model following Refs. [14,54,55,90], 
and then we discuss the results in different parameter regimes. 

Like in the other subsections of Sec. 2.6, for all values of D, e and W , we define 
the position of the front by the location of the foremost particle. For a given 
realization, let us denote its position at time t by kf(t). For an ensemble of 
front realizations, having denoted the probability distribution for the foremost 
particle to be at lattice site kf by Pk f (t), its evolution is then described by 

= (D + e)P kf - 1 + [DPffi + WP^ +1 

- (D + e)P k} - [DP™ ipty + WP°' C ] . (42) 

Here P^ c (t) and i^™ pty (t) denote the joint probabilities that the foremost 
particle is at site kf and that the site kf — 1 is occupied or empty, respec- 
tively. Clearly, P kf (t) = P° k f(t) + P™ pty (t), and £ fc/ P fc/ (t) = 1. The first 
term on the r.h.s. of Eq. (42) describes the increase in Pk f (t) due to the ad- 
vancement of a foremost occupied lattice site from position kf—1, while the 
second term describes the increase in Pk f (t) due to the retreat of a foremost 
occupied lattice site from position kf+l. The third and the fourth terms, re- 
spectively, describe the decrease in P kf (t) due to the advancement and retreat 
of a foremost occupied lattice site from position kf. 

From the definition of Pk f (t), the mean position and the width of the distri- 
bution for the positions of the foremost occupied lattice sites are defined as 
x(t) = Efe / k f P kf {t) and W 2 (t) = £ fc [ife/ - x(t)fP kf (t). 46 The mean speed 
and diffusion coefficient of the front are thus given in terms of these quantities 
as the i->oo limit of v N = dx(t)/dt and W 2 (t) = 2D f t — see Fig. 16(c). To 
obtain them, we need the expressions of P^ c (t) and P^° ipty (t). To start with, 
we have 

P^(t)=p kf - 1 Pk f (t), (43) 



where Pk f -i is the conditional probability of having the (kf— l)th lattice site 
occupied (the foremost particle is at the kfth lattice site). The set of con- 

46 One can also define the foremost occupied lattice site for a realization as the one 
on the right of which no lattice site has ever been occupied before, and obtain the 
front speed from this definition following Sec. 2.3. Both of them of course yield the 
same result due to the t — > oo limit. 
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ditional occupation densities Pk f - m for m > 1 can be thought of as deter- 
mining the front profile in a frame moving with each front realization. For 
obtaining and Df, we simply need to know the asymptotic long-time limit 
Pk f -i(t — > oo), which from here on we will denote simply as p kf -\. Given 
Pk f -i, it is then straightforward to obtain from Eq. (42) and the conditions 
P kf (t) = P^(t)+P^(t)MdE kf Pk f (t) = l 



dx 

Vn = ~dt = 6 ~ Pk f~^ W ~ D ^ and 
dW 2 

— = 2D + e + p kf -!(W - D) . (44) 

Of these, the second equation indicates that the front wandering is diffusive, 
and an expression of the front diffusion coefficient Df is therefore given by 



2D + e + p kf -!(W - D) 



(45) 



Note that if we use Eq. (43) in Eq. (42), the latter equation has the form of 
the master equation for a single random walker on a chain. Thus we can think 
of the foremost particle as executing a biased random walk, and Df as the 
effective diffusion coefficient of this walker. Moreover, if we eliminate p kf -i 
from Eqs. (44) and (45), we get the following exact relation 



v N /2 + D f = D + e. 



(46) 



With this background in mind for all values of D, e and W, we now explore 
the results in different parameter regimes. 

The mean-field limit of this model effectively yields some form of Fisher- 
Kolmogorov equation. The applicability of the mean- field result, however, is 
confined only to large values of the diffusion coefficient (while keeping the other 
parameters e and W finite). The need to have large diffusion coefficient is not 
very difficult to understand intuitively — notice that although in principle 
there are three parameters in the problem, effectively there are only two of 
them, since the third one simply sets the time scale. Let us choose these two 
parameters as D/e and D/W . The quantity = min (^JdJe, ^Jd/W^J [81] 
then can be seen to define a correlation volume, i.e., a length scale over which 
the position of a particle varies in a unit time. When £d oo, the correlations 
involving two or more particles within a finite range of lattice sites get washed 
out, and then a mean-field description of the problem becomes suitable. 

Having defined x = k/y/2D, for W = and large values of D/e, a math- 
ematical derivation of the mean-field limit of this model in the form of the 
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Fisher-Kolmogorov equation 



+ (47) 

can be found in Ref. [14] (also see the references cited therein). 47 While we 
leave the interested readers on their own for this mathematical argument, 
here we concentrate on the Refs. [54, 55] that investigated, by means of a 
"self-consistent two particle model", the role of interparticle correlations in 
the approach of the front speed towards its mean- field value 48 v* = l^pllTe 
for increasing value of D/e. We make an important observation here that this 
approach only works for W = (we will return to this point later). 

The first step of this "self-consistent two particle model" is to derive the front 
speed for small D/e at W = 0, and it appears already in [14]: the idea is 
that far behind the foremost particle, the density of particles approach the 
homogeneous equilibrium density N: lim™-^ pk f - m = N. From the master 
equation of the microscopic processes for W — 0, it is easy to show that the 
homogeneous equilibrium solution for the total probability is of product form 
(so that the probability of having different sites occupied is uncorrelated) , and 
that the equilibrium occupation density iV is simply given by iV = 1. 49 Then 
the crudest approximation for the front profile Pk f - m and in particular for 
m — 1 is just to take Pk f -i ~ N. Substitution of this approximation into Eq. 
(44) then immediately yields [14, 54, 55] 



VN = £ + D . (48) 

Equation (47) agrees well with the simulations in the limit D/e — > 0. 

While Ref. [55] is mostly about numerical studies of how the mean-field result 
for front speed v* = 2\/2De is obtained from the first half of Eq. (44) in 
the limit D/e — > oo, the transition from the front speed in Eq. (47) to v* 

In a recent paper [81], the approach to the mean-field description of the front 
has been investigated in the light of Eq. (11). The idea is to define a quantity 
N* = V De/(e + W) such that for N* — > oo (which happens for finite N and 
Id — ► oo), not only does one recover the Fisher-Kolmogorov equation, but also for 
N* » 1, one expects to observe v* — vn ~ 1/ln 2 N* and Df ~ 1/ln 3 N*. In other 
words, Ref. [81] has essentially shown that iV* simply plays the same role as N does 
in Sees. 2.2-2.5. 

48 The mean-field Fisher-Kolmogorov equation (47) has front speed lyfe. When the 
distance is measured in terms of k and not x, this means that the actual v* is 
2 V / 2~D£. 

49 For W 7^ 0, N becomes e/{e + W) [90]; we will return to this case later in Sec. 
2.6.5. 
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for increasing values of D/e in analyzed by theoretical means in Ref. [54]. 
What Ref. [54] does is to develop a theory for evaluating Pk f -i, by means of 
considering the probability distribution pj of having j empty sites between the 
foremost particle and the second foremost particle, and Pk f -i is then simply 
obtained as Pk f -i = Po- The task of obtaining p itself, however, is not easy, and 
it is obtained from solving the master equation for the steady state solution 
of {pj}, as we outline below. 

With the usage of the normalization J2jPj — 1, the master equation for pj is 
given by 

pj = 2Dp j ^ 1 + (2D + e)p j+ i - (AD + 3s) pj for j > and 

Po = (2D + e)p l + 2e - 2(D + e)p . (49) 

Equation (49) is obtained from the following considerations: the j > state 
can be obtained from the j — 1 state by a forward diffusion of the foremost 
particle, or by the backward diffusion of the second foremost particle, giving 
rise to the 2Dpj_i term. The j > state can also be obtained from (j + 1) state 
by a backward diffusion of the foremost particle, or the forward diffusion of 
the second foremost particle, or by the growth of another particle just on the 
right of the second foremost particle — this gives rise to the (2D + e)pj + \ term. 
Any j > state contributes not only to j + 1 and j — 1 states by the reverse 
processes, but also to the j = state when a particle is created on either 
side of the foremost particle, which finally explains the (AD + 3e)pj term. The 
(2D + e)pi is obtained from any j > state by the creation of a particle on 
either side of the foremost particle. In addition, the j = state from the j = 1 
state by a creation of a particle on the right of the second foremost particle, by 
a forward diffusion of the second foremost particle or by a backward diffusion 
of the foremost particle. On the other hand, state j = contributes only 
to j = 1 state when the second foremost particle diffuses backward or the 
foremost particle diffuses forward. When all these are properly taken into 
account along with the normalization condition J2jPj — 1? the second half of 
Eq. (48) is obtained immediately. The steady state of Eq. (48) is then easily 
solved by the ansatz pj — p (l — p o y for Vj > [54], yielding 



V9e 2 + 16De - e 
Po = 2(2D + e) ■ (50) 



The lesson one learns from this clever approach is actually quite instructive. 
Notice that Eq. (49) automatically yields the proper D/e — > limit of the 
front speed (48). Its D/e — > oo behaviour however yields = V De, which 
is a factor of \/8 smaller than v* = 2\j2De. This itself should not strike one 
as surprising. Instead, with the observation that the involvement of a third 
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foremost particle on the left of the second one is always ignored (in the sense 
that the site on the left of the second foremost particle is always assumed to be 
empty), the fact that even in this crude form this formalism at least provides 
the \/De behaviour of vn at large D/e is surprising. 
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Fig. 22. vn/v* for e = 2 and 7 = 2D: simulations (crosses), "correlated two-particle 
self-consistent approach" (solid curve) , and "uncorrelated n-particle approach" with 
n = 3 (open triangles), 5 (open squares) and 10 (open circles). Also shown is the 
ad hoc functional dependence vn = [1 + 7(1 + 7/8)~ 1 / 2 ]/2 of vn/v* on 7 by the 
dotted curve. 



For a remedy of this drawback, Ref. [54] also invokes a "correlated two-particle 
self-consistent approach", in which diffusive moves of the second foremost 
particle are prevented due to the presence of the third foremost particle just 
next to the second one. With p a as the probability of having the second and 
the third foremost particles next to each other, the probability of the second 
foremost particle making a diffusive hop towards the left is reduced by D(l — 
p a ). Reference [54] then obtains the front speed as a function of increasing 
D/e by means of an effective functional dependence of p a on p , namely p a = 
Bp — (B — 1)pq, where the parameter B is chosen in such a way that it leads 
to vn — > v* for D/e — > 00. However, in this process, an analytical form of 
vn as a function of D/e is not obtained, although a predicted value of vn 
can be numerically obtained for a given value of D/e. In Fig. 22, we present 
the graph from Ref. [54] for e = 2 and 7 = 2D. It shows all the results for 
vn/v* together: simulations (crosses), "correlated two-particle self-consistent 
approach" (solid curve), and "uncorrelated n-particle approach" with n = 3 
(open triangles), 5 (open squares) and 10 (open circles). Also shown in Fig. 22 
is the ad hoc functional dependence vn — [1 + 7(1 + 7/8) -1 / 2 ]/2 of vn/v* on 
7 (adopted because of its simple functional dependence having proper limits 



55 



at D/e — > and D/e — > oo) by the dotted curve. 

Having thus concluded the VF = cases studied in the literature in the above 
paragraph, we now embark on the W ^ case for this model. The first case 
of W 7^ 0, namely W = D, has been analyzed in Refs. [5] and [34]. Notice 
from Eq. (44) that with W = D, one does not require to evaluate Pk f -i any 
longer — the front speed and front diffusion are simply given by v n = £ and 
D f = D + e/2 [5,34]. 

Clearly enough, when W^O but W ^ D, one needs to evaluate Pk f -i- Refer- 
ence [90] considered this case, but only with all parameters of the same order 
of magnitude. In this case, having estimated Pk f -i hj N = ]im m _> 00 Pk f - m (see 
footnote 49), one obtains [90] 



v N 



e + W 



2(e + W) 



(51) 



The comparison of Eq. (51) with stochastic simulation data for D = e = 0.25 
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Fig. 23. Comparison of Eq. (51) with stochastic simulation data for D = e = 0.25 
are presented in Fig. 23 as a function of W for D = e = 0.25. The larger open circles 
in the two graphs denote the results of D = W case analyzed in Refs. [5,34]. 

are presented in Fig. 23 as a function of W for D = e = 0.25. The larger open 
circles in the two graphs of Fig. 23 denote the results of D = W case analyzed 
in Refs. [5,34]. Notice that close to W — and at larger values of W , the 
agreement between Eq. (51) with stochastic simulation data is not good. It is 
in fact caused by the errors incurred in estimating Pk f -i by N, as shown in 



Fig. 24 by means of the relative deviation d 
lattice sites behind the foremost particle. 



(pk — N)/N between for a few 



In principle, for W ^ 0, it should be possible to perform an analysis in the 
spirit of Refs. [54,55] to get successively more accurate expressions for Pk f -i, 
and correspondingly for the front speed and diffusion coefficient. In particular, 
such extensions might allow one to use the results in a wider parameter range, 
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Fig. 24. Relative deviation d = (pu — N)/N of the average density from 
N = e/(e + W) for the first six lattice sites to the left of the foremost occupied 
lattice site, kf, for D = e = 0.25 and three different values of W. 

such as D/W -> oo while D/e ~ (9(1), or D/e -> oo while W/e ~ 0(1). 
However, inspection of the earlier analysis suggests that such higher order an- 
alytical expressions of Pk f -i are less trivial to obtain than one might expect 
at first sight. More precisely, for W = 0, the master equation for the proba- 
bility pj that the two foremost particles are separated by j lattice sites can 
be closed in a simple manner [54,55]. For W = formalism, no particle gets 
annihilated, and as a result, the hierarchy of equations for the joint probabil- 
ity density distribution of the two foremost particles can be closed easily at 
the simplest level, since in the absence of annihilation, the third foremost par- 
ticle never becomes the second foremost particle ( "uncorrected two-particle 
approach"). At this level, the expression of Pk s -\ can then be analytically 
solved, leading to a better approximation than what we use in this paper for 
W = 0. Of course, the master equation can be closed at a higher level, by 
considering more than two foremost particles to determine Pk f -i, but then 
one does not obtain an analytical expression of Pk f -i ("correlated two-particle 
self-consistent approach"). As soon as W ^ 0, this is not true anymore: con- 
sider the following situation where the two foremost particles are next to each 
other. With annihilation of particles allowed, one of them can annihilate the 
other, and then the probability distribution function of the two foremost par- 
ticles is crucially coupled to those which involve particles further back at the 
simplest level. It is therefore clear that for W ^ 0, the master equation for 
Pj involves particles that are further back. While it is certainly possible to 
solve the master equation numerically, it does not appear to lead one to an 
analytical expression of Pk f -i that provides a better approximation than the 
one already used in Ref . [90] . 
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2.6.6 Fluctuating "Pulled" Fronts in Higher than One Spatial Dimensions 

Up until now in Sec. 2.6, we have only considered fronts on a one-dimensional 
lattice — and we have witnessed that with iV being the conditionally aver- 
aged number of particles per lattice site in the stable phase for propagating 
fronts into unstable states, the front propagates with a speed v^, and simulta- 
neously the front positions in different realizations of an ensemble undergoes 
diffusive wandering around its ensemble averaged position. In higher than one 
spatial dimensions, the situation gets more complicated, although the basics 
of front propagation in discrete particle and lattice systems remain unchanged 
- therein the interest is in the wandering properties of the front, not so much 
in its speed [80,96,116]. 
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Fig. 25. Microscopic processes in the <i-dimensional reaction-diffusion system 
X+X^ 2X along the k\ direction (Z = 1,2, ... ,d) for a given value of I. On any 
lattice site, at most one X particle is allowed at any time. The dynamics along any 
one of the ki directions is exactly the same as that of Fig. 16. 

Since all the studies on discrete particle and lattice system of propagating 
fronts in higher than one spatial dimensions have been on reaction-diffusion 
systems X+X^ 2X [80,96,116] or X+Y^ 2X [72], to understand the general 
idea, it is best to first follow the reaction-diffusion system X+X^ 2X of Ref. 
[96] on a d dimensional square lattice for various values of d and briefly discuss 
its results. The system consists of a long lattice in the ki direction (— L x < 
ki < Li) with a transverse length L = L 2 = L 3 = . . . = L d . Any lattice site 
can hold at most one X particle at any time. The microscopic dynamics of 
the particles along the fcj direction (7 = 1,2, ... ,d) for a given value of I is 
shown in Fig. 25 — the dynamics along any one of the k\ directions is exactly 
the same as that of Fig. 16. With periodic boundary conditions applied along 
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all the transverse directions &2, . . . , kd and the initial condition for the 
particle density (f>^(t = 0) = Q(ki)N, one studies front propagation in the 
fci-direction. 50 

From this perspective, it is conceivable that the front properties can be studied 
by simply projecting the quantities defined in d spatial dimensions on to the 
k\ direction. Precisely this theme was developed in Ref. [96]: a mean particle 
density p kl (t) was introduced by projecting the particle density profile in the d- 



£ =2 M*) 



/L d 



-1 51 



One 



dimensional space on to the h\ direction as p kl (t) 
of the natural consequences of having to study the front properties by means of 
this projection is that one cannot define a front position based on the foremost 
particle as we have done so far in Sees. 2.6.2-2.6.5. To bypass this problem, a 
second quantity Ki(t) = Y^kl=-Li Pki(t) /N — L\ was introduced in Ref. [96] 
for the front position for any single realization. While the ensemble average 
(Ki(t)) ~ visit at long times, it was proposed that the wandering properties of 
the front around its ensemble averaged front position (Ki(t)) as a function of 
time can be measured by tracking the time development of quantity w(t) [96], 
where [the angular brackets in Eq. (53) also denote ensemble averaging] 



w\t) = (l£ L h Pkl (t) - (^i(t)) 2 - . (53) 

With these propositions, Ref. [96] invoked the scaling hypothesis for w(t), i.e., 
w(t) ~ t a F(t/L^), and then continued further on to measure a and (3 from 
simulations: for d = 2, a = 0.272 ± 0.007 and f3 — 1; for d — 3, the data are 
indistinguishable between the two possibilities w(t) ~ t 01 and w(t) ~ Vint; 
and for d = 4, a = (all simulations were carried out with D = 0.5 and 
e/W = 0.1). 52 



50 Even in d spatial dimensions, N = e/(e+W), just like the one-dimensional model 
of Fig. 16. See footnotes 49 and 52. 

51 It is then clearly seen that p\- x (t) obeys the Fisher-Kolmogorov equation 

= D\p kl+1 + p kl -i - 2p kl ] 

+ [ Pkl + 1 + 2(d- l)p kl + p fcl _x] [e - (e + W)p kl ] . (52) 

Notice that the corresponding one-dimensional problem in Sec. 2.6.5 also yields Eq. 
(52) with d = 1 (although we never quite discussed it in Sec. 2.6.5). From Eq. (52), 
we find that N = e/(e + W) for any d, just as we found in Eq. (32). 

52 At this juncture we briefly discuss the results of Ref. [72]. This paper consid- 
ers front propagation in a three-dimensional reaction-diffusion system X+Y^ 2X, 
where initially the whole lattice is full of Y particles with an equilibrium concen- 
tration N, and at time t = one X particle is introduced at one location. Both 
the species diffuse with equal rate D on the three-dimensional nearest neighbour 
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It turns out however, that for the reaction-diffusion model X+X^ 2X of Fig. 
25, invoking the scaling relation to measure the front wandering is an ex- 
tremely tricky issue. In fact, although w{t) in Eq. (53) is a measure of the 
front wandering, the full measure of wandering of the (d — l)-dimensional in- 
terface for fronts in d spatial dimensions is already lost when the projection 
Pkj{t) of <p%{t) along the k\ direction is taken before the front wandering is 
defined. Truly speaking, if 4>%{t) is not reduced to Pkj.it) before the scaling 
relation is invoked, then on a rf-dimensional lattice, the problem of front wan- 
dering can be related to interfacial growth and associated scaling analyses of 
the interfacial roughness [10,47,53,60,61]. The scaling idea in these interfa- 
cial growth phenomena is the following: for interfacial growth in d dimensions 
spanned by (k\, &2, . . . , kd) co-ordinates [leading to a (d — l)-dimensional in- 
terface spanned by {k 2 , k 3 , . . . , kd) co-ordinates], the "height" of the interface 
in the k\ direction is described by hk 2 ,k 3 ,...,k d {t) = h{t) at time t. With overline 
denoting the projection of a d- dimensional quantity on to one (along the k\ 
direction) dimension [e.g., Pkjit) = <Pj:^}, and the angular brackets denoting 
ensemble average, the actual scaling relation is given by W(t) = t a F(t/L^), 

where W 2 {t) = {[h{t) — h(t)} 2 ). When the argument 77 of F(rj) is <C 1, then 
F(rj) behaves like a constant; on the other hand, for 77 ^> 1, F(i]) ~ r]~ a . This 
implies that for t ^> LP, W(t) saturates at W sa t ~ L 1 , where 7 = afl. 

Even with this modification, the actual values of a and (3 in different dimen- 
sions for the X+X^ 2X reaction-diffusion model of Fig. 25 are not free from 
ambiguity [80, 116]. 53 The ambiguity stems from how hit) is related to (f>%{t) 
[although it is clear that one should not use p^it) for hit) to draw analogy 
between the front wandering properties and the interface growth roughness 
scaling relations on a d- dimensional lattice]. To illustrate this point, let us 
consider the example case for d — 2. It turns out that when hit) equals a 
simple local coarse-graining over the field values <j)%(t) within a cubic box of 
size {21 + l) d around k, then for £ — 1, one obtains a = 0.29 ± 0.01 and 
7 ~ 0.4 ± 0.02 [116] (a = 0.27 ± 0.01 and 7 ~ 0.41 ± 0.02 [80]). Interestingly 
enough, these exponents are surprisingly close to the KPZ exponents for a 



lattice sites, and any Y particle coming in contact with any other X particle on the 
same lattice site is instantaneously and irreversibly converted to an X particle - 
this is the d-dimensional generalization of the one-dimensional lattice model of Sec. 
2.6.2. The object of study was to see the deviations of vn from the prediction of 
the Fisher-Kolmogorov equation, v* . The numerical finding was that deviates 
from v*. Not only that there was no study of the broadening of the front region by 
means of measuring w(t) therein, but also the front speed as a function of the total 
average particle density was not carefully fitted with any theoretical curve. In this 
context, we refer to Fig. 2 of Ref. [80], which shows, for the reaction-diffusion model 
X+X^ 2X of Fig. 25, that for e/D — > 0, front speed in d > 1 does behave as 
y/e as the Fisher-Kolmogorov equation would predict. 
53 We note here that only W = D case was considered in Ref. [80] . 
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two-dimensional interface [116] (a = 0.245 ± 0.003, 7 = 0.393 ± 0.003 [77]). 
However, it has been shown that with increasing values of £, i.e., £ = 2,3,4 
etc., the exponents a and 7 approach the KPZ exponents for a one-dimensional 
interface (a = 1/3, 7 = 0.5) [80]. So the general picture that emerges from 
these analyses is that as one looks at progressively longer length scales, the 
wandering properties of the front makes a transition from a non-KPZ to a 
KPZ behaviour [80]. 

With £ = 4, the KPZ scaling behaviour for the front wandering on a three 
( 7 = 0.393 [77]) and a four dimensional lattice (7 = 0.313 [77]) for the X+X^ 
2X reaction-diffusion model of Fig. 25 has been recovered [80]. These results 
together indicate that in this model, any generic upper critical dimension has 
to be higher than 4. More simulations are however needed to conclusively 
clarify this point. 



2. 7 Convergence of the Asymptotic Front Speed to for Fluctuating Pushed 
Fronts as N — > 00 



We had earlier pointed out that at the time of writing this review article, I have 
not seen any work on how vjy behaves as a function of N in the limit iV — > 00. 
Nevertheless, in Sec. 2.7, we will show that if one incorporated the discrete 
nature of the particles by means of an effective cutoff 54 following the ideas of 
Brunet and Derrida as described in Sec. 2.2 [19], then the corresponding front 
speed Vj\[ must behave as v * — v n ~ iV~ 7 , where 7 is a positive number (c.f. 
footnote 8). 

The route we follow to demonstrate this power law convergence of vn to v' is 
by means of considering an example of a front propagation from a stable state 
to a (meta)stable state in [56] 



54 The introduction of an effective cutoff to model the discreteness of particles, even 
for fluctuating pushed fronts, is not unusual. Consider for example the following 
deterministic equation in continuous space and time [12, 100] 

^flg + ^d-'-iirry (54) 

for b > 0. This system then has a linearly unstable state at (f> = 0, and a stable state 
at <j) = 1. In a discrete particle and lattice model of this model with <pk = ((Nk))/N, 
where N is the (conditionally) average number of particles per lattice site at the 
stable phase of the front, at the tip, the value of 4> is once again of 0(1/N). The 
effect of discreteness of particles on the front, once again, can then be mimicked by 
putting a growth cutoff for (j> at s ~ 1/N. 
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f = ^ + (l-0 2 )(0 + a), (55) 

where < a < 1. This equation, which appears with the name of Ginzburg- 
Landau equation or with the so-called Schlogl model, has two stable states, 
= ±1, and an unstable state at = —a. For a steep initial condition that 
connects = 1 to = — 1, a front propagates from the = 1 to the = — 1 
state with speed = ay/2; and the corresponding front profile in terms of 
C = x — v H is described by [56] 



tanh 



(56) 



such that A f = [u* + yV 2 + 8D(1 - a)]/{2D). 55 Notice that the Ginzburg- 
landau equation (55) does not fall in the category of fronts propagating into 
unstable state, the procedure to obtain V s — v e ~ |e|~ 7 in Eq. (55) is very 
instructive, and it can be easily applied to pushed fronts propagating into 
unstable states. 56 



As we add a growth cutoff in the Ginzburg-Landau equation (55) at = —1+e, 
we can use the technique developed by Brunet and Derrida (c.f. Sec. 2.2 [19]) 
to solve for the asymptotic front speed v £ [56]. Let us denote the location, in 
the comoving co-ordinate ( = x — v £ t, where the value of reaches —1+e 
by Co- For e — > 0, Eq. (55) can clearly be linearized on the right of Co- in 
addition, in the left neighbourhood of ( , the value of is still infinitesimally 
close to —1, which allows one to linearize Eq. (55) around the unstable state 
once again. This yields the following form of the front solution [56] 57 



55 It is not surprising at all that the front solutions in Eqs. (54) and (55) bear strong 
similarities [29] — with appropriate redefinition of and rescaling of the diffusion 
coefficient D and time t, one equation can be converted into another. However, 
one has to be careful here: despite the convertibility of one form to another, the 
difference between the nature of the fronts in Eqs. (54) and (55) stems from the 
fact that b > in Eq. (55) corresponds to a > 1 in Eq. (54), for which the front 
propagation is into a linearly unstable state. On the other hand, < a < 1 in Eq. 
(54) corresponds to b < in Eq. (55) — in this case, front propagates from a stable 
to an (meta)stable state. 

56 Notice that with = ay/2, the front profile (56) is directly obtained by solving 
Eq. (55). In general, in the comoving frame for any given asymptotic front speed 
Vas, upon linearizing Eq. (55) around = — 1 and using the front profile ansatz 
0(C) ~ — 1 + Ae~ xc - for ( — > oo, two roots of the exponent A are obtained: Ai = 
K + y/vl + 8D(1 - a)]/ (2D) > and A 2 = [ Vas - yfig + 8D(1 - a)]/ (2D) < 0. 
In the corresponding front solution 0(C) ~ —1 + ^4ie~ Al '' + ^e^ 21 *, A2 must be 
zero for the solution to be physically relevant. This condition is clearly obeyed by 
the front profile (56). 

57 The reason for expressing the coefficient of e~ A2 "> in Eq. (57) in this manner is 
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= -1 + A 1 e- x ^ + A 2 ( v t _ Ve ) e - A k for large C, but C < Co 
= -1 + ee- v ^-^ for C > Co , (57) 

where the expressions of \\ and A2 are obtained from footnote 56 after having 
Vas substituted by i>t. As the two solutions of Eq. (57) and their derivatives 
are matched at Co (with the condition that A± and A 2 are independent of e; 
see footnote 57), it is easily seen that Co ~ I hie| and [56] 58 



£ i-xl/x{ _ 



(58) 



1/2, \ 2 /\\ = -1/2, and v 



t _ 



e 3 / 2 . This scaling of — v £ 



For D = ( 

has been observed in the simulation [56], and it is shown in Fig. 26. 



10" 




Fig. 26. The log-log curve of 5v = — i> £ vs. £ for the Ginzburg-Landau equation 
(55) with D = a = 1/2 [56]. Solid line: simulation data, dashed line: visual aid for 
the Sv ~ e 3 / 2 behaviour predicted in Eq. (58). 



As already pointed out, the exercise can be repeated to obtain a power law 
behaviour of v> — v £ for the Eq. (54), where the front does propagate into 
an unstable state. Based on Brunet and Derrida's insight that the effect of 

actually quite subtle. The idea lies in the fact that just like we linearize Eq. (55) 
with a growth cutoff around —1 for large £ to obtain Eq. (57), we can also linearize 
it around 1 for £ — > —00. This linearization suggests an equivalent form of the 

front solution <^°)(C) = 1 — i?ie A i'^ — B 2 {v^ — v £ )e X ^'^ for ( — ► —00, where Aj, = 
[-■yt + vV 2 + 8D(1 + a)]/ (2D) > and \\, = bV - y^t 2 + 8L>(1 + a)]/ (2D) < 0. 
The linearized equations at £ — > ±00 involve v e as an input parameter. However, in 
this form of (f>^°\C) at £ — > ±00, Ai, A 2 , By and S2' all become independent of e. 
58 Equation (58) has also been obtained by Chomaz et al [25], in a study of the 
influence of a fixed boundary in a system with a convective instability. 
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discreteness in terms of particles and the lattice is well captured by the choice 
e~ 1/iV as iV — * oo, the above analysis indicates that in a discrete particle 
and lattice model of fluctuating pushed fronts v * — v N would indeed behave 
as A" 7 for some 7 > 0. However, at the time of writing this review article, I 
have not seen any empirical verification of this prediction at this point. 



2.8 Epilogue I 

As we discussed front propagation in discrete particle and (mostly one-dime- 
nsional) lattice systems throughout Sec. 2, we started with an asymptotically 
large value of A in Sec. 2.2, and from there onwards we gradually reduced 
the value of A till Sec. 2.6. We have witnessed that pulled fronts in discrete 
particle and lattice systems do not exist, at best they are weakly pushed. With 
decreasing value of A, the weakly pushed nature of these fronts becomes more 
and more prominent, and the front properties start to depend heavily on 
the details of the models under consideration. Because of these complications 
at intermediate values of A, in general, it is very difficult to obtain a first- 
principle based predictive expression of vjy for fluctuating "pulled" fronts. In 
a contradictory trend to this, however, at the limit of iV when there is at 
most one particle allowed per lattice site, front speed and the front diffusion 
coefficient become tractable again (c.f. Sec. 2.6.5), giving rise to exact relation 
like Eq. (46), or approximate theoretical expressions like Eqs. (46) and (51). 

Even if we exclude the small or moderately large (but not asymptotically 
large!) values of A from these concluding remarks on Sec. 2, so many different 
features of discrete particle and lattice systems of fronts propagating into 
unstable states for iV — > 00 discussed in Sees. 2.2, 2.4-2.5 and 2.7 may at 
first seem a bit puzzling to the reader. A moment's reflection, however, shows 
that the origin of such versatile phenomena can be intuitively understood 
in a unified perspective — to be more precise, from the stability spectrum 
of (deterministic) pulled and pushed fronts propagating into unstable states, 
already briefly discussed at the end of Sec. 1.1: the gap less stability spectrum 
of pulled fronts leads to its power law convergence (4) of the front speed v(t) 
to v* in time for the asymptotic front speed [38, 101, 102] has already been 
noted in Sec. 1.1. On the other hand, stability spectrum of pushed fronts is 
gapped [22,23,37,56,92,101,102], which leads to exponential convergence of 
the front speed v(t) to [56, 102]. 

What we observe in fluctuating "pulled" fronts and as well as in fluctuating 
pushed fronts for A — > 00 is intimately connected to their stability spectra. 
For pulled fronts, the gapless nature of the stability spectrum is indeed another 
representation of the linear marginal stability criterion. It therefore must come 
as no surprise that any change in the front dynamics close to the linearly 
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unstable state must affect the front properties severely. In Sees. 2.2 and 2.5 
this is manifested by the very slow convergence of and Df to v* and zero 
respectively, while the stability spectrum of fluctuating "pulled" fronts itself 
becomes gapped with a very small gap oc 1/ In 2 N (c.f. Sec. 2.4). On the other 
hand, for pushed fronts, the spectrum is gapped, and any change in the front 
dynamics close to the linearly unstable [or (meta) stable] state is expected to 
affect the front properties minimally. This is then reflected in the stronger 
power law (in N) convergence of to v^. 59 

We finally close Sec. 2 on discrete particle and lattice systems of fluctuat- 
ing fronts propagating into unstable states by pointing out a few problems, 
whose solutions are not yet completely settled. These are: (a) for fluctuating 
"pulled" fronts a first principle based predictive theory for the front speed vn 
at intermediate values of N, (b) for the model of Sec. 2.6.3, why one observes 
mean-field Fisher-type behaviour even at a reasonably small value of N, (c) 
a better analytical prediction for the model of Sec. 2.6.5 for a larger set of 
parameter values, and (d) possible existence of an upper critical dimension for 
fluctuating "pulled" fronts, when there are at most one particle allowed per 
lattice site (c.f. Sec. 2.6.6). 



3 Field-theory of Fluctuating Fronts: External Fluctuations 

3. 1 External Fluctuations, Multiplicative Noise and Novikov 's Theorem 

We have seen earlier in Sec. 1.1 how propagating fronts have been observed in 
many systems related to physics, chemistry and biology. In all of these cases, 
the evolution equations of the propagating fronts invariably involve some pa- 
rameters that relate to experimental situations, such as parameters in chemical 
or biological reaction kinetics, an externally imposed fields etc. When these 
parameters correspond to externally imposed fields, then the variation in these 
fields acts as a varying external influence on the system — such as an exter- 
nally imposed electric or magnetic field, or light intensity in a photosensitive 
reaction. Similarly, when these parameters originate within the system itself, 
then usually they are law-of-mass-action-based mesoscopic mean-field theory 
estimates of quantities that originate in the underlying microscopic dynam- 
ics. For the most part, (Langevin-type) field-theoretical studies on the effects 
of fluctuations in systems with propagating fronts have been motivated by 
these phenomenological considerations. Other significant motivating factors 
for the field-theoretical approaches have been phase transitions induced by 

59 Let us just state here that the theoretical arguments indicate that the diffusion 
coefficient of fluctuating pushed fronts made of discrete particles on a lattice should 
scale as 1/N; however, for now we will leave this issue until Sec. 4.2.1. 
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noise [11,18,44,93], and the discovery that surface roughness in various growth 
processes fall under a handful number of universality classes [10,40,53,60,61]. 



In my opinion, for systems admitting propagating front solutions, external 
fluctuations are easier to handle than the internal ones — not only will we 
find ample evidences of this later on, but also two simple examples at this 
stage will make this point clear. Consider a front dynamics under the influ- 
ence of an external electric field E = Ex that couples linearly to the front 
field 0, also propagating in the x-direction. Fluctuations in the electric field 
of magnitude SE(x, t) around its mean value E^ simply gives rise to a term 
oc 5E(x,t)<f)(x,t) in the evolution equation of the front. Contrast this situa- 
tion with a chemical reaction front: therein the deterministic front evolution 
equation is already a mesoscopic (mean-field) description of the underlying 
microscopic process — these mesoscopic descriptions collectively give rise to 
the very definition of the front field 0, and as well as to the mass-action based 
parameters in the dynamics such as the reaction rate. When fluctuations are 
to be considered in such systems, then in general, fluctuations in both the 
number of particles per correlation volume (which usually defines the front 
field) and as well as in the parameters are to be considered together. In the 
full description of these fluctuations, naturally, complicated correlations be- 
tween the fluctuating front field and the parameters are unavoidable. To get a 
reasonable grip on the theoretical handling of these processes, the mass-action 
based mesoscopic parameters for the internal dynamics are often kept intact 
in the dynamics of the system (we have seen a large number of examples of it 
in Sec. 2) and separate fluctuation terms are added by hand to the evolution 
equations to mimic the internal fluctuations. 

Precisely because of these subtleties, field theories for the internal fluctuations 
often prove to be very tricky. For now, we will leave it for Sec. 4. As for Sec. 
3, we will devote it to the field-theory of external fluctuations on propagating 
fronts. 

To begin with, let us consider the deterministic equation for a scalar front in 
one spatial dimension coupled to n external fields Oj, i = 1,2, ... ,n: 



d(j) 
~dt 



F({ai},(f)) . 



(59) 



The effect of the fluctuations in the fields 
a- -* on the front properties is then given by 



around their mean-field values 



,(») 



(x,t) 



,(°) 



(60) 
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5 k T 



is a function of the front field. 

,(°) 



In Eq. (60), the functional derivative . , 

baf 

What is the precise form of this function depends on how the fields are coupled 
to the front dynamics, but the fact remains that in the evolution equation of 
the front, the coefficients of the fluctuations in the external fields are multiplied 
by functions of the front field — Eq. (60) is thus an equation with multiplicative 
noise. 

Before we proceed further, it is however necessary to first properly interpret the 
multiplicative noise term in Eq. (60). This brings us to the well-documented 
It 6 vs. Stratonovich dilemma for multiplicative noise [45,52]. In this review 
article, we will not delve deep into the details of which interpretation is the 
right one, but it is extremely important to be clear about the fact that throughout 
Sec. 3, the noise term in Eq. (60) will be interpreted in the Stratonovich sense. 
Section 3, in fact, is a collection of works done on the field theory of external 
fluctuations on propagating fronts for the last two decades, 60 and all these 
works interpreted the noise term in the Stratonovich sense. 

We will now see in Sec. 3.1.1 that in the context of front propagation, the 
so-called Novikov's theorem [85] has proved to be a very useful tool for the- 
oretical analysis with Stratonovich interpretation of multiplicative noise. It 
should however be kept in mind that under general circumstances, Novikov's 
theorem is not an automatic gateway to success in these problems. First of 
all, the theorem can handle noise only with Gaussian statistics. Secondly, it is 
practically useful only if the spatial and temporal part of the noise factorize, 
and in the limit when both the noise correlation length and correlation time 
go to zero. 



3.1.1 Novikov's Theorem in the Context of Front Propagation in Reaction- 
Diffusion Systems 

Consider the stochastic differential equation describing a general reaction- 
diffusion process coupled to a single external field 61 



60 Field-theoretical studies on the effect of external fluctuations (as multiplica- 
tive noise) on propagating front in the context of reaction-diffusion systems date 
back to the early 1980s [79, 107]. Some threads of it were picked up in the early 
1990s [94, 108]; but all of these works considered only weak noise, i.e., the case of 

[ai — a- \/af ] \ <C 1 in Eq. (60). The late 1990s till now witnessed a surge of activ- 
ities in this field [3,4,97-99,105,106] — they considered multiplicative noise that is 
not necessarily weak. 

61 Novikov's theorem can be easily generalized to arbitrary spatial dimensions and as 
well as to the case when the front field is a vector field, see Ref. [99]. Generalizations 
to more than one external field coupled to the front dynamics is also trivial. 
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d(f> & 



D- 



l 4 + f(<f>)+e 1 / 2 g( ( f>)r } (x,t), (61) 



dt dx 2 

where er)(x,t) is the space-time dependent fluctuation in the external field, e 
is the strength of the (multiplicative) noise, and r](x,t) has zero mean with 
a correlation time r c and correlation length A c . The factorized temporal and 
spatial correlations are expressed by the functions C t and C x respectively, such 
that (the angular brackets with a subscript rj denote averaging over the noise) 



(ri(x,t)ri(x',lf)) ri = 2C t (\t-t'\;r c ) C x (\x - x'\; X c ) and 

dx'C x (\x-x'\;X c ) = J dt'C t {\t-t'\-T c ) = 1. (62) 



/ 



The difficulty with using Eq. (61) directly is that although the mean of r](x, t) 
is zero, in the Stratonovich interpretation, the stochastic term in Eq. (61) 
[oc g((f>)r)(x,t)] does not have zero mean. This is precisely where Novikov's 
theorem comes in — its purpose is to (non-perturbatively) rewrite Eq. (61) in 
such a manner that the fluctuating term in it does have zero mean. Novikov's 
theorem, however, is applicable only if rj(x,t) itself is a Gaussian random 
function with zero mean [85]. 

The exact statement for Novikov's theorem is the following: Suppose h(s) 
is a Gaussian random function with zero mean value and correlation 

(h(s)h(s')) h = H(s,s'), (63) 

then for any functional H({h}), the following result holds [85]: 

(h(s)H({h})) h = J ds'H(s, s>) (^^y) h ■ (64) 

Equation (64) does not involve any approximations. However, to make any 
practical use of it for any analytical calculation in the present context, Nov- 
ikov's theorem is useful only if the correlation time r c and correlation length 
A c become much smaller than the respective time and length scales of the 
propagating front 62 — for example when r c is much smaller than the inverse 
reaction rate and A c is much smaller than the front width. In these cases, one 
can use both r c and A c — > and reduce Eq. (62) to 

(rj(x, t)r){x', t')) v = 25(t - t') S(x - x') . (65) 



62 For finite values of r c and A c , and for non-Gaussian forms of n(x, t) [and even 
for non-factorizable spatiotemporal correlations, unlike Eq. (62)], it is possible to 
formulate a Novikov-type theorem by carrying out a cumulant expansion of the 
noise. Sec Rcf. [106]. 
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The limit of both r c and A c — > 0, however, causes certain intricacies in the 
application of Novikov's theorem, as we now illustrate below. From Eqs. (63- 
65) we have 



oo oo 



(g[^x,t)]r ] (x,t)) v = 2 J dt' J dx C t (\t - t'\;r c ) C x (\x - x'\;X c ) x 



— oo — oo 



V 



After having solved <t>{x,t) from Eq. (61) as 



oo 

(j)(x,t) = J dt" 



+ t")} + e^gMx, t")] V (x, t") 



x 



x 6(t - t") , (67) 



the expression ^ ' ) in Eq. (66) reads 
dr)(x,t') 



^4 = e 1 /^^, t')} e(t - t>) S(x - x>) , (68) 
which subsequently implies that 



(g[ ( j ) (x,t)}r ] (x,t)} ri = 2e^ 2 C x (0;X c ) J dt' C t (\t - t'\;r c ) x 

— oo 

x (g>[ ( f ) (x,t)}g[ ( f)(x,t')}) J] Q(t-t>). (69) 

One can now clearly see that the replacement of C t (\t — t'\; r c ) by S(t — t') 
in the limit of r c — > results in a perfectly well behaved r.h.s. of Eq. (69), 
but the spatial white noise limit in the case of A c — > does not, since the 
A c — > limit sends C x (0; A c ) to oo. In most papers on field-theory of front 
propagation with multiplicative (external) noise, the limit r c — > in Eq. (62) 
is often taken implicitly, while the explicit dependence on C x (x — x'\ A c ) is left 
intact. In this limit, Eq. (69) reduces to 

e 1 / 2 t)] v (x, t)) v = e {g'[4>{x, t)\ g[4>{x, t)\) , (70) 

where e = eC x (0; A c ) and 0(0) = 1/2. 

In order to avoid the difficulty associated with the remaining C x (0; A c ) in Eq. 
(70) approaching infinity in the limit of r c and A c — > 0, operationally (i.e., in a 
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computer simulation) one simply discretizes both space and time in units of Ax 
and At. These discretization units then act as natural ultraviolet cutoffs. So 
long as At is sufficiently small, the results are independent of At. On the other 
hand, when Ax is sufficiently small, the Gaussian representation C x (0; A c ) of 
5-function in space is replaced by a square representation of 5-function, such 
that A c = Ax, and one uses [3,97] 



C*(0;A c -0) = (71) 

That this is the right interpretation is confirmed by very good agreements 
between the theoretical and the simulation results. 63 

With Eq. (70), we can finally rewrite (61) as (see for example Ref. [3]) 



^ = D^ + fM+e 1/2 R(<t>,x,t), (72) 
where f g (<f>) = f{4>) + eg(<f>)g'(<f>). Notice now that 



R(cf>, x, t) = g((/>)r,(x, t) - e l ' 2 C x (0; A c ) g(<f>)g'(<f>) (73) 
has zero mean and the correlation property that 64 

(R(<f>,x,t)R(<f>,x',t')) v = (5[0(x,t)]»7(a:,Off[^O]^(^>O)^ + O(e 1/2 ). (74) 



63 Equations (66-69) are indeed a fancy way of deriving Eq. (70). Because these 
equations use Novikov's theorem, they are only applicable if rj is a Gaussian random 
variable. When r\ is not Gaussian, then there is a simpler way to derive Eq. (70) for 
the Stratonovich interpretation of Eq. (61): see for example Ref. [45] on the issue of 
the equivalence of Ito and Stratonovich integrals. 

64 We will see in Sec. 3.1.2 that the front speed is obtained from the equation 

dep d^(f> 

— - = D——^ + f g ((p), and the right parameter for the effect of the stochastic term 
at ox 1 

on the front speed is e and not i. In that sense, not only is it clear that the choice 
of Ax does affect the front speed, but also in the limit Ax — ► 0, Eq. (72) captures 
the entire effect of the external fluctuations on the front speed and not just up to 
some order in a perturbation expansion in powers of e. The diffusion coefficient of 
the so-called Goldstone mode Dq [94,109,110], however, involves the assumption 
that e € 1, so that one can neglect the 0(e 1 / 2 ) term in Eq. (74). The issue of Dq 
is however a bit more complicated — although the ^-dependence of the diffusion 
coefficient of the Goldstone mode is more subtle than a simple series expansion in 
powers of e, the fact remains that for e<1, the contribution of the 0{e 1 / 2 ) term 
in Eq. (74) indeed yields a higher order correction [3]. 
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3.1.2 Reaction- Diffusion Systems: The Effect of Multiplicative Noise on the 
Front Speed and the Diffusion Coefficient of the So-called Goldstone 
Mode 

From the discussion in Sec. 3.1.1, it is clear that in the context of front prop- 
agation, application of Novikov's theorem to rewrite Eq. (61) to (72) is not 
a series expansion in powers of e. Neither does Eq. (72) involve any approxi- 
mations. Instead, it is simply a way to separate the systematic part and the 
fluctuating part of the stochastic differential equation (61) — as we will now 
see, rewriting Eq. (61) in the form of Eq. (72) has the advantage that the 
systematic part gives rise to a steady shift in the front speed, while the fluctu- 
ating part gives rise to two other phenomena at two different time scales: the 
(random) displacement of the front from its uniformly translating position at 
long time scales, and shape fluctuation of the front around its instantaneous 
position at short time scales (see the illustration in Fig. 27). A very simple 
and instructive illustration of this appears in Refs. [79,94,97,107,108], and 
here we will follow it to make our point. 




Fig. 27. Snapshots of fronts (with multiplicative noise) in the stochastic differen- 

tial equation — = -— r + 6(1 — 6)\a(x, t) + 6] and in the deterministic differential 
at ox 2 



equation — = — — » + 6(1 — 6)[a^ + 6] at two different times [41. The (dotted) 
dt ox 1 

shape of a step at x = 10 denotes initial condition for both the deterministic and 
fluctuating fronts. The dashed lines correspond to the front shapes for the deter- 
ministic differential equation with = 0.3, and the solid lines correspond to the 
fronts in the stochastic differential equation, where a(x, t) fluctuates around 
with an amplitude e = 0.3 (Ax = 0.1, At = 10~ 3 ). The increased gap between 
the deterministic and the fluctuating front profiles at the later time is due to the 
increased front speed in the presence of the noise, as described in Eq. (76). Symbols 
have been changed from the original to keep consistency with the ones used here. 



The method of separating these two different time scales goes back to the early 
1980s [79,107], although it was carried out for weak noise and without the aid 
of Novikov's theorem. For weak noise one carries out a perturbative expansion 
in the noise strength e, but the method to extract the diffusion coefficient of 
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the so-called Goldstone mode arising out of the random displacement of the 
front from its uniformly translating position at long time scales, nevertheless, 
remains equally applicable. We will not derive the expressions for the diffusion 
coefficient of the so-called Goldstone mode for weak noise [79,94,107,108] here; 
instead we will summarize the key ideas of Refs. [79,94,107,108] in points (i)- 
(iii) at the end of Sec. 3.1.2. Enthusiastic readers are encouraged to consult 
these references on their own to further satisfy their curiosity. 

To analyze different phenomena taking place at different time scales, we first 
express in Eq. (72) as a combination of a fixed shape (p^ displaced by the 
amount X(t) from its uniformly translating mean position £ = x — v £ t and 
(time-dependent) fluctuations in the front shape 5<p around the instantaneous 
position of the front as 



<f>(x,t) = <j>W[£-X(t)] + 5^-X(t),t]. 



(75) 



It is precisely at this point where one first separates out the systematic and the 
fluctuating part of the front. As for the fluctuating part, a further separation 
of the effect of fluctuating term oc R(<f), x, t) in Eq. (72) at two different time 
scales also takes place in Eq. (75). The fluctuations at long time scale is coded 
in the (random) wandering X(t) of the Goldstone mode around £ = x — v e t, 
whereas the fluctuations at short time scales in the front shape is coded in the 
fluctuations in the front shape, expressed by <50[£ — X(t),t\. In this form, the 
shape as a function of its argument £ is obtained from the equation 



D 



0. 



(76) 



As one can now see, the effective deterministic part now has a front solution 
that propagates with a speed v £ — the important thing to notice is that v £ is 
different from the speed of the deterministic front obtained by dropping the 
oc rj(x,t) term of Eq. (61) altogether. 

Clearly, in this process, the displacement of the front X(t) from its uniformly 
translating mean position £ can be chosen in many different ways. One of the 
convenient ways is to obtain it through the so-called Goldstone mode. The 
idea is that just like in Eq. (15), if is a solution of Eq. (76), then so it 

(/>(°'(£ + a) for any a. This translational invariance of Eq. (76) indicates that 

50(0) 



G,R = 



is the right eigenvector of the linear stability operator C Ve with 



eigenvalue zero, where 



+ v £ 



d_ 



+ + 



= 0(0) 



(77) 
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The so-called Goldstone mode is simply the quantity <&g,r(Q [94,109,110], 
which also propagates with the speed v £ of the front solution </>^(£)- 



As we already saw in Sec. 2.4.2, 5>g,l(0 — e 



is the left eigenvector 



of C Ve with eigenvalue zero, and then based on the idea of the Goldstone mode, 
we now uniquely define X(t) from the relation [79, 94, 97, 107, 108] 65 



/ dCe^ D ^5^-X(t),t] = 0. 



(78) 



Thereafter, having used Eq. (75) in Eq. (72), at the lowest order of the mag- 
nitude of the front shape fluctuation, we have, on the instantaneous comoving 
frame of the Goldstone mode 



d(5(f)) 
dt 



C 



Ve 1 



+ X(t) 



(79) 



and then having left multiplied Eq. (79) with $g,l(0 an d integrating both 
sides w.r.t. £, we get 



X(t) 



_~l/2 



d^aAO R(<f>,U) 



d^oAO^aAO 



(80) 



The upshot of the entire exercise in terms of the usage of Novikov's theorem 
for propagating fronts is now clear — the stochastic differential equation (61) 
has a front speed v £ , which has to be calculated from Eq. (76). On the other 
hand, the front speed also has (time-dependent) random fluctuations around 

0]. 



v £ [notice that (X(t)) 



One measure of the wandering of the entire front that is ubiquitous in litera- 
ture on the external fluctuation effects on propagating fronts is the diffusion 
coefficient D G of the Goldstone mode. The quantity D G , as can be easily 
guessed, is borne out of the random displacement X(t) at long time scales; 
but in order to obtain an expression for Dq from Eq. (80), it is necessary that 
one replaces _R(0, £,t) by R((j)^ \^,t). This replacement makes R 5-correlated 
both in space and time in Eq. (74), and then the front speed fluctuation be- 
comes a Markov process originating from the fluctuating R term in Eq. (72). 



65 For fluctuating fronts in reaction-diffusion systems, Eq. (78) is the operational 
definition of how to extract the conditionally averaged front profile from an ensemble 
of front realizations. 
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Moreover, it also reduces the whole formalism of Langevin-type front evolu- 
tion equation (61) with multiplicative noise in the Stratonovich interpretation 
to Eq. (72) with additive and white noise, resulting in the further reduction of 
the fluctuating "Langevin force" R of variable magnitude in time to a force 
with a fixed magnitude. This reduction is absolutely necessary in order to even 
define Dq [and also Eq. (84) in this form] for Stratonovich noise. 66 

By using the Kubo formula for the diffusion coefficient as the time integral 
of the autocorrelation function of the fluctuation in the front speed, and with 
the aid of Eq. (70) and (74) [the 0{e 1/2 ) term in Eq. (74) is dropped], the 
diffusion coefficient D G of the Goldstone mode is obtained as 



D 



G 



oo 

/ 



oo 



(81) 



Before we end Sec. 3.1.2, it is worth spending some time on three important 
issues. First, only after having derived Eq. (81), we can appreciate footnotes 
64. The front speed v £ is obtained from Eq. (76), and to obtain v £ , we have 
not carried out any expansion in powers of e. On the other hand, for the 
diffusion coefficient Dq of the Goldstone mode, we have dropped the 0{e 1 ^ 2 ) 
term in Eq. (74). Despite that, the actual e-dependence of Dq is not simply 
the prefactor e on the r.h.s. of Eq. (81); in fact, there are other e-dependent 
terms already occurring in the expression of v £ that enter both the numerator 
and as well as the denominator of Eq. (81) [3]. 

Secondly, the nomenclature "diffusion coefficient of the Goldstone mode" may 
a priori appear to be a misnomer — after all, since Dq characterizes the 
wandering properties of the Goldstone mode, one would, quite rightfully, tend 
to think that D G should really be the front diffusion coefficient. Why the name 
"diffusion coefficient of the Goldstone mode" has been chosen in this review 
article is in a way related to my own personal choice of words. It turns out 
that Df for the front defined earlier is in fact conceptually totally different 
from the field-theory expression (81). While we leave further details on this 
point till Sec. 4.2, to differentiate these two quantities, we will continue to use 
the terminology "diffusion coefficient of the Goldstone mode" for Eq. (81). 

Thirdly, the front shape fluctuations around the instantaneous position of 
the front at fast time scales can also be analyzed by defining the (mutually 
ortho normal) shape fluctuation modes {^ m (0} i n the eigenspace of non-zero 
eigenvalues of the linearized operator C Ve as 



66 We will later see in Sec. 4 that when R((p, £, t) is interpreted in the ltd sense, Dq 
can be well-defined without replacing i?(0, £,i) by R(<p^ \^,t). See also Sec. 3.5. 
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8<f>(£,t) = J2 C m(t)^m, R (0, (82) 

Here ^„i,r(0 is the right eigenvector of C Ve with eigenvalue r" 1 , and \I/ m ,L(0 — 
e Ve £/ D ty m> R is the corresponding left eigenvector. In the appropriate context, 
e Ve ^^ 2D ^ m R is basically the same as ip m (^) encountered in Eq. (21). Needless 
to say, 



oo 

J d£*m,L(OSG,fl(0 = 0, Vm (83) 

— oo 

i.e., each of these modes is orthogonal to the Goldstone mode. The mode 
expansion (82) then easily leads one to 



c m (t) = - r- 1 c m + a 1 ' 2 J di * m , L (0 R{^\ f , t) (84) 

— oo 

and its corresponding solution 

t' oo 

c m (t + f) = e-^c m (t) + e 1/2 Jdt"J di e-^^^OR^i^ O . (85) 



-oo 



Notice that in Eqs. (82-85), we have first used Vm(0 = e^ /(2D) ^ mvR (£) to 
first convert C Ve to a Hermitian operator [see Eqs. (16) and (17)], of which 
VVn(£) is the eigenvector with eigenvalue r" 1 . We then implicitly made use of 
the completeness condition for {ip m (0} : 



EW0«') = *(£-0- (86) 

m^O 

This completeness condition also means that c m (t) at any time t is automati- 
cally determined from the following equation: 



oo 

C m (t) = J d^m, L (0H(U)- (87) 



We will not need Eqs. (82-87) in Sec. 3; however, they will come handy for 
Sec. 4.2. 
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We finally end Sec. 3.1 with the remark that to derive Eqs. (61-81) for weak 
noise (e <C 1) and for the cases when r](x,t) does not necessarily satisfy 
Eq. (65) [79,94,107,108], using small noise expansion formulation [45], (p^ 
comes out to be the same as <fi* or (ft (correspondingly with v £ = v* or v^), 
depending on whether the deterministic part of Eq. (61) gives rise to a pulled 
or a pushed front. This implies that in these approaches, while there is no 
steady modification in the front speed due to the presence of the noise, the 
only new phenomenon that the noise gives rise to is the diffusive wandering of 
the Goldstone mode. Nevertheless, there are a few interesting results derived 
in Refs. [79,94,107,108] for weak noise, and we summarize them in points (i) 
through (iii) below. 

(i) A derivation of the diffusion coefficient of the Goldstone mode for general 
noise correlations can be found in Refs. [79,107]. An application to derive an 
exact result for Dq for front propagation in bistable systems can be found 

\-(t -t'^ 2 ' 

in Ref. [79] when (r)(x,t)r)(x f ,t')) v oc exp 



r 2 



cxp 



characteristic time and length scale r c and A c respectively. 



r') 2 



with 



(ii) The probability density of X(t) based on the dynamics similar to Eq. 
(80) has been analyzed in Ref. [108]. An application to derive an exact result 
for Dq for front propagation in a bistable system can be found therein when 
r)(x,t) satisfies Eq. (65). 

(iii) A treatment of how to calculate the correlation function (cf>(x, t)(f)(x', t)) v 
exists in Ref. [94]. The theory was then applied to front propagation in a 
bistable system and to front propagation into an unstable state (of a Fisher 
equation type model). 



3.2 Application of Novikov's Theorem to Reaction- Diffusion Systems: Effects 
of Multiplicative Noise on Specific Models 



Equipped with the necessary groundwork in Sec. 3.1.1, we now analyze the 
effect of multiplicative noise on specific models of reaction-diffusion systems. 



3.2.1 An Example of (Multiplicative) Noise- Induced Front Transitions 
Consider the reaction-diffusion equation 



f = 4? +«i -*)[«<"+*]. (88) 
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with three homogeneous stationary states, = 0, — and 1. To provide the 
state = 1 with global stability, we confine within the range [—1/2, 1] - 
with this restriction on the value of a^°\ one arrives back at Eq. (54). 

The interesting point to note in this reaction-diffusion equation is that it 
admits front solutions whose characteristics are different for different values 
of a(°\ as we illustrate below in (a)-(c). 

(a) —1/2 < a(°) < 0: In this range, = state is metastable. The linear 
marginal stability criterion does not hold in this range, and a front solution 
propagating with speed = \^D[2a^ + 1]/a/2 emerges uniquely. Let us 
denote this phase of the front by M. 

(b) < a(°> < 1/2: In this range, = is linearly unstable, but for sufficiently 
localized initial conditions, the front speed is still given by v'. For less localized 
initial conditions, front speeds higher than are also accessible. We denote 
this phase of the front by NL. 

(c) 1/2 < < 1: Finally, linear marginal stability holds in this range. For 
sufficiently localized initial conditions, a pulled front solution propagating with 
speed v* = 2y/ DaW emerges. For less localized initial conditions, front speeds 
higher than v* is also accessible. This phase of the front is denoted by L. 

How multiplicative noise can affect the characteristics of a front in a reaction- 
diffusion system is beautifully illustrated when a*- -* in Eq. (88) is replaced by 
a fluctuating quantity a(x,t) = + e l ^ 2 rj{x,t) [3,4], yielding 



& f t = + 0(1-0) [a{x,t) + </>]. (89) 

Here, r](x, t) is a Gaussian random function — in the limit of correlation length 
and time of the noise going to zero, rj(x,t) satisfies Eq. (65). 

The application of Novikov's theorem to this model yields the equivalent of 
Eq. (76), expressed as [3,4] 



D 8 ^ + v £ ^ + 0(1 - 0) [a<°> + e + (1 - 2e)0] = . (90) 



In the presence of noise, this modification of the front equation from Eq. 
reorganizes the characteristics of the front — effectively, it is this reorgani- 
zation that ends up being termed as the noise-induced front transition. The 
idea behind this transition is simple: just like the characteristics of the uni- 
formly translating front solution of Eq. (88) depend on the value of a^°\ the 
characteristics of the front solution (90) for a given value of depends on e. 
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Fig. 28. Noise-induced front transition in Eq. (90) for D = 1 [3, 4]. Solid 
lines: theoretical values for the front speed v £ = 1\J D(a^ + e) (L phase) and 
v £ = (2a(°) + l)y/D/[2(l - e)] (NL and M phase). Filled circles with error bars: 
simulation data. Symbols have been changed from the original to keep consistency 
(also the figure has been modified for greater clarity), and a should be read as a^°\ 
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Fig. 29. Wandering behaviour of the Goldstone mode for D = 1, = 0.1, 
Ax = 0.5, At = 0.1, and e = 0.125 [3]. The graph has been modified from its 
original for clarity. The notation Ax appearing in the y-axis should not be confused 
with the discretization unit Ax. 

This implies that by adjusting the value of e, it is possible to convert a front 
solution of Eq. (88) that belongs to one of the M, NL or L phases to a front 
solution (90) belonging to another one of these phases. For sufficiently steep 
initial conditions, in the L phase, the front (90) propagates with the speed 

2JD(a^ + e), while both in the NL and in the M phase, the front propagates 
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with the speed (2a^ + 1)^ D/[2{\ — e)]. In fact, the boundary between the 
NL and the L phase for the front solution (90) lies at e = 1/4 - of® /2. 67 

The v £ vs. e plots for Eq. (90) [3,4] are shown in Fig. 28. Simulation results 
for the wandering behaviour of the Goldstone mode for D = 1, =0.1, 
Ax = 0.5, At = 0.1, and e = 0.125 [3] appear in Fig. 29. 



3.2.2 An Example of (Multiplicative) Noise-Induced Fronts 

In Sec. 3.2.1, we considered a system that already admits a front solution in the 
absence of noise, and witnessed how the characteristics of the front properties 
can change when multiplicative noise is introduced in it. We now consider a 
system that does not admit a propagating front solution in the absence of 
noise, but when the amplitude of the external noise exceeds a threshold value, 
fronts start to propagate (see Fig. 30). In the front propagation literature, such 
a front is known as a noise-induced front. An example system that exhibits 
such noise-induced fronts is [105] 



f = ^-#*OM) + 2 ]- (91) 

Clearly, when a(x, t) is replaced by > in Eq. (91), the equation does not 
admit a propagating front solution — any perturbation around the state = 
simply decays in time. Elsewhere in the literature, this phenomenon is known 
as propagation failure. Similarly, the emergence of propagating front solutions 
in Eq. (91) when a(x,t) becomes a sufficiently strongly fluctuating quantity is 
known as the breakdown of propagation failure (see e.g. Ref. [57]). In fact, the 
situation is comparable to the so-called stochastic resonance paradigm [2]. 

With a(x,t) = a(°) + e 1//2 ?7(a;, t), where r)(x,t) is a Gaussian random function 
whose correlation properties reduce to Eq. (65) in the limit of infinitesimal 
correlation time and length, the equivalent of Eq. (76) obtained with the 
application of Novikov's theorem is expressed as [105] 



It is clear from Eq. (92) that when e > a^°\ fronts propagate from the stable 



67 For example, a^ = 0.1 and e = 0.2 in Eq. (90) lies exactly on the boundary 
between the NL and the L phase. Later in Sec. 3.3, we will see that this boundary 
between the NL and the L phase marks a change in the wandering properties of the 
Goldstone mode as well. 
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state st = yD[e — a(°)] to the linearly unstable state = with the (pulled) 
speed v £ = 2yjD[e-aW] (see Fig. 31). 



0.5 




Fig. 30. Snapshots of front configurations for noise-induced front propagation in 
Eq. (91) [with D = 1] when a(x,t) fluctuates around = 0.1 with an amplitude 
e = 0.15 [105]. The snapshots have been taken at t = 50, 100, 240 and 450. 




Fig. 31. Noise-induced front propagation in the reaction-diffusion system (91) [105]. 
Left figure: v £ vs. e for D = 1 and two different values of a(°b namely = 0.1 
(triangles), and a( ) = 0.3 (squares). Right figure: theoretically predicted value 
of v £ /(p st vs. simulation for D = 1. Filled symbols: simulations of the stochastic 
differential equation (91) for space and time discretization units Ax = 0.5 and 
At = 0.01 respectively; open symbols: Ax = 0.1 and At = 0.001. Circles represent 
front solution (92), and lines denote the theoretical prediction v £ = 2\J D\e — aW\. 



3.3 Fluctuating Pulled Fronts with Multiplicative Noise in Reaction-Diffusion 
Systems: Diffusive vs. Subdiffusive Wandering of the Goldstone Mode 

In footnote 67, we had commented that for the reaction-diffusion model (88) 
at a(°) = 0.1 and e = 0.2, the wandering of the Goldstone mode is subdiffusive, 
and also that this phenomenon is caused by the fact that the — 0.1, £ = 
0.2 sits precisely on the boundary between the L and the NL phase. Truly 
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speaking, the above argument is misleading — contrary to what may appear 
from footnote 67, this subdiffusive behaviour has nothing to do with noise- 
induced front transitions. Instead, the root of this subdiffusive behaviour turns 
out to go deeper — as soon as a fluctuating pulled front 68 with multiplicative 
noise belongs to the L phase, the understanding of the wandering properties of 
the corresponding Goldstone mode calls for an entirely different approach [98]. 

The reason for the requirement that a different approach is needed to study the 
wandering properties of the fluctuating pulled fronts is not difficult to trace; 
after all, for these fronts it is pretty clear that one runs into trouble with Eq. 
(81). The point is that for them the (linearized) semi-infinite leading edge of 
(O) (O in Eq. (76) behaves as £ exp[—v e £/(2D)]. This implies that for £ — > oo, 

# (0) (0 

at the leading order, the corresponding Goldstone modes $g,r = — -77 — and 

dq 

$ GL = e veZ/D d< t>^(£) behave as ~ f exp[-u e f/(2.D)] and ~ £ exp[u e f/(2£>)] 
dq 

respectively, which then send the magnitude of both the numerator and the 
denominator in Eq. (81) to oo. 

In the next paragraph, following Ref. [98], we will argue that the wandering 
width W(t) = \J (X 2 (t)) of the Goldstone mode for fluctuating pulled fronts 
with multiplicative noise scales as t 1 / 4 at long times. A more rigorous deriva- 
tion of this result by means of converting Eq. (88) for a*- -* = 1 to a KPZ 
equation for a one-dimensional interface in terms of the variable h(£,t), de- 
fined by the Cole-Hopf transformation (f)(x,t) = exp[h(£,t) — v e £/(2D)], can 
also be found in Ref. [98] . We will summarize this rigorous method later, but 
will not elaborate on it. 

The argument for the t 1 / 4 scaling behaviour of W(t) at long times [98] is 
the following: notice that before even one arrives at Eq. (81), a more severe 
problem already occurs at the level of Eq. (75) for fluctuating pulled fronts 
- namely that when is pulled, then there is no inherent time scale for 
the convergence of the front speed to its asymptotic value [c.f. Eq. (4)]. We 
had earlier seen that this behaviour originate from the gapless spectrum of the 
stability operator (77) for pulled fronts. It implies that separating the effect of 
the fluctuations into a "slow" wandering of the Goldstone mode from the "fast" 
modes of fluctuations in the front shape, upon which the whole underlying 
philosophy of Eqs. (75) through (81) is based, does not make sense any longer. 
Nevertheless, the t 1 / 4 scaling behaviour of W(t) emerges beautifully when the 
long-time behaviour of (/>^(^t) is successfully married to Eq. (81), and that 



68 What we mean by fluctuating pulled front here is a fluctuating front with multi- 
plicative noise belonging to the L phase. Here our usage of the phrase fluctuating 
pulled front is chosen to distinguish these fronts from the fluctuating "pulled" fronts 
that we encountered in Sec. 2. 
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is carried out in the following manner. 

With the result that the leading edge of the pulled front (f>^°\^) relaxes asymp- 
totically as [38] 



(O) (^) ~ ^ e -^/(^)-e/^Dt) /t 3/2 for u ^ >1 (g3) 

one defines a time-dependent upper cutoff £ c ~ y/ADt for the £ integral in the 
denominator of Eq. (81). When this upper cutoff is inserted in Eq. (81), one 
obtains a time-dependent diffusion coefficient for the Goldstone mode [98] 



3e 
l ypM 



(94) 



for t > 1. The implication of Eq. (94) then is that for t > 1, W(t) could 
possibly scale as t 1//4 . 




Fig. 32. Log-log plot of W(i) vs. t showing subdiffusive wandering of the Goldstone 
mode for fluctuating pulled fronts with multiplicative noise (89) [98]. The data 
correspond to = 1. Solid line: simulation data, dashed line: a line with slope 
1/4. 

To obtain the t 1//4 scaling of W(t) in terms of the rigorous Cole-Hopf transfor- 
mation of Eq. (89) to a one-dimensional KPZ equation [98], one invokes the 
scaling relation W(t) = t a F{t/L^). When the argument \i of F(fi) is 1, 
then F(n) behaves like a constant, and when \x ^> 1, then F(fi) ~ fi~ a . The 
proper way to implement the above scaling relation is to use a time- dependent 
length scale L(t) ~ £ c ~ y/ADt. For a one-dimensional interface, the KPZ 
exponents are a = 1/3 and (5 = 3/2, which indicates that for t — > oo, the 
argument of F also approaches oo. The long time behaviour of W(t) is then 
clearly given by t 1 / 3 [t/( v / t) 3/2 ]- 1/3 = t 1/4 . 

The long time t 1 / 4 scaling behaviour of W(t) for the model (88) with = 1 
has been confirmed by numerical simulations [98] (see Fig. 32). We note here 
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that although the wandering properties of the Goldstone mode for the model 
(91) was not studied in Ref. [105], from the argument based on the upper 
cutoff £ c for the integrations in Eq. (81), in this model too we expect the 
same t 1 / 4 scaling behaviour for W(t) at long times. A rigorous derivation of it 
based on the Cole-Hopf transformation in the spirit of Ref. [98] is left for the 
enthusiastic readers. 



3.4 Fluctuating Fronts with Multiplicative Noise in Reaction-Diffusion Sys- 
tems and Kinetic Roughening 

While the Cole-Hopf transformation allows one to rigorously derive the t 1//4 
scaling of the wandering properties of the Goldstone modes of fluctuating 
pulled fronts with multiplicative noise in one spatial dimension, questions can 
certainly be raised regarding its applicability to higher dimensions. Surpris- 
ingly enough, it turns out that the Cole-Hopf transformation is equally suc- 
cessful for reaction-diffusion systems in two spatial dimensions, only if for 
£ — > 00, the multiplicative noise couples linearly to the front field. Such a 
reaction-diffusion equation 

^ = DV 2 + [1 + Vhj{x, t)W - 2 ] , (95) 

in two spatial dimensions has been considered in Ref. [115]. The variable h then 
corresponds to a two-dimensional interface (see Fig. 33). The front wandering 
properties in Eq. (95) has then been shown to obey the KPZ scaling for a 
two-dimensional interface 69 — see Fig. 34. 70 

69 In the scaling W(i) ~ t Q F(t/L^) for two-dimensional interfaces, existing accepted 
values are a = 0.245 ± 0.003 and a/3 = 0.393 ± 0.003 [77]. One has to be somewhat 
careful in how to interpret "long time" wandering exponent in Fig. 34. For fluc- 
tuating pulled front with multiplicative noise in one spatial dimension, we earlier 
saw that L itself scaled as y/i. Since for a one-dimensional interface, j3 = 3/2, at 
large times, the argument of F then behaved as t 1//4 — > 00, and that in turn im- 
plied that as t — > 00, F(t/LP) behaved as [t/L@]~ a . However, it turns out that for 
fluctuating pulled fronts with multiplicative noise in two spatial dimensions (95), 
even for long times, the argument of F remains small, resulting in the behaviour of 
F(t/L,P) ~ const, and correspondingly W(t) behaves simply as t a [115]. 

70 At this juncture, one cannot help but wonder whether the fluctuations in the 
discrete particle and lattice model of Ref. [116] have really anything to do with 
fluctuating pulled fronts with multiplicative noise that behaves as e 1 / 2 ^ as <p — ► - 
after all both observe the same scaling for W(t) in one and two spatial dimensions. 
While we are leaving the major question of how well the effect of the internal 
fluctuations due to the discreteness effects of the particle and and the lattice are 
represented by fluctuations in field-theory for Sec. 4, at this short interlude, we can 
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Fig. 33. Fluctuating pulled front propagation with multiplicative noise in two spa- 
tial dimensions for the model (95). The front propagates in the x direction, and 
£ = x — v £ t. Left figure: an instantaneous snapshot of the front configuration. Right 
figure: the corresponding configuration in terms of the two-dimensional interface 
"height" h(£, t) obtained by the Cole-Hopf transformation. The solid curve denotes, 
in both figures, the contour for position (£, y) where the value of <p is 1/2. 
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Fig. 34. Log-log plot of W(£) vs. t showing the wandering of the Goldstone mode for 
the fluctuating pulled front with multiplicative noise in Eq. (95) [115]. The with and 
without nonlinearity respectively indicate whether nonlinear terms oc exp[/i(£,i)] 
generated through the Cole-Hopf transformation of Eq. (95) were retained in the 
simulation or not. 

The success of the Cole-Hopf transformation to reduce reaction-diffusion equa- 
tions with multiplicative noise in d spatial dimensions to a KPZ equation for 
a d dimensional interface is certainly extremely intriguing. Although in any 
spatial dimensions, this transformation works only if for £ — > oo, the fluctua- 
tions in the external parameter linearly couples to the front field, one wonders 
if there exist further general connections between the evolution equations of 
fronts and those describing interfacial growth phenomena. From this perspec- 
tive, a beautiful analysis for fronts (even with vector front fields) in reaction- 
diffusion equations in arbitrary dimensions has been presented in Ref. [99]. 

say that the results of Ref. [80] points the finger to a negative answer. We will further 
see in Sec. 4.1 that the internal fluctuations are usually modeled by a fluctuating 
term in the stochastic differential equation that behaves as y/4> for <f> — > 0, which 
too, rules out an affirmative possibility. 



— With Nonlinearity 




84 



The analysis is not rich in mathematical rigour, and we will not describe its 
mathematical procedure here either. The idea therein is the same as what 
has been already presented in Sec. 3.1.2; namely that one needs to perform a 
separation of length and time scales. As the short time and length scales are 
averaged out by means of a coarse graining method, leaving behind the long 
length and time scales in a "renormalized" evolution equation of the front, 
one does recover the dynamical equations for growing interfaces, standardly 
seen in interfacial growth literature. 

Applied to specific systems of Sec. 3.2.1, this formalism yields the same results 
for the expression of the shift in front speed and the diffusion coefficient of 
the Goldstone mode in one dimension, as expected [99] . It also reproduces the 
front fluctuation characteristics in Fig. 34 for the reaction-diffusion system 
(95) in two dimensions. 



3.5 Epilogue II 

Generally speaking, (Langevin-type) field-theoretical methods have been used 
ubiquitously by chemists to describe the effect of external noise on chemical 
reaction-diffusion systems for decades. In many of these chemical systems, one 
does observe front propagation. From that point of view, it has been difficult 
for me to decide, for Sec. 3, what to significantly include, what to mention in 
appropriate context and what to leave out. In fact, in Sec. 3, I have adopted 
the philosophy of discussing only those papers in detail that have contributed 
significantly to the field-theoretical framework on the effect of external noise 
on propagating fronts. In this process, I have strived to provide a unified 
perspective of these field-theoretical methods. 

While all the works that I am aware of on these field-theoretical methods have 
exclusively been on reaction-diffusion systems, quite expectedly, there exist a 
few scattered results that concerned weak noise, or simply additive fluctua- 
tions. These results, in addition to some of the real laboratory experiments 
that in the first place motivated the field-theoretical studies of multiplicative 
noise on propagating fronts can be traced by following the cited references in 
the works that we have discussed in Sec. 3. 

There is one further point to that I would like to bring the reader's attention 
to. That is the fact that effectively majority of the works on field-theoretical 
treatments of propagating fronts with multiplicative noise considered only 
spatiotemporally white noise — after all, spatiotemporal white noise in the 
ensuing analytics is a lot easier to handle than coloured or structured noise. 
Although we started with a fancier form of the noise (62) in Sec. 3, we reduced 
it to spatiotemporal white noise (65) very soon in the limit of both A c and 
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This however does not imply, in any way, that structured noise forms have 
not been considered in the literature to study the effect of multiplicative noise 
on propagating fronts. Few examples for weak coloured noise forms have been 
mentioned in points (i)-(iii) at the end of Sec. 3.1.2. A recent work has ap- 
peared with spatiotemporally structured (and not necessarily weak) noise [106] 
- wherein an approach along the lines of Novikov's theorem 71 [106] has been 
adopted. The central theme of this line of thought is, once again, to rewrite 
Eq. (61) in the form of Eq. (72), such that the stochastic term does have a 
zero mean. In the absence of the applicability of Novikov's theorem, in this 
case, one carries out a cumulant expansion of the noise term to evaluate the 
l.h.s. of Eq. (70) [106] (the precise mathematical procedure in this cumulant 
expansion is non-trivial; but naturally, in the limit of both r c and A c — > one 
recovers the white noise results). The limitation of this technique, however, is 
that it is a perturbative expansion in (small) correlation length and time, and 
therefore, one should not expect more than some limited success. 

Before we end Sec. 3, let us return to the paragraph below Eq. (84). We 
had commented therein that Dq is defined only when R(4>, £, t) is replaced 
by R(<j>(°\£,t) in Eq. (79). While Eq. (80) is still well-defined without this 
replacement, Eq. (81) is not. First of all, note that without having _R(0,£,t) 
replaced by R(^°\^,t), R is no longer white both in space and time [c.f. Eq. 
(74)], and that certainly gets one into trouble. Secondly (and more impor- 
tantly), taking the average over the noise alone to obtain an expression of D G 
is not enough. 

To make the second point more articulate, let us take a note of the fact 
that the average over the noise, such as in Eq. (74), implies that for t' > t, 
one averages over all possible noise realizations after t. With this average 
alone, the expression for Dq [without having R((j), £, t) replaced by R((f>(°\£, t)] 
depends on the precise front configuration at time t. Thereafter, one still has 
to also average over the ensemble of realizations of the front configurations at 
time t [just like in Eq. (28)] to properly define D G . In that case, to obtain a 
corresponding theoretical expression of Dq for general noise terms is not an 
easy task (we will see more of it in Sec. 4.2.1. 

Nevertheless, whether one can replace R[(/>(£, t), £, t] by R[(f>(°\£, t), £, t] under 
general conditions is a natural question. A priori, it seems to be a trivial 
issue, because for weak noise, one is always accustomed to think that such a 
replacement is always possible. The point however is that since we are dealing 
with fluctuating fronts that propagate into unstable states, i.e., f g ((f>) ~ (f> for 
0^0, one must always bear in mind the following note of caution: in the limit 

71 Recall that Novikov's theorem is applicable only when the noise is Gaussian. This 
requirement renders it unusable for coloured noise with non-Gaussian statistics. 
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of weak noise, so long as R[cj>(£,t),£,t\/(f> for — > (all the cases that 
we discussed in Sec. 3 follow this property), R[(f>(£,t),£,t] can be replaced by 
R[<f)( \£,t),£,t]. Indeed, in Sec. 4, we will see this condition violated, namely 
that for the sFKPP equation (96) R[<j>(£, t), £, t]/<f> for ~ l/N and 
N — > oo. In that case, if one does replace -R[0(£, t), £, t] can be replaced by 
i), £, t], one gets a front diffusion coefficient scaling as l/N, which is 
certainly wrong! 



4 Field-theory of Fluctuating Fronts: Internal Fluctuations 

In first three paragraphs of Sec. 3, we have already got a taste of the compli- 
cations associated with internal fluctuations. It therefore should not come as 
a surprise that (unlike what we experienced in the field-theory of fluctuating 
fronts for external fluctuations) describing the effect of internal fluctuations on 
propagating fronts is a real challenge for the field-theorists. Before we plunge 
into the details in Sec. 4.1 onwards, let us demonstrate the major intricacies. 

To start with, let us remind ourselves that all the works on the field-theoretical 
treatments of the effects of internal fluctuations on propagating fronts that I 
am aware of have been exclusively in the reaction-diffusion systems. However, 
as far as Sec. 4 is concerned, the readers must not overconclude the reference 
to reaction-diffusion systems in the context of field-theory for the effects of 
internal fluctuations on fronts in general. To make this point clearer, I want 
to draw the reader's attention to the fact that there exists a vast literature 
on field-theoretical treatments of the effect of internal fluctuations on fronts 
in various reaction-diffusion systems, where two opposing currents of reactant 
species of particles meet from two sides and reactions between them take place 
in a reasonably localized volume of space known as the reaction zone. In the 
reaction zone, one can model the fluctuations in particle densities by means 
of a field-theory, and one obtains various kinds of spatiotemporal behaviour 
of particle densities or correlations between reacting particles etc (see for ex- 
ample, Refs. [64,65]). In these systems reaction fronts do develop, but they do 
not propagate (asymptotically) with a constant speed. The dynamics of these 
fronts are extremely interesting on their own merit, and their understanding 
has contributed significantly to renormalization group techniques, but as far 
as this review article is concerned, they are beyond our purview. 

On the other hand, as far as internal fluctuations in propagating fronts are con- 
cerned, Langevin type field-theories of react ion- diffusion fronts were of interest 
first to chemists [67,68] and mathematicians [83]. To be more particular, the 
interest was in a Langevin description of stochastic Fisher equation. It would 
still be another few years until physicists found field-theories for propagating 
fronts to be of considerable interest [35,36,95]. 
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So why are field-theories of internal fluctuations in propagating fronts so dif- 
ficult? The answer to this question, in fact, lies at two levels. The difficulty at 
the first level stems from the complication associated with the description of 
fluctuations in a real experimental situation. In the second paragraph of Sec. 
3, we discussed this briefly using the example of a chemical reaction — that 
in general, in a chemical reaction, not only are there local fluctuations in the 
density of reactants, but also complicated correlations between fluctuation in 
particle densities and in the mesoscopic parameters such as reaction rate or 
diffusion coefficient of particles exist. Taking all these into account in a reason- 
able and satisfactory degree of detail is indeed very difficult, and other than 
Ref. [69] in the chemistry literature, I am not aware of any treatment where 
the effect of local fluctuations in the mesoscopic react ion- diffusion parameters 
on the propagating front has been considered. 

To get a good theoretical grip on reaction-diffusion processes, just like we 
have seen all along in Sec. 2, in all the field-theoretical approaches based on 
Langevin type stochastic differential equations, it is customary to consider 
fluctuations only in the density of particles while the usual mesoscopic param- 
eters such as diffusion coefficients or reaction rates are held constant. Even 
then, these field theories are not free from intricacies. A glance at the results 
(a)-(c) below in the chemistry literature [67,68,70] will make this clear: 

(a) For the react ion- diffusion process X+Y^ 2X described by the Fisher equa- 
tion (1) with n — 2, where the Langevin forces have been obtained from the 
master equation, a front speed higher than v* has been observed in Ref. [67]. 
It appears from the results of the numerical simulation that in the range 
N ^ 10 2 -10 3 - 5 , vn behaves as v n — v * ~ 

(b) In a Ginzburg-Landau or Schlogl model (55), a similar Langevin formalism 
deduced from the master equation also yields a front speed higher than [68]. 
Numerical simulations in this case yields vn — ~ JV -1 - 38 . 

(c) Direct Monte Carlo simulation techniques applied to the reaction-diffusion 
process X+Y— > 2X described by the Fisher equation (1) with n = 2, however 
yields a decrease in the front speed from v*, and as a function of N in the 
range N ~ 10-10 3 5 , v*-v N ~ iV" 1/3 [70]. 

Note that given our prior knowledge of how vn behaves as a function of N 
from Sec. 2, it is difficult to reconcile (a) and (b) with anything else [even 
with (c), which in itself is not the asymptotic behaviour of vn for N — > oo]; 
and this is precisely where the second level of difficulty comes in. The point is 
that the Langevin type field-theories derived from the master equation effec- 
tively are expansions over the "small parameter" l/((iV/ c )), where ((Nk)) is the 
conditionally averaged number of particles on lattice site k (see for example 
Chap. VII of Ref. [52], or Ref. [16]). Clearly, the l/((N k )) expansion is only 
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possible if ((Nk)) are large quantities — so while the expansion is valid in the 
bulk phase of the front, it break down at the tip of the front, where there are 
a very few particles per lattice site. At least for a fluctuating "pulled" front 
such as the stochastic Fisher equation, we already know that its properties 
are tremendously sensitive to the fluctuation dynamics at the tip of the front. 
For this reason alone, such Langevin type field theories has to be applied with 
extreme caution, as we illustrate in the next two paragraphs below. 

For fluctuating "pulled" fronts, the reason why the results from the chemists' 
side [67,68,70] contradict the l/ln 2 iV convergence of the front speed v N to 
v* from below is not difficult to trace — it lies in the fact that the form 
of the stochastic terms that they consider amounts to an additive noise to 
the deterministic equation. If we compare the resulting evolution equation of 
the front with the microscopic reaction-diffusion process X^== 2X for example, 
the contrast between the physical implications of the stochastic differential 
equation and what actually happens in the microscopic process becomes really 
clear. In this microscopic process, lattice sites on which there is no particle to 
start with cannot grow a particle on its own, meaning that the first particle 
on any lattice site has to diffuse from one of its nearest neighbours. When this 
idea is applied to the lattice sites that lie ahead of the front, one immediately 
recovers the idea of the i.f.o.l.s. (see footnote 36). The lattice sites on the 
right of the i.f.o.l.s. cannot generate particles on their own, and that was the 
motivation of Brunet and Derrida behind using a growth cutoff for < 1/N. 
The growth cutoff then predicted that the leading edge of the front would be 
of finite extent, and it also gave rise to the l/ln 2 iV convergence of the front 
speed v N to v* from below for fluctuating "pulled" fronts. 

An additive noise term totally changes this scenario. In that case, any place 
ahead of the actual front region where value is identically zero at a given 
time can generate a non-zero <f) within an infinitesimal time interval simply 
due to fluctuations. These fluctuations then have the capability of pulling the 
rest of the front faster than even the leading edge of the deterministic pulled 
fronts! It is therefore no wonder that these models [67, 68, 70] do yield front 
speeds that are higher than v*. 



4-1 Fluctuating "Pulled" Fronts in Stochastic Fisher-Kolmogorov-Petrovsky- 
Piscunov (sFKPP) Equation 

From the discussions in the last few paragraphs above, it is clear that as far 
as a stochastic differential (Fisher-type) equation is concerned, additive noise 
to model fluctuating "pulled" fronts is out of the question. The noise value 
has to be zero where the value is zero and this requirement then points the 
finger to multiplicative noise. 
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While this argument is generally valid, question remains regarding what form 
of multiplicative noise one should choose. Clearly, one should certainly avoid 
a noise term that behaves ~ e 1 ^ 2 ^ for 0—^0 — for example the case of 
fluctuating pulled fronts with multiplicative noise as we encountered in Sec. 3. 
For such noise terms, no matter howsoever small is, the noise amplitude is 
always an order of e <C 1 weaker than 0, and one does not expect such noise 
terms to yield a leading edge of a finite extent for the front (although I have 
not seen any work proving or disproving it). 
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Fig. 35. A snapshot of the fluctuating "pulled" front with D = 1 in sFKPP equation 
(solid curve) demonstrating compact support property [35,36]. The corresponding 
plot for the deterministic pulled front in Fisher equation (1) [n = 2] is shown with 
dashed lines for comparison. 

It turns out, however, that for the front in Fisher-type (reaction-diffusion) 
equation, the most frequently used noise term that does yield a leading edge 
of a finite extent in fact behaves ~ e rl//2 \/0 for — > (of this, the e < 1 and 
e ^> 1 limits are respectively referred to as the weak noise and the strong noise 
limits). The corresponding stochastic differential equation is then known as 
the stochastic Fisher-Kolmogorov-Petrovsky-Piscunov, or in short, the sFKPP 
equation [35, 36, 83, 95] : 




f = + - 2 + V (x, t) . (96) 

Here r](x,t) satisfies the properties that {r](x,t)) v = and {r](x, t)r](x', t')) v = 
25(x — x')5(t — t'). The stochastic term is interpreted in the Ito sense. The 
factor of 2 in the noise correlation is really unnecessary, but it is introduced 
to maintain consistency with the corresponding form (65) in Sec. 3. 

For the microscopic reaction-diffusion process 2X however, the sFKPP 
equation (96) has been derived in terms of a path integral formalism [95], 
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which avoids the subtleties associated with very few particles per lattice site 
at the leading edge of the front. In this derivation, one first rewrites the full 
master equation in a path integral form, and then proceeds with the idea 
that the fluctuations in the front around the mean front shape are always 
small and retains the effects of the fluctuations only at the leading order. 
The compact support property for sFKPP equation — not only that there 
exists an x r (t) for which <p[x > x r (t),t] = 0, but also there exists an xi(t) 
for which (f>[x < xi(t),t] = 1, and that [x r (t) — xi(t)] < oo Vi — has been 
rigorously derived too [83] (see Fig. 35). Although none of these derivations 
will be discussed here any further, note that the compact support property 
of the fronts in sFKPP equation is not a surprise. Both around ~ e and 
~ 1 — e 1 / 2 , the strengths of the stochastic term are respectively of the same 
order as and 1 — 0, and that helps to abruptly cut down the growth of 
ahead of and the decay of 1 — behind the front region. 
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Fig. 36. Speed for the fluctuating "pulled" front in the sFKPP equation on a discrete 
lattice with D = 1 and e = 1/2N [95]. Circles: simulation data (with error bars), 
solid line: the theoretical expression (11). The symbols have been changed from the 
original to maintain notational consistency. 

In the weak noise limit, e is replaced by 1/2N (the factor of 2 is put in to 
counter the 2 appearing in the expression of (r)(x, t)r)(x', t')) v ) for N — > oo to 
make contact with the react ion- diffusion system 2X. Even in this limit, 
due to the non-local nature of sFKPP equation, it is not possible to obtain a 
theoretical expression for the speed of the corresponding fluctuating "pulled" 
front (although some simple local forms of it can be solved exactly at all noise 
strengths [95]). Nevertheless, one can put it in the computer [95]. For N — > oo, 
the speed of the fluctuating "pulled" front is then seen to confirm the 1/ In N 
convergence (11) of the front speed vn to v* from below (see Fig. 36). 
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4-1.1 Fronts and Duality between the Weak and the Strong Noise Limits in 
the sFKPP Equation 

The sFKPP equation, in fact, is more versatile than what appears from above. 
What is required is that one puts it in appropriate perspective. 



Fig. 37. Log-log plot of vjy/v* in the sFKPP equation on a discrete lattice as a 
function of noise strength e 1 ! 2 [35,36]. Open circles: simulation data for the sFKPP 
equation, solid lines: the theoretical expressions (11) at weak noise £«1 and (97) 
at strong noise £>1. 

Consider once again the reaction-diffusion process 2X. We had earlier 
seen, in Sees. 2.3 and 2.6 that N, the number of particles per lattice site at 
the stable phase of the front, dictates the strength of the fluctuations. The 
smaller N is, the (relatively) stronger the fluctuations become (c.f. Sec. 2.6). 
On the other hand, since N is essentially the ratio between the rates of the 
forward and the backward reactions in the 2X process [16,35,36,87], in 
Eq. (96), the noise strength e ~ 1/N is decided by this ratio. 

The formal connection between the noise strength and the ratio between the 
forward vs. backward reactions in the reaction-diffusion process X±=? 2X has 
recently been rigorously proved through the notion of duality [35,36]. This 
connection then enables one to predict the front speed at strong noise e ^> 1 
in the following way: the strong noise limit reduces the X^ 2X system to 
the model of Sec. 2.6.5 in the e — > and W — > oo limit, 72 while keeping 
D ~ 0(1). In this limit N = e/(e + W) ~ e/W, and then very easily, Eq. (51) 
predicts a front speed 

72 Although the 2X system formally allows more than one particles per lattice 
site, effectively no lattice site ever gets occupied by more than one particles when 
the backward reaction rate is much higher than the forward reaction rate. The 
condition for the model of Sec. 2.6.5 that each lattice site can be occupied by at 
most one particle at any given time is then automatically satisfied. 
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vn = DN ~ — 

e 



(97) 



The rescaled front speed vn/v* = vn /(2yDe) is plotted as a function of 
the noise strength i 1 ' 2 in Fig. 37 [35,36]. The open circles denote simulation 
results for the sFKPP equation and the solid lines denote the corresponding 
theoretical results for weak [Eq. (11)] and strong [Eq. (97)] noise. 



4-2 The Dj and D G Dilemma 



From what we learnt in Sec. 4.1, it seems that the sFKPP equation, interpreted 
in the It 6 sense, does provide a lot a of input to the subject of fluctuating 
"pulled" fronts — it correctly predicts the front speed behaviour as a function 
of the noise strength for the limiting cases of both the weak (very large values 
of N) and strong noise (N <C 1). 

There is, however, a missing piece of the puzzle in the success story of the 
stochastic differential equations for internal fluctuations, namely that the front 
diffusion coefficient in one spatial dimension has been explored in neither for 
the strong and nor for the weak noise. As the issue of front diffusion now 
becomes our focus, at the very first step, it becomes a necessity that for 
general reaction-diffusion systems, we reconcile Df of Sec. 2.5 and D G of Sec. 
3, both introduced for the measure of front diffusion. 

4-2.1 Df vs. D G : What Exactly is the Difference? 

We start with the ltd stochastic differential equation for e <C 1 



^ = + f(<P) + iV*R[<l>{x,tlx,t\ (98) 

[which is of the same form as Eq. (73)]. Here R[(j)(x, t),x, t] = g[4>(x, t)]rj(x : t), 
and in the Ito interpretation, 



{R[(j)(x,t),x,t])r, = g[(p(x,t)} (rj(x,t)) v = and 

(R[<Kx, t),x, t\R[(t>{x\ t'),x', t']) v = g[<j>{x, t)\g[<j>(j> , t')\ ( V (x, t) V (x', t')) v 

= 2g 2 [<i){x,t)}5(x-x')5{t-t'). (99) 

In Eq. (98), Dq is defined from X(t), the random displacement of the Gold- 
stone mode of the front [see Eq. (80)]. On the other hand, if one measures 
the position of the front by its centre of mass location using the continuum 
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version of Eq. (24), then from Eq. (78) one obtains, for the fluctuation in the 
speed of the centre of mass of the front, 



S(t) = X(t) + J d£6fa,t), (100) 

— oo 

where 6<f>(t) is measured on the instantaneous comoving frame of the Goldstone 
mode. 

Next, having simplified S(t) with the aid of Eqs. (77) and (79), and with the 
inherent assumption that 6(f>(£, t) — > for £ — > ±oo as 



S(t) = 1 di 5 -^ 5<j>(U) +e 1/2 ]dtR[<f>(U),U], (101) 

J 0(f) 6(0) J 

— oo — oc 

we immediately notice something interesting - - namely that through Eq. 
(101), we have essentially recovered Eq. (27). Remember that Eq. (27) was 
obtained through a heuristic argument for the clock model! For the uniformly 
translating front solution of the deterministic react ion- diffusion equation 

oo 

d t (f) = Dd 2 x <\) + f(4>), the front speed is given by J dx f[(p^}, and naturally, 



oo 

deviation in the front speed is given by J dx- 



5<fr(x, t) just like we had 

0(0) 



when the distribution of (j>(x,t) deviates from 0' O) , at the leading order, the 

seen in Eq. (26). Secondly, in Eq. (101), we have not replaced R[(j)(£, t), £, t] by 
R[(f>(°)(£, t),£, t] as we did in Sec. 3. In fact, the reason we had to do so in Sec. 
3 is to make the stochastic term in Eq. (72) ^-correlated both in space and 
time in the Stratonovich interpretation. In the Ito interpretation (99) this re- 
placement is no longer necessary. The choice of not replacing R[(j>(£, t), £, t] by 
i?[0(°)(£, t), £, t], however, implies that to obtain the front diffusion coefficient, 
we would first have to average over the noise realizations (subscript rj) and 
then average over an ensemble of initial front realizations at time t (subscript 
f), just like in Eq. (28): 



D S = \ Mm J dt' ((S(t) S(t + 0>„>« • (102) 
o 

Thirdly, with (R[cf>(£, t), £, t]) v = 0, the average of Eq. (101) over all possible 
noise realizations at time t does not yield (S(t)) v = 0. Indeed, this is the reason 
why we had a non-zero 5v Ttm {(t) in Eq. (25). Only when a further average over 
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all possible front realizations at time t is taken, by virtue of (S(f>(£,t))t = 0, 
one gets (S(t)) = 0. This is also the same as in Eq. (27). 

Now to simplify Eq. (102). The r.h.s. of Eq. (101) contains a sum of two terms, 
and the product of S(t)S(t+t') that appears in the Green-Kubo formula (102) 
contains four terms. Of these, an average over all possible noise realizations 



after time t kills the 



5/(0) 



0(0) 



s<f>(u) J deR[<f>(?,t+t f ),e,t+t f ] but 

— oo 

the rest three survive. Among these three surviving terms, the easiest one to 

oo oo 

handle is J d£R[0(f J d^R[<f>(^,t + t'),£,t + t'}. With Eq. (99), its 

— oo — oo 

contribution to Df is seen to be 



D (i) 



dt(g 2 m,t)}) 



t ■ 



(103) 



To obtain the contribution D { f of fdZR[<f>(£,t),Z,t] Jd£ 

— oo — oo 

we formally solve S(j)(^,t + t') first by expanding it as in Eq. (82), and then by 
using Eq. (85) (without having replaced R[<j>(£, ?'),£, t"\ by i?[0 (o) (£, *"]). 
Thereafter, Eq. (99) and the integration over t' yield 



oo oo 



-oo -oo 



• (104) 

4,(0) (e) 



This is still not simple enough. To reduce Eq. (104) further, we return to Eqs. 



(16) and (17) and use the fact that [d^mAO^Tp- = r' 1 [d^mAO 

— oo — oo 

With the completeness relation (86), we then get a simple expression: 



Df = e J dt(g 2 mt)]) t 



(105) 



Finally, proceeding exactly along the lines of Eqs. (104) and (105), D f , the 



contribution of the / d£ — - — S4>(£, t) [ d£ 

J 0(f) 6(0) J 

— oo — c~ 

Df can also be simplified to 



Jf{<t>) 



4,(0) 



5<j)(€', t + t') term to 
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oo oo 

D f = E T m l (c m (t) Cm ,(t)) t J di J dC'^mAO^m'AO- (^6) 



To obtain Df = Dj + Df + Df , one now has to calculate the aver- 
age of the quantities () t in Eqs. (103), (105) and (106), i.e., over an en- 
semble of (initial) realizations at time t. When R[cj>(£,t),£,i\ is replaced by 
i?[0(°)(£, t), £, t] in the Stratonovich interpretation of Sec. 3, this is an easy 

oo 

task: while J dt; <7 2 [0^(£)] * s ver Y eas Y to evaluate directly, Eq. (88) effec- 

— oo 

tively becomes a Langevin equation (with additive noise) for a Brownian parti- 
cle in a viscous fluid. One can then simply use fluctuation-dissipation theorem 
to evaluate (c m (t)c m >(t)) t . In the Ito interpretation, however, one does not 
require to replace R[tfr(£,t),£,t] by R[(j)^(^,t),^,t] anywhere between Eqs. 
(100) and (106), and then strictly speaking, one cannot use the fluctuation- 
dissipation theorem for multiplicative noise to evaluate (c m (t)c m '(t)) t from 
Eq. (88) any longer — is is a well-known problem in Stochastic processes that 
there is no fluctuation-dissipation theorem for multiplicative noise (see, e.g., 
Ref. [52]). 

Nevertheless, what is the problem if we simple-mindedly replace R[<f>(£, t), £, t] 
by t), £, t\? It certainly makes calculations easier. As for the answer 

to this question, we note that one really has to be very careful to do this. 
The integrand in Eq. (87) contains x 3/ m ,L(0 = e" e ^ 2D Vm(0; an d due to the 
presence of the exponential weight factor e Vs ^^ 2D \ it is the fluctuations at the 
tip of the front that matter the most. This further implies that c m 's are strong 
functions of time in general — after all, in case of internal fluctuations of 
fronts, the fluctuations at the tip of the front are of the order of the front 
value itself. 73 However, when the nonlinearity /(</>) makes the front pushed, 
i.e., for fluctuating pushed fronts, all fluctuation modes have a finite lifetime 
comparable to the time scale set by l/v £ . In that case, the dependence of the 
precise value of 5<p in one snapshot of a front realization at any time t at 
the tip of the front does not matter so much on the precise noise realization 
that has been used to update the front profile before t. This then justifies 
the replacement of R[(j)(£,t),£,t] by R[(j)(°\£,t),£,t] [and further on the us- 
age of fluctuation-dissipation theorem to obtain (c m (t)c m >(t)) t ~ 0(e)}. For 
fluctuating "pulled" fronts however, such as the sFKPP equation the fluctu- 
ation modes can have very long decay times (c.f. Sec. 2.4.2). Then the the 
replacement of R[(f>(£,t),£,t] by R[<f>(°)(£,t),£,t] and the subsequent usage of 
fluctuation-dissipation theorem to obtain (c m (t)c m i(t)) t is simply incorrect! 



73 This actually depends on the functional form of g(<fi). For example, when g{4>) ~ 
y/4> for (f> — > 0, the fluctuations at the tip of the front are certainly as strong as 4>. 
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For fluctuating pushed fronts, when R[<f>(£,t),£,t] is replaced everywhere by 
R[<f>(°\£, t), £, t], one can easily show that = D^p = 2D^p [after using 
fluctuation-dissipation theorem in Eq. (84) and thereafter with the complete- 

oo 

ness relation (86)], yielding Df = 5e J d£ (? 2 [0^(£)]/2. This expression is 

— oo 

clearly different from that of D G in Eq. (81). 74 As for fluctuating "pulled" 
fronts, the scalings of Df and D G will now be dealt with in Sec. 4.2.2. 



4.2.2 D f andD G for the sFKPP Equation 

From the point of view that the sFKPP equation belongs to the family of Eq. 
(98), a question that naturally arises is whether it is possible to apply our 
collected wisdom of Sec. 4.2.1 to the sFKPP equation. 

A priori, one can expect such an attempt to fail, or at best expect it to be 
full of subtleties, for by now, we have already learnt from Sec. 3.3 that in the 
limit of zero noise strength, whenever Eq. (98) gives rise to a pulled front, D G 
[Eq. (81)] does not remain well defined any more. This simply has to do with 
the fact that the integral in the denominator of Eq. (81) diverges due to the 
semi-infinite leading edge of pulled fronts. 

We can however ask ourselves the following question: from Brunet and Der- 
rida's analysis, we know that an effective description of fluctuating "pulled" 
fronts in Fisher equation is in terms of a growth cutoff below <p ~ 1/JV. So 
how about using Brunet and Derrida's cutoff solution for in the whole 
formalism of Sec. 4.2.1, assuming that the shape fluctuation of the front in 
the sFKPP equation are indeed described by the shape fluctuation modes 
{\l/ m }? 75 At least in that case, any complication associated with the diver- 
gence of the denominator Eq. (81) is safely avoided! 

What we will show in this section is that such an simple-minded approach 
works [91]. It yields the asymptotic scalings D G ~ 1/ In 6 N and Df ~ 1/ In 3 N. 
Although at the time of writing this review article, neither D G nor Df for the 
sFKPP equation has been studied (analytically or numerically), we will later 
argue that the asymptotic scaling Df ~ 1 / In 3 N is a generic property of 



74 For fluctuating pushed fronts made of discrete particles on a lattice, e ~ 1/JV. 
This indicates that for these lattice models, one does obtain a Df that scales asymp- 
totically as 1/JV. From Eq. (81), D G too is seen to scale as 1/JV. 

75 Recall that in this case, {^ m } are localized within the front region and they 
vanish for £ — > ±oo (c.f. Sec. 2.4.2). This is precisely the compact support property 
of the front solution in the sFKPP equation entails; and moreover, it is the same 
requirement that underlies Eq. (100). This means that all the results we obtained 
in Sec. 4.2.1 can be safely applied to this case. 
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fluctuating "pulled" fronts. 

Between Dq and Df, the l/\n 6 N asymptotic scaling for Do is the easiest 
to demonstrate. For simplicity, we set D — 1. With e = 1/(2N) and g(4>) = 
y/cj) — (f) 2 for the sFKPP equation, from Eq. (81). we obtain 




When one uses Eq. (10) [e = 1/N] in Eq. (107), the exponential weight fac- 
tor in the integrand of the denominator drops out and the denominator itself 
asymptotically scales as 1/ In 6 N at the leading order. Dealing with the inte- 
grand of the numerator, however, is slightly more complicated — there exists 

a factor of e VN ^ in the expanded expression of e 2vN ^ ( ^r- \ ■ 76 The presence 

of this exponential weight factor in the integrand means that mathematically, 
we need to consider (0(£, t) — 2 (£, t)) t pretty much only at the tip of the front 
(i.e., around £ , where the cutoff is implemented; see Fig. 13), and physically 
it means that the fluctuations at the tip of the front contribute the most to 
Dq. At this region, due to the fact 2 (£,t) <C <f>(£,t), 2 (£,t) can be dropped 
from Eq. (107). Thereafter, with (<f>(£,t)) t = 0^(0) the contribution of the 
tip region to the numerator turns out to be of O(N) at the leading order, 
which then cancels the 1/N prefactor. 

The 1 / In 6 N asymptotic scaling of Dq is finally obtained by noticing that the 
region behind the tip cannot contribute enough to cancel the 1/N prefactor. 

The l/ln 3 iV asymptotic scaling of Df, on the other hand, is a much more 
involved process; but we can already considerably simplify its expression with 
the results of Sec. 4.2.1. For example, Df — D^p+D^+D^p is a simplification. 
We also know that Eqs. (103), (105) and (106) are exact results for the ltd 
differential equation (98). Since and D^f 1 for the sFKPP equation are 
easily seen to scale ~ 1/N from Eqs. (103) and (105), our focus naturally 
concentrates on the asymptotic scaling of D . 

Formally, it is not possible to simplify Eq. (106) any further. We have also 
demonstrated in Sec. 4.2.1 that one cannot also simply replace </>(£, t) by (£) 
everywhere for fluctuating "pulled" fronts — doing so yields the wrong form 
D*f \ namely that D^p = ~ 1/N. Instead, to evaluate , one really 
has to compute the integrals and the sums in Eq. (106). As it turns out, doing 
so is full of subtleties as well. 



We are using, at the leading order, vn — v* = 2 = 2A*. 
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The first difficulty we face is that we cannot calculate J d^ miR (^) eas- 

— oo 

ily. While we do know the solutions of ^ m ,R{0 for the low lying modes 
from Sec. 2.4.2 - - essentially all the fluctuation modes are localized within 
the leading edge of the front — there is a danger in using them to eval- 

oo 

uate J d£^ m ^ R (£) right away. This is easily demonstrated as follows: at 

— oo 

the leading edge, = A m sm[q m (£ - £i)]e~ A * 5 , where q m = (m + 

l)7r/lniV and the normalization A m ~ y^/ ln 1//2 N is chosen in such a way 

oo 

that J dC*m,L(0*m',ii(0 = 8m,m'- When this expression of * m ,i?(0 is used 

— oo 

oo 

for not very high values of m, one gets J d£ty miR (£) ~ (m + l)/ln 3//2 iV at 

— oo 

the leading order. The point to note here however, is that since A* = 1, prac- 
tically the only contribution to the integral comes from the region around £1, 
i.e., at the very left of the leading edge within a distance of 0(1/ A*) = 0(1), 
where the argument of the sin-function approaches zero [this is what allowed 
us to replace sin[g m (£ — £i)] by q m (^ — £i)]. Precisely at this region, the nonlin- 
ear term in f ((/>), too, becomes significant within a distance of (9(1) (see Fig. 
13). Therefore, although in Fig. 13, V(£) can be replaced by Vo(£) to obtain 
r m , ^ m ,i?(0) an d the normalization constant A m with a reasonable accuracy, 

oo 

^m,i?(0 cannot be used to evaluate J d^ m:R (^) so simple-mindedly. 

— oo 

(3) 

To circumvent this difficulty, we rewrite Dj in a different way: 



Df = i lim Jit' L 7< « () £M W , ( + A ) 

1 iT^coJ J J \\ 56 0(O)( O (50 0(0) // 

£ OO OO ' £ 

=£ T .iaMi)), LM ^ K) L, W) s (0 . (108) 

2 J 56 (h J 56 sio) (th 



oo 

In this form, / d£ 



0(0) (0 



^m,n(0 is the contribution of the m-th shape 



fluctuation mode to S [Eq. (101)]. The rationale behind writing Dj- in this 
form is to articulate the fact that for the uniformly translating front solution 



the front speed is given by J d^ /[0^(£)]> and naturally, when the front 

— oo 

shape (f>(x,t) deviates from 6^ by an amount ^m,R,(0i the contribution of 
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the shape fluctuation to the to the fluctuation in the front speed is given by 
8v m f m = / dC, — :— — ^mii(0 at the leading order [just like in Eq. (26) for 

J 0(f) 0(0) (f) ' 

— 00 

the clock model; the notation 5i> m f,m, too, is motivated by Eq. (26)]. 

From the paragraph above Eq. (108), we learnt that to evaluate 5i> m f,m exactly, 
we need the full solution of Eq. (17) [this we do not have]. Nevertheless, what 
we are really interested is how <5i> m f,m scales with IniV. To determine this, 
using A m = y/2/ In 1 / 2 N, we first rewrite the front shape only due to the 
shape fluctuation mode ^ m ,R at the (linearized) leading edge as 



•#>({.«) = 4> m (Q + *™,h(0 

^ e -«{ sinM+V 2,||f}. (109) 

In this form, it becomes immediately clear that for any pure normalized fluc- 
tuation mode, the fluctuations in the front field are always an order 1 / In 3 / 2 N 

weaker than the front field itself. Thereafter, as we further realize that for 

00 

Brunet and Derrida's cutoff solution, J dt; f[<j>^°\^)] = vn — v * = 2 is actu- 

—00 

ally of O(l), we conclude that the presence of the pure normalized m-th shape 
fluctuation mode introduces a fluctuation <5t> m f,ra in the speed of the centre of 
mass of the front, such that 

Sv m{ , m = { J^ N , (HO) 
at the leading order, and s rn of 0(1). Equations (108-110) then imply 



D? = ± E T m S fp£(c m (t)c m ,(t)) t . (Ill) 

Furthermore, since there are 0(\nN) number of bound states within the po- 
tential well Vo(£) of width In and depth (9(1) [see Fig. 13], the sum over m 
and ml in Eq. (Ill) runs from 1 to IniV. 

The last tricky part now is how to obtain the dependence of (c m (t) c m /(t)) t on 
IniV, and how to evaluate the sums in Eq. (111). To this end, notice that for 
a given realization, c m (t) is obtained from Eq. (87), and therein, the presence 
of the factor e A *^ in the ^ m ,L{0 implies that c m (t) for a given realization is 
practically determined from the fluctuation characteristics at the tip of the 
front. Thus, in Eq. (87), we keep the integral over £ only over the leading edge 
of the front and write 
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c m (t) = j-rj^ J ^e A ^sin[g m (£ - 6)] S<f>^,t) . (112) 

Although Eq. (112) yields (c m (t))t = as it should, in the absence of any 
statistics of the shape fluctuations of the front, one cannot obtain an expression 
of (c m (t) c rn i{t)) t from it. Therefore, as far as any formal derivation of the 
1/ln 3 N asymptotic scaling of D^p (or Df) is concerned, one simply cannot 
progress beyond Eq. (112). 

Nevertheless, we can still proceed with two approximations. The first one 
stems from the fact that although it is clear from Eq. (112) that c m (t) and 
c m i(t) (m 7^ m') are correlated in general [after all, for a given realization, all 
the c m (t)'s are determined through the same 5(f>(t)], these fluctuation modes 
will have a finite correlation "length", i.e., (c m (t) c m r(t)) t will be negligibly 
small [compared to both (c^(t)) t and (c^/(t))^ values] when \m — m'\ exceeds 
a certain threshold a <C IniV. Based on this observation, our approximation 
is to choose a = for the extreme (and unrealistic) case to simplify the 
expression for Df to ( see later for the discussion on non-zero values of a) 



1 IniV 2 

Df ^E^W' (H3) 



)(3) _ 1 ^ S m I 9 



Then the second approximation is that due to the presence of the e A ^ in the 
integrand of Eq. (112), only the value of 5(f) within a distance ~ 1/A* = 1 
of the tip determines c m (t). This is argued in the following manner: typically 
the magnitude of 50(£, t) is of order \J(j)(£, t)/N ; at the tip, e A *^° ~ iV cancels 
5(f)(£,t) ~ 1/N, but further behind, the l/^/N factor of yJ(j>(£,t)/N can no 
longer be compensated by e A * 5 . We therefore use 

|Cm(t)l ~ h^ J d ^ Sm ^~^-]^ = li^T" ( 114 ) 

With T m = In 2 N/[ir 2 {(m + l) 2 — 1}], the asymptotic scaling of Df is finally 
obtained as 



D f ~Df ~±y „ V "^ V ?\ ^— ■ (H5) 



1 ^ (m + l) 2 5 2 
' ~ 2 ^ [(m+l) 2 -l]ln 4 iV In 3 iV 



We end Sec. 4.2.2 with four observations: (i) with more reasonable assumption 
a ^ 0, the expression for D^p certainly gets more complicated. However, so 
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long as a C In TV, which is what one expects in reality, the l/ln 3 iV asymp- 
totic scaling of Df (and hence of Df too) still continues to hold, (ii) We 
have extensively used the construction for the left eigenvector of the stability 
operator C VN for Fisher equation (with a growth cutoff below <fi = 1/N) all 
along Sec. 4.2.2. For clock model [122-126], or for the model that Brunet and 
Derrida considered [21] to demonstrate the 1/ In 3 N asymptotic scaling of Df, 
construction of the left eigenvector for the corresponding stability operator is 
not easy. In view of that I do not know how to repeat the same derivation 
for those two models with an equivalent amount of rigour. Nevertheless, since 
all the arguments for the sFKPP equation between Eqs. (109) and (115) are 
concentrated on the leading edge or at the very tip of the front, they can be 
easily repeated for these two models to derive the same 1/ In 3 N asymptotic 
scaling of Df. (iii) In the clock model, in terms of 0, one cannot create a lo- 
calized fluctuation in the front shape — any fluctuation in the front shape is 
by definition non-local. The "collisions" between clocks are also non-local in 
nature [as reflected in Eq. (25)]. From these complications, although a priori 
it may seem that the clock model may have a different asymptotic scaling of 
Df, in reality, it does not happen. In fact, Brunet and Derrida showed in a 
simplified version of their original microscopic model (which closely resembles 
the clock model) that the introduction of localized fluctuations in at the 
very tip of the front is all that is needed for the 1/ In 3 N asymptotic scaling of 
Df. In that sense [also in view of point (ii)] the 1/ In 3 N asymptotic scaling of 
Df seems to be a generic property of fluctuating "pulled" fronts, independent 
of the microscopic model, (iv) Last but not the least, through the scaling of 
the c m (t)'s, the l/ln 3 iV asymptotic scaling of Df is seen to be determined 
only from the tip of the front. This is in perfect agreement with Brunet and 
Derrida's "simplified model" [21]. 

4.3 Epilogue III 

The existence of two different expressions for front diffusion, Df and Dq seems 
confusing to say the least — after all, any front realization maintains its in- 
tegrity at all times, and thus, at no instant of time, the Goldstone mode and 
the centre of mass of any given front realization can be infinitely far apart 
from each other. This implies that whether measured through the Goldstone 
mode, the center of mass, or the any place where the value of the front field 
reaches any fixed value, say <f> , operationally (i.e., in a computer simulation) 
one should always obtain the same front diffusion coefficient! 

In fact, the present situation can be compared to the a similar one that occurs 
in gas mixtures (see for example, Chap. 11.2 of Ref. [46]) Therein, the expres- 
sion (and the value) of the the diffusion coefficient depends on its precise def- 
inition (although operationally there is exactly one diffusion coefficient), but 
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these different expressions of the diffusion coefficient are (quite non-trivially) 
related by means of the Onsager relations for the diffusion coefficients. As for 
fronts too, it should not be surprising that the precise values of Df and D G 
are not the same — conceptually they are simply two entirely different quan- 
tities. Whether they could be related by any clever means or not is left here 
for future investigation. 

It is however not enough to stop at the above paragraph; one has to identify 
what physical quantities Df and D G respectively correspond to. Although 
surprising it may seem in the first place, I will argue in the next few paragraphs 
that for fronts with 5(f) — > at £ — > ±00 and finite r m 's, Df is the long-time 
front diffusion coefficient, while Dq is the short-time (or the instantaneous) 
diffusion coefficient measured through the Goldstone mode of the front. 

The most important point to note is that from its very definition, Dq can be 
measured by tracking the spread in the position of the Goldstone modes for 
an ensemble of front realizations over an infinitesimally short interval of time; 
one does not at all need to follow the front realizations over a long period of 
time. Through its very definition, therefore, Dq completely ignores the time- 
correlations that exist in the speed fluctuations for individual front realizations 
of the ensemble (precisely these correlations exhibit themselves through the 
relaxation time scale r m of the m-th mode of front shape fluctuation). These 
correlations do affect the whole dynamics of the individual front realizations, 
and only Df captures the full flavour of it (c.f. Sec. 4.2), albeit one must pay 
a price of having to deal with more complicated algebra. Thus, it should not 
come as a surprise that Df in general yields an expression different from Dq. 
The best example that supports this scenario — rather an extreme case where 
r m 's can be arbitrarily large — is the fluctuating "pulled" front (Sec. 4.2.2). 
Therein, Dq simply yield a wrong scaling for the front diffusion as iV — > 00. 

In view of existing simulation results on front diffusion, it is worthwhile to re- 
visit the examples for front diffusion in the literature (reviewed in this article) 
with a critical eye to clarify the above viewpoint. 

1. When R((/),£,t) is replaced by R(<f>(°\Z,t) in Sec. 3.1.2: with the 0(e 1/2 ) 
neglected in Eq. (74), in this case, there is no difference between Eqs. (72) 
and (98). Then, with D — 1, one obtains 



he 



Df = { J dZg^]^D G = £ 



2 U(0)1 



-00 

The first point to note for this case is that just like D G , Df also scales as 



(116) 
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e. The second one is that for pushed fronts, the numerical values of Df and 
D G , in all likelihood, are close by (i.e., one needs good data to distinguish 
between the two). This is in fact easily seen from the following argument: 
pushed fronts are strongly localized, and therefore, their Goldstone modes 
$g,r are localized in space, say around £ = £ - Having used this information 
in Eq. (116), we see that for Dq, the e Ve ^ factors in the numerator and the 
denominator cancel, and thus, effectively apart from some numerical coef- 
ficients of 0(1), D G ~ ^ [0 (o) (£o)]/<^(£o), while D f ~ heg 2 [^\^)}/2. 
How their numerical values compare depends on $g R (£o)- 

To me therefore, the issue regarding the difference between Df and D G 
for pushed fronts requires further numerical investigation. One particular 
instance where the above argument breaks down is when the actual front 
region is very large in size. This is precisely what happens for fluctuating 
"pulled" (or weakly pushed fronts), i.e., when <p^> is the Brunet and Der- 
rida's cutoff solution (11). In that case, there is a clear difference between 
Df and D G through their respective scalings (also see 3 below). 

2. Fluctuating pulled front with multiplicative noise in Sec. 3.3: in this case, 
it has been clearly shown, numerically and as well as via the Cole-Hopf 
transformation, that the front wandering is subdiffusive. The first sign of 
the fact that there is a problem associated with the separation of time scales 
in this particular example (slow for front wandering and fast for front shape 
fluctuations) is manifested through the breakdown of convergence of the 
integrals in the definition of D G . A priori, it seems strange that Df [see Eq. 
(116)] does not suffer from any similar peculiarities. 

The reason for the expression of Df not exhibiting peculiarities, however, 
is not hard to trace. It is in fact hidden within the very derivation of Eq. 
(116) — recall that there have been cancellations of factors of r m all along 
Sec. 4.2.1. For a pulled front r m 's are infinite, and these cancellations are 
meaningless. The derivation of Eq. (116), thus, simply breaks down for 
fluctuating pulled fronts with multiplicative noise. 

Nevertheless, since W(t) has been measured operationally via the centre 
of mass in Ref. [98], it makes sense to use Df to calculate W(t) as well. 
With the introduction of the upper cutoff £ c ~ \/4Dt, which makes the r m 's 
finite, one can very easily show [having R((j), £, t) replaced by R(<f>(°\£,t)] 
that 

4:iD 5/2 3e 

for t 3> 1. Interestingly enough, in this example, one gets away with the 
usage of D G to obtain the theoretical curve (dashed line) for W(t) in Fig. 
32 — both Df and D G scales as 1/Vt for t 3> 1, and for Fig. 32 that is all 
that matters. 

3. Fluctuating "pulled" fronts, such as the sFKPP equation: although we 
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have discussed this at length in Sec. 4.2.2, it is still useful to articulately 
point out a few distinct subtleties. First, if the noise term in Eq. (96) is made 
additive by replacing by where (f)^ is Brunet and Derrida's cutoff 
solution (11), Df scales as 1/N [c.f. Eq. (116)] for N — > oo, as opposed 
to l/ln 3 iV. It is therefore important to take notice of this subtlety for 
simulation purposes — Brunet and Derrida properly take this into account 
in their own simulations [21]. For Dq however, it does not matter whether 
one uses or </>(°) for the noise term in Eq. (96) [this is quite obvious from the 
expressions of D G in Eqs. (107) and (116)], one always gets D G ~ l/ln 6 iV 
for N — > oo. 

Secondly, since Dq is the short-time (or the instantaneous) diffusion co- 
efficient of the front, if the diffusive spread (Ax G (t)) in the front posi- 
tions measured through the Goldstone modes are plotted as a function of 
time for an ensemble of fluctuating "pulled" fronts, one would see that 
d(Ax G (t))/dt = l/ln 6 iV right from the very start. On the other hand, 
when the front positions are tracked by the positions of their centres of 
mass, one expects d(Ax'j(t)) / dt to be smaller for small t, but it should nev- 
ertheless converge to 1/ In 3 iV for t — > oo. In fact, one can further argue that 
d{Ax 2 f (t))/dt\ t ^o ~ l/ln 4 iV. In order to do so, one simply has to observe 
that the l/ln 3 iV scaling of Df is obtained from IniV number of modes, 
each of which contribute an amount of l/ln 4 iV to Df [c.f. Eq. (115)]. The 
shortest and the longest surviving modes have lifetimes respectively of 0(1) 
and 0(ln 2 N). Thus, d(Ax 2 f(t)) / dt\ t ^ captures only the contribution of the 
shortest surviving modes and should behave as ~ 1/ In 4 N, while in order 
to observe the true l/ln 3 iV scaling of Df, one has to wait till 0(ln 2 iV) 
time for the contributions of all the IniV number of modes to Df to add up 
together. 

I finally end Sec. 4 (and this review article too) with the comment that further 
simulations and research are necessary to confirm the conjectures made in Sec. 
4.3. 



5 Author's Note 

This review article has been written from my own perspective, and it is clearly 
reflected by the organization of the topics that I have covered. After a brief 
introduction to the subject, I spent a longish section on the discrete particle 
and lattice models of fluctuating fronts. Needless to say, there is no unify- 
ing underlying structure for the plethora of results in that section — each 
problem or topic has been dealt with a different approach in the existing lit- 
erature. The two sections 3 and 4 I spent on field-theory — therein, we find 
an organized theoretical backbone for studying the effects of fluctuations on 



105 



propagating fronts by means of stochastic differential equations. Finally at the 
end, I tried to merge the predictions of field-theories on internal fluctuations 
and the corresponding results for the discrete particle and lattice models. 

Unless existing derivations are too involved or too complex, 1 have tried to 
outline the main line of procedure. The only exception has been the issues of 
front diffusion in one spatial dimension, specially for the fluctuating "pulled" 
fronts, where I have delved deep in the dirty world of algebra. This was done 
with keeping in mind that front diffusion in one spatial dimension is a relatively 
new topic. In any case, I have also thrived to put in physical insights wherever 
appropriate. At various places, open issues have been spelled out for further 
investigations as well. As for myself, it would be gratifying to see if the expert 
researchers, and as well as the beginners in this field, both of whom are the 
target readers, find this review article useful. 
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